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ADMINISTRATIVE COMMITTEE MEETING 


The following is abstracted from the minutes of the 
Administrative Committee meeting held at Jansky and 
Bailey, Washington, D. C., January 11, 1957. Present 
were J. I. Bohnert, Vice-Chairman, who presided, P. S. 


- Carter, Arthur Dorne, Coleman Goatley (guest), F. T. 


Haddock, J. W. Herbstreet, D. D. King, R. L. Matting- 


: ly, and D. C. Ports. 


1) Mr. Mattingly presented a brief membership and 


financial report. Membership has grown each year. As of 
_ the end of 1956, our membership totaled 2392, of which 
365 were student members; whereas, on March 8, 1956, 


total membership was only 1891. Also, as of the end of 
1956, our treasury balance was $3448.00, a deficit hav- 


; ing been avoided through the postponement of TRANS- 


ACTIONS publications. 
2) There was an extended discussion of programs for 


_ improving our financial status to ensure that in the fu- 
ture it will not again impede either the publication of 


our TRANSACTIONS or any other phase of our profession- 
al activity. It was concluded that the publication of ad- 


_ yertising in the TRANSACTIONS (as just recently permit- 
_ ted by IRE policy) constitutes the most promising ave- 


nue for use, and further that our goal should be the sale 
only of sufficient advertising to meet our actual needs. 
In view of this, it was decided that we should not engage 
a commercial sales representative, but instead should 
handle the sales ourselves. Mr. Ports volunteered to ad- 


- minister this activity. His offer was gratefully accepted 


by the Committee. He will, with the assistance of other 
members of the Administrative Committee, organize a 
program of personal letters and calls to likely prospects. 
It was the belief of the Committee that the TRANSAC- 


trons should be a uniquely valuable medium for many 


‘advertisers, such as manufacturers of test equipment 


_ which meets our particular needs and laboratories seek- 


ing personnel in our profession. 


3) The next topic was the Affiliate Membership Plan 
recently approved by IRE headquarters. In this plan 
the membership of certain selected scientific societies 
would be extended the privilege of joining the various 
professional groups on an “affiliate” basis without join- 
ing the IRE. Such a member of PGAP would incur less 
than half the cost of full IRE-PGAP membership and 
would still receive the TRANSACTIONS. The IRE Execu- 
tive Committee must approve any such societies select- 
ed by the professional groups. After discussion, during 
which it was brought out that the committee does not 
yet have a full understanding of the intent of this policy, 
it was decided to select a few societies, and submit them 
for approval, together with a request for a clarification 
of the type of society considered appropriate, and an in- 
dication of whether the total number was expected to be 
limited. The following were then selected: The American 
Astronomical Society, the American Meteorological So- 
ciety, the Optical Society of America, and the Royal 
Astronomical Society. It was noted that the complete 
list of societies which might possibly be included is a 
very extensive one, numbering perhaps in the hundreds. 
In the case of some of these there is a very substantial 
overlapping of interests with IRE. 

4) Mr. Goatley reported on the Symposium on Optics 
and Microwave Techniques held in November at the 
George Washington University, pointing out in particu- 
lar that for its success we are deeply indebted to the 
University, which cooperated and assisted most com- 
pletely. 

5) Mr. Ports reported that he has been contacted by 
John P. Hagen in connection with the URSI Spring 
meeting. The Administrative Committee passed a reso- 
lution indicating our approval of the PGAP’s joint spon- 
sorship with the URSI of this meeting. The Secretary 
will advise Dr. Hagen of this action and Dr. J. T. Boll- 
jahn (PGAP Meetings Chairman) will be instructed to 
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implement this decision. The Secretary will also send a 
letter to the Executive Committee of the IRE asking 
for their approval. It is hoped and expected that the ab- 
stracts of this Spring meeting will be published in the 
TRANSACTIONS. 

6) Mr. Dorne brought up the matter of how News and 
Views could best further the development of PGAP and 
suggested a number of alternative types of material 
which it might include. After discussion it was agreed 
that News and Views should concentrate upon reporting 
specific items of information which best convey a picture 
of the current activity and concerns of PGAP and its 
members, with the hope that this will increase the co- 
hesiveness of the organization and stimulate more active 
and widespread participation in its affairs by the mem- 
bership. Mr. Dorne also suggested that the society might 
benefit along these lines if it sponsored social events 
which would acquaint the members with each other. It 
was agreed to hold an informal gathering for cocktails 
sometime during the National Convention. 

7) Mr. Ports brought up the special IRE issue dedi- 
cated to Radio Astronomy. C. M. Jansky was originally 
scheduled to handle the task of procuring and editing 
the material for this issue. Unexpected commitments 
have made this impossible and Frederick Haddock is 
now in charge. Dr. Haddock pointed out that there are 
many ways of handling an issue of this type depending 
upon the type of reader to which it is most particularly 
directed. The papers might, for example, emphasize 
either the techniques of radio astronomy or alternatively 
astronomy itself. He asked for an expression from the 
committee of the direction he should take. It was the 
consensus that the issue should be directed toward non- 
specialist electronic people and that it should be bal- 
anced to include papers on both the electronic tech- 
niques and astronomy, of the type which will give the 
best over-all picture of the field, despite the fact that 
this would necessarily exclude much of the type of new 
detailed material which would be of most interest to in- 
dividuals working in radio astronomy. 

R. L. MATTINGLY 
Secretary-Treasurer 


OTHER ADMINISTRATIVE NEws 


Delmer Ports reports that initial response to the ad- 
vertising sales campaign is favorable. To date, as a re- 
sult of Committee activity, five tentative commitments 
have been received for the July issue. Note: The dead- 
line for copy for this issue is June 1. All PGAP members 
are urged to help in this activity. Prospective advertis- 
ers should be referred to: 


Delmer C. Ports 

c/o Jansky & Bailey 

1339 Wisconsin Avenue, N.W. 
Washington 7, D. C. 
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The Radio Astronomy Issue of the PROCEEDINGS is 
now definitely scheduled for December of this year. 

S. M. King is the PGAP representative on the Na- 
tional Convention Record Committee. 


FutTurRE MEETINGS AND SEMINARS 


The 1957 WESCON will be held in San Francisco's 
Cow Palace on August 20, 21, 22, and 23. Four sessions 
of the technical program are expected to be devoted to 
antennas and propagation. Your assistance in making 
these PGAP-sponsored sessions of the highest possible 
caliber will be greatly appreciated. 

Please submit 200-word abstracts of proposed papers 
before May 1 to: Dr. J. T. Bolljahn, Stanford Research 
Institute, Engineering Division, Menlo Park, Calif. (An 
tenna Papers), or Dr. L. A. Manning, Stanford Elec- 
tronics Laboratory, Stanford University, Stanford, 
Calif. (Propagation Papers). 

A two-week summer statistical seminar will be held at 
the Endicott House in Dedham, Mass., beginning July 
29, 1957. It is understood that the content of this course 
will make it of interest to propagation people. Further 
information can be obtained from Dr. M. E. Terry, Bell 
Telephone Laboratories, Murray Hill, N. J., or the sec- 
retary, Dr. Geoffrey Beall, Gillette Safety Razor Com- 
pany, Boston, Mass. 

Ohio State University has again been invited to con- 
duct a classified Radome Symposium for the benefit of 
the Armed Forces. The meetings are scheduled for June 
3 through June 5 in the Ohio Union at Ohio State Uni- 
versity. 


CHAPTER NEWS 
Akron 


We are pleased to announce that the formation of the 
Akron Chapter was approved by IRE, effective Febru- 
ary 6, 1957, We have not received word at this time of 
the officers appointed. 


Syracuse 


We are also pleased to announce the formation of the 
Syracuse Chapter, which was approved effective Febru- 
ary 6, 1957. We have not been advised to date of the of- 
ficers appointed. 


Washington 


A meeting of the Washington Chapter was held Feb- 
ruary 14, in Tompkins Hall, the new building for the 
School of Engineering at George Washington Univer- 
sity. John E. Jackson, of the Naval Research Labora- 
tory, presented a paper on the latest results of Rocket 


Studies of lonospheric Propagation in the 4 to 50-mega- 
cycle range. 


1957 
Personnel Notes 


Robert Elliott, formerly head of antenna research at 
Hughes’ Aircraft, has joined with Robert Krausz, for- 
merly head of engineering research at SRI, Southern 
California Laboratories, to form the RANTEC Corp. 
in Calabassas, Calif. Other members of the new com- 
pany are Alvin Clavin, formerly with Canoga Corp., 
Lou Kurtz, formerly with Hughes, and Bill Harrison, 
formerly with SRI. 

Recently-formed Systems Laboratories Corp., Sher- 
-man Oaks, Calif., a group of some 36 scientists and en- 
gineers under President John Barnes has announced the 
election to Vice-President of James Marsh, formerly 


News and Views 
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head of the Autonetics Radar Group. Also at SLC in the 
antenna and radar field is Don Margerum, formerly 
with SRI. 

George Underberger, formerly with Bendix-Pacific, 
has joined Kearfott’s microwave group in Van Nuys, 
Calif., which is headed by W. A. Hughes, formerly with 
Hughes Aircraft. 

Bob Blaine has formed Blaine Electronics, 14757 Kes- 
wick, Van Nuys, Calif., a company devoted to fabrica- 
tion of models and electronic equipment for antenna 
measurements. 

L. Cutrona is on a leave of absence from the Univer- 
sity of Michigan while heading an antenna group at 
Ramo-Wooldridge Corp. 
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contributions 


Some Electromagnetic Transmission and Reflection 
Properties of a Strip Grating* 
ROBIN I. PRIMICH} 


Summary—Some theoretical and experimental results for the 
reflection and transmission of a uniform plane electromagnetic wave, 
normally incident on an ideal strip grating, are presented. The theory 
is based on the variational method, and the measurements were 
made at normal incidence in a parallel plate region operating in the 
8- to 10-cms wavelength range. A few results for oblique incidence 
are also included; the experiments in this case were carried out at 
1.24-cm, wavelength in a parallel plate spectrometer. Resasonable 
agreement between theory and experiment was obtained throughout. 


INTRODUCTION 
\ STRIP GRATING of the form shown in Fig. 1 


is considered in this paper. The strips are as- 

sumed to consist of perfectly conducting, in- 
finitely thin metal and are of uniform width, ¢, and 
infinite in extent. Each grating contains an infinite 
number of equally spaced, identical strips. The problem 
is of interest for the following reasons. 

1) A number of such gratings in sequence will form 
an artificial medium of refractive index less than unity. 
However, it has no advantages over the rodded medium! 
and it would be inferior mechanically. Theoretically, 
this medium illustrates very clearly one of the factors 
which gives rise to the interface phenomena in artificial 
dielectrics.® 

Consider Fig. 2 in which 6 and ¢ are constant and 
do. (Ao = free space wavelength.) 

Suppose that, initially, d is sufficiently large for no 
appreciable coupling to exist between adjacent gratings. 


* Manuscript received by the PGAP, September 1, 1956. 

t Radio Physics Lab., Defence Res. Board, Ottawa, Canada. 

1J. Brown, “Artificial dielectrics having refractive indices less 
than unity,” Proc. IEE, vol. 100, part 4, pp. 51-62; 1953. 
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Then each grating may be replaced by an equivalent 
shunt reactance? and the behavior of the medium may 
be determined to a high degree of accuracy by using the 


> G. G. Macfarlance, “Quasi-stationary field theory and its appli- 
cation to diaphragms and junctions on transmission lines and wave 
guide,” Proc. IEE, vol. 93, part 3a, pp. 703-719; 1946, 
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transmission line analogy. However, for d of the order 
of Xo/ 10 or less, coupling between gratings can no longer 
be ignored and the analogy breaks down. The limiting 
case, d=0 (or db) is that of a system of parallel 
metallic plates of finite thickness, which possess the 
interface effect to a marked degree; 7.e., in particular 
the reflection coefficient at the interface (s=0) of the 
medium no longer bears a simple relationship to the re- 
fractive index of the medium. It is evident that interface 
effects will be present in all similar medium in which 
d<KXo <0. . 

2) A-single grating can be used for the determination 
of wavelength and will act also as a polarization filter. 
A study of the higher-order waves diffracted by the 
grating serves as an accurate measure of the exciting 
_ wavelength. This is a standard optical procedure, but is 
now of interest in this millimeter wave region.‘ An inter- 
_ esting application is one in which a grating is used to 
separate out a particular magnetron harmonic.’ 

By restricting the spacing between adjacent strips, 
the grating is made opaque to waves whose polarization 
is parallel to the strips. This enables a linearly-polarized 
wave from an arbitrarily-polarized source to be 
produced. 

3) A strip grating has similar properties to the inter- 
face of a semi-infinite system of metallic plates of finite 
thickness.® This is of interest in microwave lens applica- 
tions so that the possibility of a similar behavior in re- 
gard to the surface reflecting properties arises. In fact, 
an empirical relation involving the grating reactance has 
been used to determine the effect of plate thickness in 
a metal plate system.’ If a rigorous basis for this analogy 
could be found it would be extremely useful, but a brief 
consideration of this matter, given at a later stage, 
indicates that if this analogy does exist, it would be in 
a form too complicated to be of value for parallel plate 
lens design. 

Basic to the above applications is the requirements 
_ that the properties of a single grating be known. There- 
- fore the following sections will review grating theory as 
applicable to metal strips and it will include the applica- 
tion of the variational method which was used to calcu- 
- late the reflection coefficient of an ideal strip grating. 
Experimental results will be included. 


THEORY OF STRIP GRATINGS 
Review 


In the field of optics, attention has been mostly con- 
fined to the angular distribution of waves arising from 


8 J. Brown, “The design of metallic delay dielectrics,” Proc. IEE, 
vol. 97, part 3a, pp. 45-48; January, 1950 ; 

4J. L. Farrands, “The Dielectric Properties of Some Ketones 
at Cm and Mm Wavelengths,” PhD. dissertation, University of 
London, England; 1953. : 2 

5 J. A. Klein, et al., “Magnetron harmonics at mm wavelengths, 
Rev. Sci. Instr., vol. 23, pp. 78-82; February, 1952. ; 

6 R. I. Primich, “A Study of the Properties of Artificial Dielec- 
trics,” Ph.D. dissertation, University of London, England; 1954. 

7D. J. Epstein, “Phase Shift of Microwaves in Passage through 
Parallel Plate Arrays,” M.I.T., Lab. for Insulation Res., Tech. Rep. 
no. 42; August, 1950. 
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the grating. These are determined by the periodicity of 
the grating. The detailed calculation of the scattered 
fields has been virtually ignored, mainly because of the 
role of the reactive fields in the vicinity of the grating 
has not been fully understood (see Wood’s anomalies’). 
Palmer has been reexamining these so-called anomalies 
recently.® 

In the microwave field, the complex nature of the re- 
flection caused by an arbitrary obstacle is of major in- 
terest and it is necessary to solve the complete field 
problem. 

The first attempt to deal with the strip grating was 
made by Macfarlane.” It should be noted at this point 
that by the use of Babinet’s principle!® the solution of 
the problem indicated in Fig. 1, for any angle of inci- 
dence, is also the solution of the complementary grating 
in which the metal and aperture regions are inter- 
changed and the polarization is rotated through 90°. 
Macfarlane, among others, pointed this out and pro- 
ceeded to solve, using the approximate quasistatic 
method, the #-plane (or capacitive iris in a waveguide) 
and H-plane gratings in free space. In a separate paper?! 
dealing with a grid of circular wires, he gave the condi- 
tions for the existence of higher-order propagating 
waves. 

The “Waveguide Handbook,” using an approximate 
treatment based on an integral method, gives numerical 
results for 6<Xpo (see Fig. 1). In both of the above works 
the range b>Xp is excluded. 

Booker’® has given a more general qualitative dis- 
cussion. 

Miles appears to have been the first to consider all 
values of 6/Xo, but only for normal incidence. This is an 
extension of his earlier work. He obtains the reflection 
and transmission coefficients for the strip grating from 
a rapidly converging expansion of the integral equation. 
The resonance in either coefficient at the point b=NXo, is 
specifically mentioned by him for the first time, but 
surprisingly no recognition is given to the fact that 
higher-order propagating waves are present and, in 
particular, that the resonance is due to one of the 
evanescent waves starting to propagate freely. This does 
not effect the validity of the results for the main propa- 
gating waves, but as he uses the variational method to 
check the accuracy of the results, it is dubious whether 
the limiting procedure described has much significance 


8 R. W. Wood, “Physical Optics,” The Macmillan Co., New York, 
NepyelO1i: 

9C. H. Palmer, “Parallel diffraction grating anomalies,” J. Opt. 
Soc. Amer., vol. 42, pp. 269-276; April, 1952. ; 

10H. G. Booker, “The elements of wave propagation using the 
impedance concept.” J. IEE, vol. 94, part 3, pp. 171-202; 1947. 

1G. G. Macfarlance, “Surface impedance of an infinite parallel 
wire grid at oblique angles of incidence,” Proc. IEE, vol. 93, part 3a, 
pp. 1523-1527; 1946. ; 

122.N. Marcuvitz, “Waveguide Handbook,” McGraw-Hill Book 
Co., Inc., New York, N. Y.; 1951. 

13 J, W. Miles, “The diffraction of a plane wave through a grat- 
ing,” Quart. Appl. Math., vol. 7, pp. 45-64; January, 1949. 

4 J. W. Miles, “The analysis of plane discontinuities in cylindrical 
tubes,” J. Acoust. Soc. Amer., vol. 17, pp. 250-285; 1945. 
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in the presence of higher-order propagation. 

A. E. Heins (in a personal communication) has used 
the Weiner-Hopf technique! to obtain a rigorous solu- 
tion for the special case of the half-open grating and the 
significance of the higher-order waves is fully realized. 

Recent work by Muller!® seems to have the greatest 
generality combined with actual numerical computa- 
tions. In treating an arbitrary strip grating for any 
angle of incidence, the usual Fourier expansion for the 
fields is transformed to a pair of Fredholm integral equa- 
tions of the first order, 7.e., 


+B ar 


Ais 


P(e — x')a(x')dx’ B= Bue a= 
—Br b 


b 


+ar ‘ 
ropagation constant 

B — ee daal , / =, p So 
dix EVIE IEL |e of mth order wave. 


a 
P(x — x’) =— >> —exp(jmx— x’). 

2m m= JY¥m 
u(x’) is related simply to the transverse electric field in 
the aperture for the E-plane case, v(x) is related to the 
strip current in the H-plane case, and A and B are 
simply related to the corresponding field amplitudes. By 
a simple shift of the coordinate system two further 
equations are obtained: 


1 aa 

Aa I) Var —{ x:u'(x)dx  E-plane case. 
Te arg 
1 +B 

B= —j/y+— x-v'(x)dx  H-plane case. 
2a —Br 


u'(x) is now associated with the strip current and 
v(x) with the aperture electric field. An approximate 
solution (which was promised for a future publication) 
is valid for (b/No) <4; when [(b—a)/b| <3 the first pair 
of equations are used, while for (a/b) <4 the second pair 
are used. It was demonstrated that the integral equa- 
tions are encountered in the study of capacitive and in- 
ductive irises in rectangular waveguides!® and that the 
numerical solutions are equally applicable. Muller has 
then given a general practical solution although the 
question of higher-order wave propagation is not spe- 
cifically discussed. 

Twersky" has obtained a general rigorous solution for 
an arbitrary grating, but numerical computation has 
proved to be too difficult except in some special cases. 


The Variational Method 


Miles used the variational method to check the ac- 
curacy of his results, but apart from the general formula- 


15 J. F. Carlson and A. E. Heins, “The reflections of an electro- 
magnetic plane wave by an infinite set of plates,” Quart. Appl. 
Math., vol. 4, pp. 313-329; January, 1947. 

1 R. Muller, “Eine strenge behandlung der beugung elecktro- 
magnetischer wellen am streifengitter,” Zeits. Naturforschung, band 
8a, pp. 56-60; January, 1953. 

17V. Twersky, “Multiple scattering of radiation by an arbitrary 
configuration of parallel cylinders,” J. Acoust. Soc. Amer., vol. 24, 
pp. 42-46; January, 1952. 
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tion no details are given. The alternative development 
to be given here is believed to be more convenient for 
computation and moreover leads to summations which 
are identical to those encountered in a previous prob- 
lem. 

An ideal strip grating will be considered. This consists 
of an infinite number of similar, infinitely thin, perfectly 
conducting, metallic strips each infinitely long. The 
strips are equally spaced and have the same width. A 
uniform plane wave, incident normally to the plane of 
the grating is assumed. (See Fig. 1.) 

Because of the periodicity of the grating plane, it fol- 
lows that the free space fields may be written** for 
z=0: 


E, = exp (—jBoz) + Ao exp (JGoz) 


2 Q2mrx 
+ >) Am cos exp (jamz). 
m=1 


OF ‘ ; 
ns = — 7B exp (—JBoz) + 78040 exp (7602) 
Z : 
a 2mm x . 
+i > AmAm COS exp (jamz). 
m=1 
2=0: 
S 2mrx : 
E, = Co exp (—jBoz) + >) Cm cos exp (—jam2). 
m=1 
OE . . 
— = — jBoBo exp (—jBoz) 
Oz 
cs 21x 
— 7 > amCm COS exp (—jam2). 


m=1 


where am is the propagation constant of the mth order 
wave in free space, a») =» and A,, and C,, are the ampli- 
tudes of the mth order waves on either side of the grat- 
ing. The time variation exp (jw#) is assumed throughout. 

Using the boundary conditions in the grating plane 
(2=0), 


cs "T2m3rx 
E, = 1+ Aot+ >) Amcos 
m=1 
5 |x| <a. 
& marx 
= Co + 2 Ca cos 
m=1 
=O | a | < (b — a), 
so 
1 
1 -- Ay ca Co = —{ Ee ide: 
bidt ass 
2 2mm x 
Am = Cn = —{ Ey, COS dx. 
Dw Gp 


18F, Berz, “Reflection and refraction of microwaves at a set of 
parallel metallic plates,” Proc. IEE, vol. 98, part 3, pp. 47-55; 
January, 1951. 


195% 


Equating the magnetic fields in the aperture regions, 


i ae 2 
—JBo(1 — Ao) +7 DS amAm Cos sl 
m=1 
7 rac Qmarx 
=4f8iCo— 4 y, On m COS 
m=1 
or 
Bo(1 — Ao) = Bo(1 + Ap) 
4 = 2mrx 2 
+ —= D7 am COs if E,, cos ges dx. 
De at b ae 


Multiplying by £,,, integrating over period “a” and 
using 


1 
1+ 4o=— fo Eqdr. 
b ap 


2m wx 2 
dx 


matty 


Ey, €0s 


p b 
a Ey, &x 
ap 


A similar expression can be obtained in terms of the 
magnetic field in the grating plane, but for present pur- 
poses the above will be sufficient. For computational 
purposes it is convenient to consider two regions of in- 
terest, depending on the value of b/Xo. 

1) (b/Xo) <1. 

If am/Bo= —jQm, Qm will be real and positive. 

Then, if Yin is the input admittance at the grating 
plane, 


po At, c : 
= 4 se 
md, Pot 4 Dia 


Bo- 


PRE iB) 
1+ A 
and 
. Eeyg COS vt x 
Boat) 0a ||. (2) 


m=1 


af fi 0% 
ap 


where B, is the susceptance of the grating. 

This expression for B, can be shown to be stationary 
and consequently any error in an assumed trial field 
E,, will result in only second-order errors in B,. Further- 
more B, has an absolute minimum (for the correct E,,) 
and any error in the assumed £,, will result in a B, that 
is always too large. 

From an expression for the aperture magnetic field 
(which is derived in a manner similar to that above), a 
value of B, which is too small can be found and conse- 
quently it is possible to limit the error in B,. 

The susceptance of the grating is then known within 
definite limits. 

2) (b/Xo) > 1. 

For any given b/A, a certain number p of the a, will 
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now be real. Using a trial field of the form (which satis- 
fies the boundary conditions) 


nr 
Bjlist S21) Gog et |x| <e 
n=1 
in 
f e 2m1rx 2 
cos x 
1— A, Cs) On, Yo 

Vin/= c= Ue et ip aS (GHEE S 4 
1+ Ao x 0 


will result in complex D,’s. This means that the grating 
can no longer be represented by a pure shunt suscept- 
ance but has, in addition, a resistive component of im- 
pedance. The expression for Y, is still stationary® and 
the convergence of Y;, towards the real value can now 
be studied best by the methods discussed elsewhere,® 
1.€., convergence towards a single point in the complex 
plane. 


Computations 


The assumed trial field used for numerical work is 


Tx Tx rx 11% 
Ey, = cos Psi yes Rogier iste C08 = eae CPs 


| «| Saas 


Insertion of this in (1) and use of the stationary condi- 
tion 


OV in 
=) n= 1.3, 5, and 7 
OD, 


results in explicit expressions for the grating admittance. 

1) (6/0) <1. 

The “Waveguide Handbook” predicts that at the 
point b=Xo, the magnitude of the wave reflected by the 
grating is zero whereas the variational method showed 
that the reflection coefficient at this point depended on 
a ratio a/b. Consequently, the formulas from the former 
source should not be used in the vicinity ’=No. Experi- 
mental results to be given later will give some idea of the 
range of validity of the “Waveguide Handbook” 
results. 

The variational expressions were used to compute the 
magnitude of the reflection coefficient for two gratings, 
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for (b/Xo) <1. Experimental work, to be described later, 
shows that the results are very satisfactory. 

2) (b/Xo) > 1. 

In principle the variational method is still applicable 
for this region but the computations are laborious. When 
this study was undertaken originally there was no time 
available to extend the computations into this range. 


EXPERIMENTAL WORK 


A number of strip gratings were measured in the strip 
transmission line, which is essentially a parallel plate 
region of finite width. The instrument and the method 
of measurement have been fully described elsewhere,!*”° 
so that for the present it is sufficient to say that it may 
be regarded as a large waveguide of about ten wave- 
lengths in width in which the field conditions approxi- 
mate to those in free space. 

The gratings were constructed by cutting out the 
apertures from a single sheet of 0.001-inch tin foil and 
supporting this between two blocks of polyfoam held 
together at their ends by adhesive tape (see Fig. 3). 


OO!" tin foil POLYFOAM 


TIN FOIL BENT AROUND 
POLYFOAM 


By wedging this structure tightly between the plates of 
the line, the contact with the line conductors was found 
to be very good, as judged by the consistency and repro- 
ducibility of the results. 

The measurements were made with a short-circuit 
termination to the line and then obtaining the usual 
Weisfloch curve” by moving the grating relative to the 
short circuit and noting the corresponding movement 
of the minimum in the interference pattern. The com- 
plex reflection coefficient follows at once.?° 

Figs. 4 and 5 show the magnitude and phase of the 
reflection coefficient of some typical gratings. The 
theoretical curves in these figures were obtained from 
the “Waveguide Handbook,” with allowances for the 
polyfoam support. The phase is referred to the grating 
plane. It is worth noting that the two curves which in- 
clude the points b=Xo reveal the inadequacy of the 
theory used. 


19M. M. Z. El-Kharadly, “Some experiments of artificial dielec- 
trics at centimetre wavelengths,” Proc. IEE, vol. 102, part B, pp. 
17-25; January, 1955. 

20 R. J. Primich, “A general experimental method to determine the 
properties of artificial media at centimetre wavelengths applied to 
an array of parallel metallic plates,” Proc. LEE, vol. 102, part B, pp. 
26-36; January, 1955. 

21 A. Weissfloch, “Ein transformationssatz tiber verlustlose vier- 
pole und seine anwendung auf der experimentelle untersuchung von 
dezimemter- und zentimeterwellenschaltungen,” Hochfrequenz. U. 
Elektroak., band 60, pp. 67-73; September, 1942. 
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Figs. 6 and 7 show the magnitude of the reflection co- 
efficient for two gratings, but the variational method 
has now been used. The “Waveguide Handbook” re- 
sults would give zero reflection coefficient at the reso- 
nance points b=Xo, and it is concluded that the varia- 
tional method is quite adequate for the region b=). No 
results are shown for 6> Xp as the strip transmission line 
is restricted to normal incidence only, due to the adverse. 
interference of higher-order propagating waves. Fig. 8 
shows typical Weissfloch curves*! in the presence and 
absence of these higher-order waves. It is evident that 
the second curve is of little use. 


Arbitrary Incidence 


At the author’s request, some preliminary measure- 
ments were made on several gratings by Sollom, who 
used a parallel plate spectrometer. This is essentially a 
parallel plate region in which provision is made to rotate 
the two horns about a fixed center.” Figs. 9 and 10 show 
two typical angular spectra. The theoretical curves were 


* P. H. Sollom, “An Accurate Parallel-Plate Spectrometer,” 
Ph.D. dissertation, University of London, England; 1954. 
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6;=angle of incidence 

6, =angle of emergence of higher-order wave. 

Fig. 11 shows a typical power spectrum for which the 
computation still has to be done. The peak in trans- 
mitted power at an angle of incidence off normal is of 
considerable interest as this shows the beginnings of a 
phenomena well known for optical gratings, 7.e., the 
series of transmission peaks at various angles of 
incidence.”3:4 


d= 88 cms 
diffracted waves present 


COMPARISON BETWEEN A STRIP GRATING AND A SEMI- 
INFINITE SYSTEM OF THICK PLATES 


It has been conjectured’ that the results for a strip 
grating could be used to compute the reflection proper- 


23 W. C. Meecham, “A Variational Method for the Calculation of 
the Distribution of Energy Reflected from a Periodic Surface,” Eng. 
Res. Inst., Univ. of Michigan, Ann Arbor, Mich.; November, 1955. 

% Since the writing of this paper, some work which has a direct 
bearing on this phenomena has come to the author’s attention. 
Meecham, ibid., has used a variational method to compute the re- 
flection from a blazed grating. A specific anomaly of Wood (loc. cit.) 
6 4 2 0 is treated and it is clearly evident that a peak in the main reflected 
l2 10 8 ) wave is associated with the first-order diffracted wave changing to 

tome an evanescent condition. Fig. 11 is an excellent experimental illus- 
Fig. 8. tration of this situation. 
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ties of a system of thick plates,® owing to the identical 
cross sections. This has been examined in a very general 
way® and no simple relation appears to exist. For in- 
stance for b<Xpo. the admittance of a thick plate system 
can be written 


Yin = [By +F 


where B, is the grating susceptance and F is a compli- 
cated function which bears no simple relation to the 
input admittance of a system of infinitely thin plates. 
Although the formal separation of the grating suscep- 
tance has been achieved, it is quite pointless as the func- 
tion F requires as much, if not more calculation than the 
original Vin. 

The general conclusion is that the grating results can- 
not be used to advantage in obtaining numerical results 
for a system of thick plates. However, owing to the simi- 
larity in the free space field structure of the grating and 
a system of thick plates, the definite possibility does 
emerge of using a grating to match a parallel plate lens, 
especially in the wavelength region where the first dif- 
fracted wave is present. This would permit such a Jens 
to be operated over a wider range of angle of incidence 
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with less interference from high-order waves, but this 
will not be pursued any further here. 


CoNCLUSION 


The variational method has been found to be very 
useful to obtain the properties of a strip grating. Its ac- 
curacy has been verified experimentally for all values of 
the parameter (b/Ao) <1. There is still a definite need 
for results for arbitrary incidence and in principle the 
variational method is still applicable. In fact, expres- 
sions have been obtained for the grating admittance for 
any angle incidence and it is hoped that these will be 
presented later. It would be instructive to compare these 
with Muller’s numerical computations, which, as far as 
the author is aware, are not yet available. 


ACKNOWLEDGMENT 


This work was carried out at Imperial College, in 
partial fulfillment of the requirements for the Ph. D. de- 
gree at the University of London, England. The author 
wishes to acknowledge the assistance of Prof. Willis 
Jackson and Dr. J. Brown, both formerly of the above 
college. Thanks are due to Dr. P. H. Sollom who did 
the measurements for arbitrary incidence. 


1957 


a 


_ IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


183 


The Dependence of Microwave Radio Signal Spectra 
on Ocean Roughness and Wave Spectra* 
C. I. BEARD} anp I. KATZt 


Summary—This paper is an extension of previous work on reflec- 
tion of microwaves from an ocean surface. The present analysis, 
dealing with spectra, is based on data obtained in a one-way X-band 
Propagation experiment performed across the Golden Gate, San 
Francisco. Two paths of 9000 and 15,000 feet were used. To describe 
the ocean surface, wave gages were mounted on a piling driven into 
the Golden Gate channel. 

Radio-signal spectra are found to be broader than the ocean-wave 
spectra and the spectral breadth a function of ocean roughness. The 
important result of this analysis is the establishment of a linear rela- 
tionship between ocean roughness and the spectral breadth of the 
radio signals. Ocean roughness is measured by the product of the 
standard deviation of the wave height and the grazing angle divided 
by the radio-wave length. Radio spectral breadths are determined 
by the frequencies at which each spectrum drops to the 0.9-, 0.8-, 
0.7-, 0.5-, 0.25-, and 0.1-power points. The breadths are then ex- 


_ pressed as ratios of these frequencies to the frequency of the peak 


in the simultaneous ocean-wave spectrum. 

The analysis now enables one to predict the approximate shape of 
the spectrum of the radio signal received in a one-way transmission 
path given only a knowledge of the geometry, radio-wave length, 
ocean-wave height, and the peak frequency in the ocean spectrum. 


INTRODUCTION 


HIS PAPER is concerned with the propagation 
le short-wave electromagnetic energy in over- 
water line-of-sight paths. In general, the aim is to 
obtain knowledge of the mean value of the signal 
strength, the amplitude of the fluctuations about the 
mean and spectra of the fluctuations. A certain measure 
of progress has been made in a previous paper’ to predict 
the mean signal level and the amplitude of its fluctua- 
tions from a knowledge of ocean-wave height and ge- 
ometry. The present paper reports on subsequent 
analysis along similar lines, but stressing primarily the 
relationship between ocean-wave spectra and spectra of 
received radio signals. It will be assumed in this paper 
that atmospheric inhomogeneities are of minor im- 
portance in the spectra since the measurements were 
'made over relatively short paths. 

Since the interpretation of the spectra is based on the 
vector model representation! of the process of micro- 
wave reflection from the ocean surface, a brief summary 
of the model is given first. The subsequent section then 
describes the experimental conditions. A section on the 
general features of the spectra serves as an introduction 


* Manuscript received by the PGAP, July 9, 1956; revised manu- 
script received, October 26, 1956. ‘ F ; 

+ Sylvania Electronic Defense Lab., Mountain View, Calif. 
Formerly at Appl. Phys. Lab., Johns Hopkins Univ., Silver Spring, 
Md 


t Appl. Phys. Lab., Johns Hopkins Univ., Silver Spring, Md. 

1C, J. Beard, I. Katz, and L. M. Spetner, “Phenomenological 
vector model of microwave reflection from the ocean,” IRE TRANS. 
vol. AP-4, pp. 162-167; April, 1956. 


for the section on the details of the quantitative relations 
which were found experimentally. A qualitative vector 
model interpretation is followed by comparison of re- 
sults of an approximate theoretical calculation. 


THE VECTOR MopDEL 


In the formulation of the vector model it is assumed 
that the radio signal may be represented as the vector 
sum of a direct signal and a water-reflected signal, each 
with certain properties. The direct ray, D, is assumed 
constant in magnitude and phase, whereas the reflected 
signal is composed of a coherent part, C, whose ampli- 
tude and phase are fixed by geometry and surface 
roughness, and an incoherent part, J, of random ampli- 
tude and phase. Except under special conditions, which 
will be mentioned later, J may be assumed to be the 
resultant of vectors from a large number of independent 
random scatterers. The incoherent part, J, may be de- 
composed into two orthogonal independent vectors each 
representing a Gaussian process with zero mean. It is 
convenient to orient these two orthogonal vectors along 
and normal to C; in addition, their standard deviations 
and spectra are assumed to be equal. Both C and J are 
related to geometry and sea conditions by use of the 
roughness parameter hy/\, where h is the standard de- 
viation of the water surface, W is the grazing angle, and 
d\ the wavelength of the electromagnetic radiation. The 
quantitative dependence of both C and J on hy/d has 
been presented.! 


THE GOLDEN GATE EXPERIMENT 


During April, 1953, a propagation experiment was 
performed by the Applied Physics Laboratory, in con- 
junction with the Electrical Engineering Research 
Laboratory of The University of Texas, across the 
Golden Gate at San Francisco. In the experiment two 
paths were used as shown in Fig. 1. At one end of the 
path marked “slope,” a transmitter was placed at vari- 
ous heights from 92 feet to 900 feet on the slope of a 900- 
foot mountain adjacent to the Golden Gate (the twelve 
sites lie along the line denoted by the numbers 1, 11, and 
12 on the map). The signal was received at a point about 
85 feet above water on the other side, a transmission 
distance averaging about 9000 feet. On the other path 
to Bonita Cove, 15,000 feet long, transmitter heights of 
75 feet to 360 feet (denoted by A, B, C, and D on the 
map) were available. This second path was chosen to 
obtain a greater variety of conditions, 7.e., a different 
angle to the ocean swells and a different path length. 
The Golden Gate has an unblocked westerly exposure 
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to the open ocean and the prevailing swell and wind 
directions are northwesterly. For the slope path the 
predominant wave directions varied from southwesterly 
to westerly. 

The transmitting antenna pattern was chosen to 
minimize illumination of the mountainside and yet to 
illuminate the water adequately. It was a 12-inch diame- 
ter paraboloid, at X band, with an 8° beamwidth at 
half-power points and side lobes down 24 db. The trans- 
mitting antenna was set on tripods at each site and was 
pointed optically at the receiver. Under these conditions 
more than 20 Fresnel zones were illuminated within the 
half-power points. 
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Fig. 1—Map of Golden Gate experiment area. 


The receiving antenna was an 18-inch diameter 
paraboloid with a 5.5° beamwidth between half-power 
points and side lobes down 30 db; it was placed on an 
85-foot high ledge where the ground dropped sharply 
away in front and to the sides of the antenna. To vary 
the illumination, this antenna was pointed at angles 
from zero to 5° above the horizontal. Both vertical and 
horizontal polarization were used in the experiment but 
only the vertical polarization spectra are analyzed in 
this paper. 

Ocean-wave heights were measured by three types of 
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gages mounted on a piling driven into the channel in 
about 35 feet of water off the northern end of Baker 
Beach.? Water-surface fluctuations between 0.3- and 
1.1-feet rms were experienced during the course of the 
experiment. 

Of interest in this paper is the method of data re- 
duction and analysis. Ocean-wave heights were recorded 
on magnetic tape. The radio signals were recorded on 
Brush strip-chart recorder, then curve followed later in 
the laboratory to put the information on magnetic 
tape. Because the fluctuations were of such long periods 
it was necessary to rerecord twice to obtain a speed-up 
of about 500. Such speed-up permitted utilization of a 
harmonic analyzer for obtaining spectra.* Tapes were 
spliced into continuous loops and played into the 
analyzer as the analyzer frequency was varied. The 
spectra had a minimum of 38° of freedom corresponding 
to a 4-minute length of record and a 0.04-cps wide band- 
pass filter. Many of the wave records were 15 minutes 
in length with a proportional increase in the degrees of 
freedom. 


GENERAL FEATURES OF OCEAN-WAVE AND 
RADIO-SIGNAL SPECTRA 


Ocean-Wave Spectra 


The observed ocean-wave spectra may be placed in 
two general classes: single-peaked and double-peaked; 
spectra are considered to be double-peaked when a 
second peak is within 3 db of the first maximum. Some 
spectra show incipient second peaks on the higher fre- 
quency side which undoubtedly indicate the presence of 
wind waves of insufficient power relative to the swell to 
show up as definite peaks in the spectrum. Figs. 2 and 
3 are examples of the single- and double-peaked spectra. 
Fig. 4 is illustrative of a spectrum in which the presence 
of wind waves is indicated at about 0.13 cps. 

Of the data runs resolved into spectra, there are cases 
in which significant peaks occur at frequencies as low as 
0.065 cps and as high as 0.12 cps. Half-power “band- 
widths” of the wave spectra range between 0.03 cps and 
0.17 cps. 

A useful analytic form for wave spectra developed by 
Neumann is? 


A*(f) = f° exp —2(g/2nfuv)? 


where A*(f) is the power spectrum, f the frequency, g 
acceleration of gravity, and v the wind speed. This ex- 
pression may be rewritten 


A*(f) = f-* exp (—K/f?) 


* K. H. Jehn, J. R. Gerhardt, D. F. Metcalf, and S. J. Prosser, 
“Some Meteorological and Oceanographic Characteristics of the 
Golden Gate, California, Area,” Elec. Eng. Res. Lab., Univ. of Texas, 
Rep. CM- 760; February 26, 1954. 

s wee E Staake, oN \ Low. Feeauency pecs and Amplitude- 
istribution Analyzer, s. Lab., Johns 

CM-833; March, 1088. me 4 : Heenan ears 
4W. i Pierson, Jr., “Advances in Geophysics,” 


Academic Press, 
Inc., New York, N. Y., vol. 2, pp. 140-146; 1955. 
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Fig. 2—Example of single-peaked ocean-wave spectrum. 
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Fig. 3—Example of double-peaked ocean-wave spectrum. 
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Fig. 4—Example of an ocean-wave spectrum indicating presence 
of wind waves at about 0.13 cps. 
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in which only the frequency dependence is explicitly 
utilized. The maximum of A®(f) occurs at a frequency 
fue = [K/3} 

The sample single-peaked ocean-wave spectrum in 
Fig. 2, normalized for unity at the peak, is shown again 
in Fig. 5 for comparison with its corresponding Neu- 
mann spectrum. The Neumann spectrum was “broad- 
ened” by computation to simulate its passage through 
a 0.04-cps bandwidth filter in order to compare it ade- 
quately with the observed wave spectrum. It may be 
seen that the measured spectrum fits the Neumann 
spectrum except for the somewhat higher energy content 
at the low frequencies. 
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Fig. 5—Comparison of single-peaked ocean-wave spectrum with 
corresponding Neumann spectrum. 


Of the many ocean-wave recordings obtained at the 
Golden Gate, 24 are used in this analysis. They are di- 
vided as indicated below. 


1) Single-Peaked Spectra: The observed ocean-wave 
spectra were compared at the half-power levels 
with broadened Neumann spectra with the follow- 
ing results: 5 of 12 have half-power widths within 
+10 per cent of the corresponding Neumann spec- 
tra, 7 of 12 within +20 per cent, and 10 of 12 with- 
in +30 per cent. For some analytical applications, 
then, it may be adequate to use the Neumann 
spectrum rather than a measured spectrum when 
single-peaked spectra are needed. 

2) Double-Peaked Spectra: In these cases a reexami- 
nation of the original strip-chart record discloses 
two separate easily identifiable “trains” of waves 
at the frequencies corresponding to the two major 
spectral components. 

A possible mechanism for this type of spectrum 
is that the swell corresponding to the lower fre- 
quency peak comes a comparatively long distance 
from a generating area where the effective wind 
speed is Vj. In traversing this distance, the swell 
may pass through a nearby generating area where 
local winds of speeds V2 superimpose a higher fre- 
quency set of wind waves, resulting in the higher 
frequency peak of the spectrum. 
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Radio-Signal Spectra 


Radio-signal spectra are found to be characteristically 
broader than the wave spectra, half-power points in 
these spectra extending out to 1 cps or more. In these 
broader spectra there appear many minor spectral 
peaks, apparently harmonics of fundamentals in the 
ocean-wave spectrum. For purposes of analysis, to ob- 
tain the equivalent of more degrees of freedom, smooth- 
ing by computation has been performed on the radio 
spectra to simulate passing the signal through a 0.2-cps- 
wide band-pass filter. A sample radio spectrum is shown 
in Fig. 6. The original spectrum replotted from the 
spectrum-analyzer record is shown as a solid line and the 
smoothed spectrum is dashed. Additional features of 
radio spectra will be discussed in the next section. 
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Fig. 6—Sample radio spectrum as obtained from spectrum analyzer 
tes 0.04-cps filter) and smoothed spectrum (with 0.2-cps 
ter). 


QUANTITATIVE RELATIONSHIPS BETWEEN OCEAN-WAVE 
AND RADIO-SIGNAL SPECTRA 


It has been demonstrated by the Electrical Engineer- 
ing Research Laboratory’ that the peaks in radio spectra 
generally occur at the same frequency as the peaks in 
the simultaneous ocean-wave spectra. In the present 
analysis, it is found that maxima in the radio spectra 
appear on some occasions at frequencies which are 
second and possibly third harmonics of the ocean-wave 
spectra maxima in addition to occurring at the ocean- 
wave fundamental. The smoothing process, however, 
blurs the fine structure around the maxima making their 
positions rather indefinite. 

Broadening of the radio-signal spectrum with in- 
creasing roughness was quite evident. When runs were 
made during one day in which the grazing angle, y, was 
progressively increased but wave height, 4, was almost 
constant (yielding increasing MW/d), the spectra broad- 
ened successively. On different days of the experiment 
runs were made at the same grazing angles but under 
different wave-height conditions; again broader radio 


5 W. J. McKune and H. W. Smith, “Comparison of Power Spec- 
trum Estimates of Overwater Microwave Signals and Associated 
Water Waves,” Elec. Eng. Res. Lab., Univ. of Texas, Rep. No. 68; 
May 6, 1953. 
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Fig. 7—Radio spectra showing broadening with increasing 
roughness. 


spectra resulted when hy/) was larger. Fig. 7 shows four 
smoothed radio spectra, under conditions of M/d from 
25 to 590. The equivalent smoothing filter used to obtain 
these curves was 0.2 cps wide except for the hy/A=25 
curve for which the bandwidth was cut in half. For this 
latter spectrum, where the spectral width is narrower 
than the filter bandwidth, using an 0.2 cps filter width 
would have masked the true spectrum; hence a narrower 
filter was used. 

A set of empirical curves has been obtained which in- 
dicates the relationship of signal-strength spectra and 
peak frequencies of ocean-wave spectra to apparent 
ocean roughness hW/\. These are shown in Figs. 8, 9, 
and 10. The ordinate in Fig. 8 is the ratio, f,/fw, of the 


400 
ROUGHNESS, hw/) (mils) 


Fig. 8—f-(0.5)/fwaves vs hy/d. f,(0.5) =frequency of 4-power point 
of radio spectrum. fwaves=frequency of maximum of wave spec- 
trum. ©: Represent ratios when ocean-wave spectra are single 
ac et A: Represent ratios when ocean-wave spectra are double 
peaked. 


frequency at which the radio-signal spectrum has de- 
cayed to 3 power to the frequency of the lower frequency 
peak in the ocean-wave spectrum. The frequency at the 
3 power pojnt on the radio spectrum is here denoted as 
f-(0.5). The circles are values of the ratio f,/f, when the 
ocean-wave spectra are single peaked while the triangles 
are for double-peaked ocean spectra. The lower fre- 
quency peak of the double-peaked spectrum is used to 
compute the points. In Fig. 9 the 1/10-power frequency 
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400 
ROUGHNESS, hw/A (mils) 


Fig. 9—f,(0.1) /fwaves v8 hw/d. f-(0.1) = frequency of 1/10-power point 
of radio spectrum. Swaves=frequency of maximum of wave spec- 
trum. ©: Represent ratios when ocean-wave spectra are single 
ems A: Represent ratios when ocean-wave spectra are double 
peake 


of the radio spectrum is used. In these figures the }- and 
1/10-power frequencies are on the high-frequency side 
of the peak in the spectrum. The plots in Figs. 8 and 9 
giving all the individual data points illustrate the be- 
havior shown by the data for other power levels. The 
results of such plots for power levels from 0.9 to 0.1 are 
summarized in Fig. 10, which shows the straight lines 
fit to the data points for each power level. Data on the 
values of f, and f, for all the runs used in this analysis 
are given in TableI,on page 188. A discussion of the possi- 
ble sources of error in the data presented here, including 
the effects of sample length, representativeness of par- 
ticular wave spectra, and reliability of the spectra at the 
low-frequency end, has been reported.® 

The effect of the parameter Ay/d overrides the in- 

fluence of other possible factors: 

1) The points for double-peaked wave spectra fall 
close to those for single peaks when the lower fre- 
quency peak, or swell frequency, is used to, com- 
pute the points. This behavior indicates that the 
width of the ocean-wave spectrum is of less im- 
portance in determining the width of the radio 
spectrum than the swell frequency. 


6 C, I. Beard and I. Katz, “The Dependence of Microwave Radio 
Signal Spectra on Ocean Roughness and Wave Spectra,” Appl. 
Phys. Lab., Johns Hopkins Univ. Rep. CF-2456; January 17, 1956. 
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Fig 10—Composite curves of fradio/fwaves VS hy/d for various 
power points of the radio spectra from experimental data. 


2) Illumination of the water was varied radically by 
changing the pointing of the receiving 18-inch- 
diameter antenna from zero to 5° above the hori- 
zontal without changing the spectrum width. 
Under-illumination of scatterers far from the spec- 
ular point (by pointing the antenna upwards) 
should have little effect on the spectra according 
to the next section on theory which considers only 
the area near the specular region. 

Although the low hW/d ends of the lines in Figs. 8, 9, 
and 10 have been left open, lower limits of the frequency 
ratios as h/\—0 can be specified in the following man- 
ner. One assumes that the radio spectrum does not be- 
come narrower than the wave spectrum as hy/\-0. 
Then for the Neumann wave spectrum the ratio of the 
upper half-power frequency to the peak frequency is 
1.55; the lowest point in Fig. 8 is 1.8 which does not 
violate the lower limit. For Fig. 9 the Neumann ratio is 
2.3 and the lowest point is 4.4. Hence there is no dis- 
agreement with this criterion in the curves. 

An analytic expression for the shape of the smoothed 
radio spectra was sought. The portions of the curves 
above the maxima of the radio spectra, fn, follow the 
form exp [—R(f,—fm)"]. Seventy-five per cent of the 
values of » fall between 1 and 1.5 with an average value 
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TABLE I 


i = 
ee SSS Sa 


Wave Spectra 


Radio Spectra 


yk eee mils) | Frpauency of | (0.9) | (0-8) | f(0.2) | (0.5) | (0.25) | f,(0.1) 
tes te ae Moe (cps) (cps) (cps) (cps) (cps) 
Single Peak 
0.63 0.84 
11 1328 175 0.090* 0.19 0.21 0.28 0.42 
11 1339 200 0.090 0.24 0.27 0.30 0.36 ee Pe 
11 1410 225 0.090* 0.19 0.21 0.38 0.53 oe pee 
18 1134 130 0.120 0.18 0.20 0.21 0.25 ae ai 
21 1102 210 0.115** 0.38 0.46 0.54 0.76 a oe 
21 1001 235 0.115 0.16 0.20 0.30 0.48 ae a 
21 1444 290 0.100 0.45 0.51 0.56 0.80 ai Aes, 
22 1008 390 0.100 0.52 0.59 0.65 0.86 Be ret 
22 1157 305 0.075 0.35 0.42 0.47 0.66 oe ‘io 
oD 1218 340 0.080 0.42 0.47 0.52 0.76 i He: 
22 1235 360 0.075 0.36 0.44 0.58 0.83 ae a 
22 1301 320 0.090 0.39 _ 0.94 1.06 Bs ee: 
22 1315 410 0.085 0.38 0.47 0.80 0.88 f ne 
23 1458 90 0.105 0.16 0.19 0.22 0.28 0.38 = 
24 1016 130 0.100 0.19 0.25 0.28 0.34 0.42 0. 
Double Peak 
034 235 0.085 0.32 0.36 0.39 0.44 0.56 0.92 
if 1055 115 0.085*** 0.21 0.23 0.25 0.28 0.38 0.47 
18 1107 140 0.105 0.18 0.21 0.24 0.32 0.52 0.72 
18 1229 170 0.095 0.19 0.26 0.32 0.41 0.59 0.73 
22 1100 205 0.075 0.25 0.31 0.38 0.56 0.84 1.10 
ig) 1144 270 0.080 0.39 0.43 0.46 0.61 0.94 1.29 
22 1412 410 0.080 0.68 0.78 0.89 tei 1.86 2.46 
DD 1448 590 0.090 0.74 0.88 27 1.58 2.20 2.95 
D7, 1525 435 0.100 0.54 0.65 0.77 0.94 1.51 2.38 
24 1009 125 0.070 0.18 0.20 0.23 0.31 0.50 0.65 
24 1041 155 0.080 0.20 0.23 0.27 0.36 0.59 0.80 
24 1106 175 0.070 O17 0.23 0.30 0.47 0.66 0.83 
25 1502 25 0.065 0.09 0.10 0.10 0.12 0.20 0.29 


* The closest wave spectrum was at 1340-1356", 
4k The closest wave spectrum was at 1111-1115. 
*4# The closest wave spectrum was at 1020-1036". 


of about 1.25. There appears to be no dependence of n 
upon hy/d. The parameter k, however, varies inversely 
with hy/d as given by the approximate relation: 


» 90 < hb/r < 590. 


Hence, an average expression for 75 per cent of the 
radio spectra on the high-frequency side of the peak is 
given by: 


W(f) = exp — |(-05 + aan) = ihe 


As h/d increases outside of the limits given above 
towards infinity, it appears reasonable to assume that k 
tapers off to zero (rather than to the straight line inter- 
cept value of —0.5). Thus as hy/A—> ~, 


WY elL— kG. — fa) 3A, 


or the spectrum approaches a white-noise spectrum. 
The explanation for the spectral broadening with in- 

creasing roughness may be found in concepts inherent in 

the vector model. Earlier work by the Electrical Engi- 


neering Research Laboratory’ has shown that harmonics 
of the ocean-wave fluctuations are found in the radio 
spectrum, having been generated by the vertical motion 
of the ocean surface. Various ratios of fundamental to 
harmonics were possible. In their treatment the entire 
surface was assumed to move up and down as a unit, 
moving the interference pattern past the receiving 
point. The writers of the present paper conclude that 
the entire surface can be considered to move as a unit 
only under the following conditions: The transmission 
path is parallel to the swell crests with short ranges and 
beams narrow enough so that the entire illuminated 
area and important Fresnel zones are contained within 
one swell wavelength. This case may be achieved across 
narrow inlets from the sea. 

In the more general open-ocean situation, however, 
waves arrive at a point from various directions, up to 
+ 30° is common; this diversity of direction of arrival 
causes the waves to be short crested. As the Fresnel 
zones become longer than a few thousand feet, the short 
crestedness results in the existence of many scatterers 


7 A. H. LaGrone, A. W. Straiton, and H. W. Smith, “Synthesis of 


radio signals on overwater paths,” IRE Trans., vol. AP-3, pp. 
48-52; April, 1955. 
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within the Fresnel zones. Irregularity of the wave struc- 
ture as indicated by the width of the frequency spec- 
trum also aids in breaking up the Fresnel zones. Using 
the vector model as a guide one may interpret the 
phenomena as follows: Fluctuations in the radio sig- 
nal are caused by the vector addition to the direct and 
coherent signal of signals reflected from several or many 
reflecting areas. The model in its simplest form postu- 
lates that there are many independent, random scatter- 
ing areas. [In the special case of the narrow inlet cited 
above where the reflecting areas are greater in size than 
the Fresnel zones, the “incoherent term” loses its ran- 
-domness and the postulate of the model of many (>5) 
independent random scatterers is no longer satisfied.] 
As each small scatterer moves vertically through any 
height it changes its contribution to the phase of the re- 
sultant signal [for example, the phase of the return from 
a patch near the specular point will vary as 4r(hy/X) as 
h changes]; if it moves through twice that height in the 
same time its contribution to the frequency of fluctua- 
tion is doubled. Thus as wave height increases one 
would expect an increased spectrum bandwidth gener- 
ated by harmonics of all the frequencies in the wave 
spectrum. Actually minor peaks in the radio spectra 
have been found which are in harmonic ratios—up to 
the sixth have been resolved—of the frequency of the 
peak in the wave spectrum. Since the integrated effect 
of signals from all the scatterers, especially after being 
passed through the band-pass filter, might be expected 
to smear out the effects of individual wave frequencies, 
this effect was looked into. In many cases it was found 
- that clear-cut trains of waves of 1 to 14 minute duration 
of practically constant frequency existed near the time 
of the data run. The persistence of one-wave frequency 
would account for the predominance of its harmonics in 
the radio spectrum. In this picture all the fluctuations 
are considered to stem from the incoherent portion of 
reflected energy; the coherent part, and consequently 
the interference pattern, are fixed by geometry, as 
stated previously. 


APPROXIMATE THEORETICAL TREATMENT 


The qualitative explanation of spectral widening 
given in the preceding section can be made more quanti- 
tative, although still not exact, as follows. 

The reflected signal (coherent-plus incoherent) from 
each individual scatterer has been shown by Spetner® to 
be proportional to 


etku(yy + Yo) 


where k=27/\, y=height of scatterer above mean 
water level, ¥,=grazing angle of incidence from source, 
and 2.=grazing angle of reflection to receiver. 


8, M. Spetner, “An Independent Scatterer Model for F orward 
Scattering Rice acces Appl. Phys. Lab., Johns Hopkins Univ. 
Rep. CF-2533; July, 1956. 
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The assumptions involved in this simplified form are: 

1) Uniform illumination of the water surface by both 

transmitting and receiving antennas. 

2) Each scatterer has a constant scattering cross sec- 

tion and radiates isotropically. 

3) Each scatterer has the same reflectivity. 

The total value of reflected signal is obtained by inte- 
grating over the area. In order to obtain a solution, it is 
desirable to assume that the principal contribution is 
from scatterers near the specular region. One justifica- 
tion for this assumption is that it leads to reasonable 
results. 

The spectrum is obtained by first calculating the cor- 
relation function R(r) =(J,(¢)I,(¢+7)), where J, is the 
component of J along c. In calculating the averages, use 
is made of the known Gaussian probability distribution 
of the water surface. The result Spetner obtains is:8 


R(t) = exp [(4rhp/d)*ou(r) — yz] 


where pw(7) = (y(t)y(E+7)) is the correlation function of 
the water surface and y¥ is the probability per unit time 
that a given scatterer disappears. 

The power spectrum W(f) of J,(¢) is then computed 
by the cosine transform relation 


W(f) = 4 f Re cos wrdr. 


The integrations were performed numerically taking 
y=0.1 cps, the predominant wave frequency, and an 
assumed expression for py(r) of a damped cosine form 
which decayed to e! in two water wavelengths.® The 
I, spectra resulting from these calculations are given in 
Fig 11, 

Before discussing these J, spectra it is necessary to 
consider which spectra should be compared to the ob- 
served total-signal spectra. This can be seen by the 
series expansion of the signal strength, 7, derived by 
Beard, et al.,! within +90° of maxima of the interference 
pattern: 


Ia pis Wage 

Rape (a +) 
ive 
~ 8D 


where D,=|D+C|. 

To the first order of approximation, the spectrum of 
the fluctuating part of 7 is that of Jy. 

The frequencies of the 0.5-, 0.25-, and 0.1-power 


" 2 

9 The expression used in the computations was pw(r) =e Ow") 

cos wr. The wave frequency fw was selected as 0.1 cps and 6 was 

chosen as 0.5 to cause a decay to e! in two periods. For a theoretical 

derivation of this form of pw(r), see Carl Eckert, The generation of 

wind waves on a water surface,” J. Appl. Phys., vol. 24, pp. 1485- 
1494; [eq. (49)]; December, 1953. 
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Fig. 11—Power spectra of I, for hy/A=200, 300, and 600 mils 
calculated from theory. 


ROUGHNESS, hw/d (mils) 


Fig. 12—Broadening of power spectrum of J, with increasing 
roughness, as computed from theory. 


points of the J, spectra (each divided by f,,) obtained 
from Fig 11 are plotted in Fig. 12 vs hW/. These plots 
show the same linear behavior as found experimentally. 
However, the theory does not predict a continuation of 
this linear behavior at power levels closer to the spectral 
peaks. This can be seen by plotting these points from 
Fig. 11. But since the behavior for these levels is not 
linear, reliable curves could not be drawn with points at 
only three values of M/) available. 

In addition to this linear broadening with hy/d, the 
relative shapes of the calculated J, spectra are approxi- 
mately the same (between the 0.1- and 0.5-power points) 
as the experimental spectra, as shown in Table II. 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


TABLE II 
Ip /\=600 ny /r=300 
fai Calc. Exp. Calc. Exp. 
ae 1.35 1.39 143% 1.38 
f-(0.25) ; 
7(0.25 
fee eed 15538 1.45 1.49 
f.(0.5) : 
The ratio 
fr(0.1) 
rm ay ea, 
f,0.25) 
for example, is actually the ratio 
f-(0.1)/fo 
wy erp G sie erst 
JA O22) ile 


but fy cancels. 

Although the trends agree, the widths’ magnitudes of 
the J, spectra are only about one-half the experimental 
widths. This is mainly a result of choice of the wave cor- 
relation function.’ Inspection of the wave correlations 
determined by the Electrical Engineering Research Lab- 
oratory!” shows that a simple exponential damping 
e— fut) fits better although it is still inadequate." Values 
of b range from 0.7 to 1.4. With 6=1, the correlation has 
decayed to e—! in one water wavelength instead of two 
as originally assumed. This would widen the spectra by 
about a factor of two, placing the theoretical spectra in 
agreement with the experiment. 

As will be noted from Fig. 11, the theoretical spectra 
peak at zero frequency rather than drop towards zero 
at zero frequency as evidenced by experimental data. 

The theoretical calculations were made only for cor- 
relation functions corresponding with single-peaked 
wave spectra. When double-peaked wave spectra oc- 
curred, the data indicated that swell frequencies were 
more important in determining the width of the radio 
spectra than the higher frequency wind waves. In order 
to establish the quantitative relative effectiveness of 
swell vs wind frequencies the next step is to apply Spet- 
ner’s expression for R(r) using correlation functions 
corresponding to double-peaked wave spectra. 

The present simplified theory agrees with the experi- 
mental trends of the broadening of the radio spectra 
with hy/d and also of the relative shapes of the radio 
spectra between the 0.5- and the 0.1-power levels. 


10 F, E. Brooks, Jr., G. P. Dubose, Jr., and C. W. Tolbert, “Pre- 
liminary Summary Report of Radio, Meteorological and Oceano- 
graphic Data Obtained During 1955 Gulf of Mexico Propagation 
an HS Eng. Res. Lab., Univ. of Texas, Rep. TG-260; January 

1. There is a wide variation among correlations. The exponential- 
ly-damped cosine form above would reasonably describe two, whereas 
the sin «/x form would describe another correlation. Since the Neu- 
mann spectrum fits many observed wave spectra, the cosine trans- 
form of this would be a reasonable correlation function to use. 
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CONCLUSION 


If one assumes a similarity of wave and radio spectral 
shape on the low-frequency side and that the spectra 
have primary peaks at the same frequency, there re- 
mains only the region on the high-frequency side to be 
predicted. The graphs of f,/f,, in Fig. 10 may be used for 
this purpose. Within limits of accuracy as indicated by 
the spread of the points in Figs. 8 and 9, one is now ina 
position to predict the approximate shape of the radio 
spectrum from a knowledge only of geometry, electro- 
magnetic wavelength, ocean-wave rms height, and the 
(swell) peak frequency in the ocean-wave spectrum. 
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Step Discontinuities in Waveguides* 
W. ELWYN WILLIAMS} 


Summary—The method developed by Wiener and Hopf for solv- 
ing a certain class of integral equations is applied to the problem of 
a discontinuity in cross section of a rectangular waveguide. We con- 
sider the problem of two waveguides of infinite width and different 
heights joined together with a short step when the only nonvanishing 
component of the magnetic field is that parallel to the edges of the 
step. It is assumed that the fundamental mode is incident in the 
larger channel traveling towards the step. For values of the channel 
height such that only the fundamental mode is propagated without 
attenuation, this problem has been treated in the ‘‘Waveguide 
Handbook”!; our results are in very good agreement with the results 
given there. The present analysis is extended to cover the case where 
the larger channel is capable of supporting the mode above the funda- 
mental without attenuation. Some numerical values are given for 
the ratio between the reflected and incident energies for various val- 
ues of the channel and step heights. 


INTRODUCTION 
Ape filters 0 involving discontinuities in wave- 


guides are of great practical importance and a con- 

siderable amount of attention has been devoted to 
them in recent years. The two types of discontinuities 
that have been most extensively treated for a rectangu- 
lar guide are the bifurcation of the guide by a semi-infinite 
plane and a step discontinuity in the guide cross section. 
The problem of the bifurcation may be solved in terms 
of known functions by applying the well-known Wiener- 
Hopf technique for solving integral equations. The 
methods for solving the problem of a step discontinuity 


* Manuscript received by the PGAP, May 17, 1956. The research 
reported in this article was done at the Inst. of Mathematical Sciences, 
New York Univ., and has been made possible through support and 
sponsorship extended by the U. S. Air Force, Air Force Cambridge 
Res. Center, under Contract No. AF 19(122)-42. 

+ Inst. of Mathematical Sciences, New York Univ:, New York, 

Y. 


NAY 
1N. Marcuvitz, “Waveguide Handbook,” McGraw-Hill Book 
Co., Inc., New York, N. Y.; 1951, 


which have appeared in the literature are, however, in 
no way based on the Wiener-Hopf approach. Two dis- 
tinct methods have been employed; the first consists of 
expanding the field components on either side of the 
discontinuity in terms of the appropriate eigenfunctions 
for the region considered, and then matching the com- 
ponents across the plane of discontinuity. The second 
method used is the equivalent static approach as de- 
scribed in the “Waveguide Handbook.”! Since there 
is obviously a physical relationship between these two 
problems it is desirable to solve them both by the same 
method in order that the solution should display clearly 
the relationship between them, and this is the object of 
the present work. 

We therefore consider the application of the Wiener- 
Hopf technique to solve the problem of a step discon- 
tinuity in a parallel-plate region (or, equivalently, the 
coupling of two parallel-plate waveguides of different 
height). It is assumed that the only nonvanishing com- 
ponent of the magnetic field is that parallel to the edges 
of the step. It is assumed that the fundamental mode is 
incident from the larger region on the step, and, initially, 
that the dimensions of both regions concerned are taken 
to be such that only the fundamental mode may propa- 
gate without attenuation. This restriction is later re- 
moved and the case of the larger region being capable of 
sustaining two nonattenuated modes is treated. The 
advantage of this method of attack is that it affords a 
clearer physical picture of the structure of the solution 
than that obtained by the other methods. In particular, 
it is possible to separate the reflection coefficient for the 
fundamental mode into two parts, one of which is ex- 
pressible in terms of known functions and another which 
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may be obtained from the solution of an infinite set of 
equations. That part of the reflection coefficient which 
is known exactly is by far the dominant part of the com- 
plete expression, and thus the percentage error in deter- 
mining an approximate value of the reflection coefficient 
is much less than the corresponding error made in the 
solution of the equations. The analytic part of the solu- 
tion may be obtained by simple physical arguments 
from the known solution of the bifurcation problem and 
thus the present form of the solution gives a distinct con- 
nection between the two types of problems. 

For the case in which only the fundamental mode is 
propagated without attenuation in any of the regions 
considered, it is possible to compare the present results 
with those given in the “Waveguide Handbook,” which 
were obtained by the equivalent static method, and the 
results of both methods show remarkably good agree- 
ment. The method used here may be easily extended to 
solve the problem of a bifurcated waveguide when the 
material composing one of the bifurcated regions differs 
in dielectric properties from that of the remainder of 
the waveguide. The approach may also be used to con- 
sider the problem of radiation from a dielectric wave- 
guide, and it is hoped that both of these problems will 
be treated in future papers. The method can, of course, 
be applied to the general electromagnetic case since it 
is known that this problem can be reduced to the solu- 
tion of two scalar problems,? but no attempt is made to 
consider this more general case. The author has been 
informed by Prof. N. Marcuvitz that Prof. J. Schwinger 
has considered the application of the Wiener-Hopf tech- 
nique to problems similar to the one treated here. An 
account of his work may appear in a forthcoming book 
by Marcuvitz and Schwinger. 


STATEMENT OF THE PROBLEM 


The waveguide is assumed to occupy the region 
—«o= X<0, d<y<a; and 0<x<~, 0<y<a. (see 
Fig. 1), where x, y are Cartesian coordinates with their 


Fig. 1. 


origin at O. Oz, the third axis of the system, is not shown, 
and the solution is assumed to be independent of z. A 
time variation e*”¢ will be taken as understood through- 


2 L. Lewin, “Advanced Theory of Waveguides,” Iliffe and Sons, 
London, England; 1951. 
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out the work, and it is assumed that an incident wave 
e** is approaching the step from the right, k being 
27 / (wavelength). 

The total magnetic field will be written as e#*+YV, 
and hence a function V is sought which satisfies the 
following conditions: 


pt (V2 + k)Y = 0, (1) 
—= — ik on + = 0, O0<y<d (2) 
Ox é 
ow 
——$ = on y= 4, —-n<4< 0, 
dy 
on y = d, —«o<x<0,and (3) 
on y = 0, OS Ee, 


It is convenient to assume that k=k,—7k;, where k;, 
k,>0; then the form of the radiation condition for V is 
v—Ofe-*l) as | «| —o. In order to specify the prob- 
lem uniquely, a condition has to be imposed on V at the 
edges of the step.* This condition is that VW is O(1) as 
the point of observation approaches (0, d) or (0, 0), 
while |VW| is O(r—/%) as the point of observation ap- 
proaches (0, d), and is O(1) in the neighborhood of the 
origin, here —r is being the distance from the point ap- 
proached. The above conditions now specify the prob- 
lem completely and uniquely. 


SOLUTION FOR V 
We shall now define a function ¢(s, y) by 


6,9) = fo we, de 
where s is the complex variable o+77. From the as- 
sumption on V as |x|, it is clear that the above in- 
tegral is convergent for —k;<o0 <k;. Hence @ is a regular 
function of s in the region |o| <k;. We also define ¢_(s, 
¥), b+(s, 9) as 


¢—(s, y) = i Ua, y)dx 


—eo 


eC nAS ff e0 (x, y)da. 
0 


Clearly, from the behavior of V as | «| —o,¢_(s,y)isa 
regular function of s for o<k;, and $4(s, y) is a regular 
function of s for o> —&;. In the region.0<y<d, $_(s, y) 
is identically zero and hence we have the following 
equations. For y>d 


ao 
dy? + «Ko = 0 (4) 
For 0<y<d 
dh 
eh ete Sate ay) (5) 


5D. S. Jones, “A note on diffraction by edges,” 
and Mech., vol. 3, pp. 189, 1950. bier rae Sinead aie 
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where f(y) =Y(0, y) and «= (s?+2)!/2, and where that 
branch of the root is taken which reduces to & when 
s=0, 

The solution of (5) which satisfies the boundary con- 
dition on y=0 is 


ik 4 
$4(s, y) = Acosxy — =| sin x(y — #)dt 
K Jo 
s son 
+ =f f(d sin K(y — Adt 
Kk Jo 


tk 
mA BOS gy one (t — cos Ky) 
K 


+ - f peAy in Reh = De (6) 


where A isa function of s but not of y. 
If primes are taken to mean derivatives with respect 
to y, then differentiation of (6) and then substitution of 
-y=d in the resulting equation gives 


tk 
$,'(s, d) = — xA sin xd — — sin xd 
K 


d 
+ sf f( cos «(d — #)dt. (7) 
0 
If the value of A from (7) is substituted in (6), then this 
latter equation may be rewritten as 


x sin kd¢x(s, y) 
ik 
— 164", d) + —sin cd — s 
k 


: [1 cos k(d — iat} cos Ky 
0 


“ {= f “f(b sin uly — Hat 


ik 
tel = (1 — cos oo x sin kd. (8) 
K 


By definition ¢,(s, y) is regular for ¢>k;. Hence, since 
x sin xd has zeros in o> —&i, the right-hand side of (6) 
must vanish for those zeros of x sin kd ino > —k,. These 
zeros will be at s=Kkm (m=0, 1, 2, -- - ) where 


m1? 1/2 

Km = ( 2 = #) ; 
Km iS positive when real, and positive imaginary when 
not real (for k real). It is also easy to show that Rkm <ks. 


Thus, in order that the right-hand side of (8) vanishes 
for s=Km, the following relations must hold 


d mrt 
$1'(m 2) = (—)"am ff 008 —* at eo 
0 


"eg . F} 
ga Rat ted f f(t. 
tk 0 
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If the expression for f(#) as a Fourier cosine series as ob- 
tained from the above expressions is substituted in (8), 
then application of the well-known results for term by 
term integration of Fourier series, together with some 
slight manipulation, yields the following equation for 


p+(s, y) 


d 
«sin kdo4(s, y) = — $4'(s, y) cos ky + te 
s+ ik 
mr y 
Ny, f cos 
K. 
+ — sin kd a 2 
: K Leg prs (9) 
where 
(—) bz! (Km, @) tk, d 
£m === sats Pe ae m ca 0, = o4! . x 
hon 21k 


We now consider the region y>d. The solution of (4) 
which satisfies the boundary condition on y=a is 


o(s, y) = B cos x(y — a) 


where B is a function of s only. Since 0V/dy vanishes on 
y=d for — 0 <x<0, ¢_’ (s, d) =0 and hence 


— Bx sin x(d — a) = ¢,4"(s, d). 
Thus we obtain the following es for #(s, y) 


cesar EL MELe a heh 
Asp ye pan (dd — 4! 
Hence 
cotk(d—a) | 
o+(s, d) sl o_(s, d) pe aie Ce tes iS: py iss d). (10) 
Substitution for ¢4(s, d) from (9) into (10) gives 
Gi eS eis.d 
_(s, d) = —————_ s,d) — - 
sis sin xcsin xd” s+ ik 
as Sy (=) 8m (11) 
d m=0 5: — Kaye 

On writing 

sin Ka 

K sin xc sin kd 

as 

N4(s) 

N_(s) 
it is seen that (11) may be rewritten as 
N_(s)¢_(s, d) 

N_(s) Caw ie 
= ’ cst eee = Ns 
Wa(9)4"(6, #) — Oneaanes 


The formulas for V,(s), N_(s) are given in the ne 
but a brief summary of their essential properties will be 
given here. N,(s) is regular and nonzero in the region 
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o>—k;; N_(s) is regular and nonzero for ¢o <k,. Also, as 
| s| 0, Ni(s)~s?2 and N_(s)~s1?, 

It is now convenient to rewrite the above equation 
for d_(s, d) as follows: 


N_(s) — N_(—ik) 2 <A) 8m 
N_(s)o_(s, d) + ro a — W-(s) yD ; = 2 
pros St WP ere Sie) a 
Caso (s ats Km) 
N_(—ik) 
= N(s)$+"(s, d) — neaeee® 

ge wine a eat one (12) 

ip ee (s ste Km) 


Clearly the left-hand side of (12) is regular for o<k; 
while the right-hand side is regular for o> —k;; hence 
both sides are equal to an integral function in the com- 
mon strip |o| <k;. This integral function may be de- 
termined from examination of the growth of both sides 
ot (12) as | s| —o., In order to do this we require the 
behavior of ¢,’(s, d) as | s| —o ing>-—k, and the edge 
condition shows that $,/(s, d)~s~/3 as |s| >. The 
integral function is now easily seen to be zero and we 
have 


N+(s)64"(s, d) 
Bee een 
s+ ik ih Sate an 


From (13), on substituting s=x,, the following set of 
equations is obtained to determine the $4' (km, d) 
1 1. 22 N (=k En 
2 kali Kn + tk d m=i Km(Km — Km) 
a,N_(—ik) 
Dikd(kn + ik) 


where d,=$4' (kn, @)/N_(—7k). 

For the case in which only one mode is propagated 
without attenuation in both regions of Fig. 1, it is pos- 
sible to reduce the above set of equations to a real set as 
follows. We define L,(s) = Ni(s) (s+7k), L_(s) = N_(s)/ 
(s—ik). From the formulas for N4(s) and N_(s) it is 
clear that L(x) is real for n>1. Hence (14) may be re- 
written as 


(13) 


(n= 0,1,2,---) (14) 


1 = OmL_—(— Km) 
n)Jan = _ a > ; ae ET mKkn — k? 
ni 3 ; (ee Keak Ken =F Kn) ; 
Lk , Woe hon ce a LR 
ee Sli dae ls 3 hes 
n=1 Km d 


From the asymptotic behavior of ¢+’(s, d) and N_(s) it 

is seen that the terms of the second series in (15) are 

O(m-!8) as m—  ; hence this series is convergent. 
Putting 7 =0 in (15), we obtain 


21k 2. Wrrdby —™ Km 
Be ey teal 85 ber aleted 


d m=1 Km 


__ agk_(— ik) 


mcly) 
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If the value for the series 


ie} 


i, ey Oe Km) / Km 


m=1 


as obtained from (16) is substituted into (15) then the 
following system of equations is obtained 


1 0 
Ly(kn) = x1! + aoLy(ik) — 7 Lin} 


1° 2 Onl — Kn) Coke — FB?) 
= 
p= 1, 2) eee 


Defining a new set of constants ¢, for n<1 by 


2dn 


Ch = 


Aeictegh ape i 130288) 


it is seen that these will satisfy the following set of equa- 
tions: 


1 2 CmL_(—Kkm)(Kmkn — k?) 
Ls(ka)én = 1 — — 
ats errant 


m=1 
(m = 1, 2,---). 


(18) 


The set of (18) is real only when the fundamental mode 
is propagated without attenuation. Finally substitution 
for d, into (16) yields the following equation for ao: 


L,(tk)ao — 1+ (ao/d)L_(—1k) 


= 9 19 
—Ly(ik)ao — 1 + (ao/d)L_(—ik) “ oe 
where 
g=— SD oot 
d m=1 Km 


CALCULATION OF REFLECTION COEFFICIENT 


In problems involving discontinuities in waveguides, 
the quantity which is of physical importance is the re- 
flection coefficient. In this section, we calculate this for 
the fundamental mode when ka/z <1. 

The reflection coefficient is the residue of ¢,(s, y) at 
s= —?k and, from (9), it is given by 


We have from (13) 
$+ (—ik, d) 
ie 1 {nia it si N_(— km) b+" (Km o 
s>—ik Ni(s) | s+ ik [ORE Km(S ++ Km) 


and 
—a N_(—ik) 
2ikde s+ik 


Nils) ~ as s—— ik, 


1957 


Thus the reflection coefficient R is given by 
d (ik, d) 
eel t hi! : 
a ( Oe 2ikd ) 
Substitution in (19) for 6,’(ik, d) gives 
aR Ls (ik) 


ad Fae) S if ia 


1fi= . (20) 


SSG 


For simplicity it is now convenient to introduce the 
nondimensional expressions H,(s), H_(s) defined by 


a | tf 
ies / ~1.6),  H() = y/ = 1.6). 


If this change of notation is made in (20), then, writing 
a for c/a, we obtain 


R H?(ik) dj H?(ik) 
(1 — a) {1 ‘a a \ a 

R H?(ik) H.?(ik) 
(1 — a) {i a \ Fs a 


Eq. (21) may be inverted to give 


ae aS ct one) , 22) 
a+ H,(ik) — igla — Hy?(ik)) 


where 6 is written for 1—a. The expressions in the Ap- 
pendix clearly indicate that | H,.(ck)| =1 and hence 
that H,(ik) may be written as e* (with @ real). 
; It is instructive to consider the physical significance 
of the expression defined by (22) when g=0. This ex- 
pression may be derived from physical considerations, 
as will now be shown. We assume that the step is absent, 
and superimpose a field Ae~** incident from the region 
0<y<d, x<0 on the original field. This incident field 
will cause a wave A R,et+*** to be reflected back in the 
negative x direction (Ri may be calculated by use of 
the Wiener-Hopf technique) and we thus find that the 
total velocity component of the potential in the funda- 
mental mode for y<d, x <0 is 


ARy,e*** + etka oh Ae ikz 


ag.. (21) 


and we choose A so that the partial derivative of this 
expression with respect to x vanishes on the line «=0, 
0<y<d. If this is done, then the transmission coefficient 
for the problem of a wave Ae~‘*** incident in the region 
formerly occupied by the step is seen to agree with the 
expression of (22) for g=0. Thus we see that the ex- 
pression (22) illustrates the connection which was 
alluded to in the first section between the bifurcation 
and the step problems. The assumption g=0 of course 
merely implies that V is constant over the short vertical 
part of the step. For ka <r, it is possible to compare the 
present results with those given in the “Waveguide 
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Handbook,” and, since the solution given there is an 
equivalent circuit formulation, the present results have 
to be recast in this form. The equivalent circuit is 


Yo Y,’ 


where Yo, Yo’ are the characteristic admittances of the 
two guides and B is purely real. B/Y,'=C may be ex- 
pressed in terms of the reflection coefficient R, and we 
have 


{(a+ 1)R+(a—1)} 
a(R — 1) 


iC = — (23) 


If the expression for R in terms of g is inserted in (23) 
then we obtain 


6 (sind cos 0 
Cu oe 


24 
cos @ — gsin@ oe 


Inspection of this expression for C shows that it may be 
obtained by applying the above simple argument for 
determining the reflection coefficient to the case in 
which the bifurcation extends a distance tang to the 
right of the origin, and we then obtain the exact form of 


_C. Thus the reflection coefficient for the step problem is 


the same as the transmission coefficient for the bifurca- 
tion problem mentioned above, when the bifurcation 
extends a distance tan! g further towards the positive 
x direction than the original step structure. 


COMPARISON WITH THE RESULTS OF THE 
“\WWAVEGUIDE HANDBOOK” 


Some simplification will be afforded if (18) is rewritten 
in dimensionless form by making the substitutions 
given in the preceding section and defining a new set of 
constants d, as c,d~!/2, The following equation is then 
satisfied by the d,: 


Am(Hmbin — €°) al ayes a (25) 


Sy 


m=1 im bin => a 


H4(kn)dn <a qil2 


where e=kd/m and pp2=m?—e. The function g may 
now be written in dimensionless form as efa!/? where 


f= ys mH _(— Km) / bm: 


The formulas are now in the simplest form possible and 
are suitable for comparison with the results of the 
“Waveguide Handbook.” First of all we shall consider 
the special case of neglecting all powers of e>2. To this 
order of approximation it may be shown that for m>0 


A™m!(m(1 — 6)/6)! 
(m/6)! 


Hy(km) = 
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where 


5 \-a)/8 
(eat) 


61/6 


The argument of H,(ik) may be shown to be ¢ log A, 
and we thus have 


C 5 
— = — (log A + a/*f), 
€ a 

The amount of computation involved in this special 
case is fairly small, since to the present order of approxi- 
mation the wm’s reduce to m. Eq. (25) has been solved 
on the assumption that all the coefficients, apart from 
the first five, vanished. This gives for 6=} the result 
C/e («=0) =1.556, while the corresponding result in the 
“Waveguide Handbook” is 1.364. For 6 =0.8 the present 
method gives C/(ka/m) =2.45 while the value in the 
“Waveguide Handbook” is 2.48. 

In order to test further the accuracy of the present 
approach, the general expressions developed in the 
Appendix were used to calculate H_(—x,,) for certain 
values of e, and an approximate solution of (25) was 
obtained. The calculations were carried out for e=4, 
6=4; «=0.64, 6=0.8; and the values obtained for C 
were 0.954 and 3.67. The values given in the “Wave- 
guide Handbook” are 0.96 and 3.73. The values of ¢€ and 
6 were chosen so that ka/z ranges from 0 to 0.8; hence 
if there were any large deviation between our results 
and those in the “Waveguide Handbook,” it would have 
been evident in the numerical results. It is thus evident 
that there is very good agreement between both meth- 
ods, and that the solution of the first five equations is 
sufficient to give a very accurate answer. In fact, com- 
plete neglect of g produces a percentage error of less 
than 30 in the expression for C while the percentage 
error in R will of course be considerably less (2.e., a 
very large error in determining g does not produce an 
equally large error in the expressions for R and C). 
This indicates that the assumption of a constant value 
for VW on the vertical step is a fairly satisfactory one. 


THE CASE oF ka/r>1 


We now consider the modifications which have to be 
made to the above analysis when ka/r>1 butkc/a kd/r 
<1. One difference is that there will be a term in V which 
behaves as exp (—7inx), where 


n= (8 t/a)”, 


as x— 0 and 7 is now real. The coefficient of exp (—inz) 
may be found, after some manipulation, to be 


cos Ty 2V/ 26/08 
a (6B + e) 


f,-— ; cote F Row 
6(1 + tefol/?) 21) bm — 18 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


April 
where €m =dml_(—Km) dm, B = (2 — 6?)!/?, and 


~~ {56,40 + 5\(™, «,0) 


A id 
-si(4, 1,0) — Si(e, 0,0) — 3,(S—, 0, 0) 


0, 0) — fe z 8) (1+ aloga + 6 logs) 


and 


— ; x x 
Sv(4%, a, B) = {sio- Ses hie 
nb eae) n=N SE pe ne n 
The Sy’s are tabulated in the “Waveguide Handbook.” 
The above expressions are now in the most convenient 
dimensionless form for numerical computation. 

Another different feature of the case ka/r>1 is that 
the functions H_(—x,) are no longer real and we have 


1 
Hs(km) a re 2 1B) G+(km) 


where the G,’s are real. The equation for the new un- 
knowns eé, defined above is 
2. m(umiin — €)? 
UnénEly2(kn) = al? — @ py pikes e th A ATL) 
(Lim or Hn) 


Eq. (26) may be written as two sets of real equations as 
follows: 


* (26) 


m=1 


Mnfn bn *gnG?(Kn) 
G5" (Ka) tnt rp ee Oper re ae 
ad? ad? 


ss (Banker ia e”) 
3 pes > PB Ra oh 
as (Um Hn) 
G42(Kn) Hn 2B ptn°G2(Ien 
4°( Kn) Me pray? Bun?G?(kn) vi 


ad? ad? 


© 


: (Umbin — €”) 
m=1 (im =+ Ln) 


_— 


(28) 


where én =fn+tgn- 

It is possible to reduce (27) and (28) to two sets of 
linearly-independent equations, but since in the present 
case they were solved on the digital computer, there 
does not appear to be any advantage’ gained by their 
reductions. The above equations also have the advan- 
tage that they are in a form from which it is fairly simple 
to estimate the error involved in their solution. It is 
possible to estimate the error made in solving only a 
finite number of the above equations. The analysis, © 
which requires the establishing of various inequalities 
for the functions occurring in the equations, will not be 
presented here, but it is available in the author’s report.! 


‘W. E. Williams, “Step Discontinuities in Waveguides,” Res. 
Rep. No. EM-77, Inst. Math. Sci., New York University, New York, 
N. Y., April, 1955. 
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NUMERICAL RESULTS 


The set of equations occurring for the case ka/7>1 
was solved by assuming that fr, g.=0 for n>5. Since 
there is no particular physical significance in the actual 
numerical values of the coefficients f, and g,, they have 
not been tabulated. One quantity which is of particular 
interest is the reflected energy crossing a plane normal 
to the x axis which will be proportional to 


ae 


lim w* — dy, 
Ox ‘ 


z~0 J 9 
where the asterisk indicates the complex conjugate. If 
the coefficient of e-*", referred to in the foregoing sec- 
tion, is written as R; cos ry/a, then the above integral 
may be evaluated to give 


~ia{t| RE +49] Rf. 


Hence the (reflected energy/incident energy) ratio is 
B 
[| +— | Rl 
2e 
The two quantities 
[ap [RP+ = [R= [RP 


are given in Table I. The numerical values show that 
the ratio of the reflected energy to the incident energy 
increases, for constant step height, as the waveguide 
height decreases. This, of course, was to be expected on 
physical grounds. The percentage of the reflected 
energy which is due to the | R,|? tends to increase as the 
waveguide height decreases for constant step height. 
This percentage is, however, fairly constant, ranging 
between 19 and 25 per cent for all the values given in 
Table I. 
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APPENDIX 


In this appendix, formulas are developed which may 
be used for computing the functions N;(km). If the same 
method is employed as that used by Carlson and Heins® 
then it is fairly easily shown that 


Wye treee lhe a P(8) 

TOS AB) Vrtie 6 
P(——) pay 
1-6 
s{ aloga —clogc — dlogd} 
-exp ——___________—>_ 
Tv 
N_(— s) N,(s) = — 1, where 


2 e \ 1/2 SE SE 
PO) TE a) ee 

It is not really necessary to retain the exponentials in 
the above infinite products when the complete expres- 
sion for NV4(s) is being considered, but they must be re- 
tained for evaluating each separate product. For gener- 
ality we shall consider the evaluation of each of the 
products separately since products of a similar nature 
occur in the other problems of diffraction theory and 
only slight modification of the following analysis will be 
required to cover other cases. 

It is clear from the expressions occurring in the text 
that we must compute expressions of the form 


© B2q2 1/2 y ie 

EA eines bere 

n=l nx? nT nw 
Though most of the computations involved in this work 
were carried out on a digital computer it is also possible 
to obtain fairly convergent expansions in powers of & 


for the above products. These expansions are conver- 
gent for all the powers of € and 6 that occur, and have the 


TABLE I 
e=0.9 e=0.8 e=0.7 e=0.6 
ka/«w 1.8 1.6 1.4 lie 1.6 1.4 12 1.6 1.4 122 11 15 1.4 
R|? 0.25 0.32 0.41 0.23 0.26 0.3 0.47 0.20 0.28 0.39 0.49 0.17 0.20 
R2|? 0.31 0.40 0.52 0.29 0.33 0.4 0.58 0.25 0.35 0.31 0.63 0.21 0.26 
ee ee ee 2 Se SS a ee 


The values obtained for | R|? were compared with 
those obtained by assuming that all the coefficients of 
the set of equations vanish, and the two results were 
found to differ by only a few per cent in all cases. The 
difference between the two values obtained for | R:|? 
was greater, but this contributes only a small amount 
to the total energy. It thus appears that a very good 
estimate may be obtained for the reflected energy by 
assuming a constant value of the magnetic field on the 
step. In all cases the values obtained for f; and g; did 
not differ by more than 5 per cent from those obtained 
by assuming that f, and g, were the only nonvanishing 
coefficients. 


advantage of being very much more suitable for numeri- 
cal computation than straightforward multiplication of 
the terms of the product. In order to obtain the numeri- 
cal values for a large number of values of ¢ and 6, the 
digital computer has a distinct advantage over any other 
form of computer. The program developed for the pres- 
ent problem has also been used for computing infinite 
products for which the approach used here may not be 
used (in these products the mm in the infinite product 
was replaced by the nth zero of a Bessel function or of a 
transcendental equation involving Bessel and Hankel 


5 F. Carlson and A. E. Heins, “Diffraction by a pair of parallel 
plates,” Quart. Appl. Math., vol. 4, pp. 313; 1947. 
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functions). The following formulas are suitable when not 
too many values of « and 6 are required, and the com- 
putations for ka/x <1 were carried out using these for- 
mulas. Expansion of 

k?q? 1/2 

=) 


and km in powers of ¢, together with some algebraic 
manipulation, gives the following result for m0 


Cs) ka? 1/2 Kava ena 
0 soar asec feo an 
et na nT nr 


mi) ~* 
— p-(imis) {= exp P,,(6) 


where x}, is the factorial function, j is Euler’s constant 
and 


Pm(6) = ayet + aoe® + aze® + aye? + age}? 


where 
8md2a, = — £(3), 48m3d°ag = — 3m?E(5) — 26°¢(3) 
2*m®57a3 = — Sm*¢(7) — 3m6?¢(5) — 364¢(3) 
5-29m75%a, = — 35m®>(9) — 140m462¢(7) + 36m754e(5) 
— 208%¢(3) 


3-2 m% "a, 


I 


— 63m%¢(11) — 35m'S?¢(9) — 31m454¢(7) 
— 30m76*°¢(5) — 3568¢(3) 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


April 
where {(s) is the Riemann Zeta function, 1.e., 


21 
s(s) = DU —. 
n=1 N 

The other two products in N;(km) are of the same form 
as the above with a replaced by ¢ and d respectively, 
and hence the above formulas can be extended to cover 
the other two products by replacing 6 in the expressions 
for a, ---, a5 by 6/(1—4) and 1 respectively. These re- 
sults are valid for (e/6) <1 but may be extended to in- 
clude the case of 1 <(e/6) <2, €<1 as follows. 

For (e/65)>1 the first term in the above infinite 
product will be complex, and hence if we consider the 
integer m as now ranging from 2 to infinity this new 
product may be treated exactly as before. The ¢(s) in 
the formulas for a; - - - , ds should now be replaced by 
f(s). 

The case of m=0 requires special consideration. In 
this case the products occurring in N,(ik) may be writ- 
ten in the form 


= ka ka 
exp if > (sin — =) : 
1 nT nT 


The functions occurring in this exponent have been ex- 
tensively tabulated in the “Waveguide Handbook.” 
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The Transient Behavior of the Electromagnetic 
Ground Wave on a Spherical Earth* 
JAMES R. WAIT 


Summary—Some calculations are presented to show the nature 
of the transient ground wave radiated from an electric dipole which 
is situated over a spherical earth. The moment of the dipole is con- 
sidered to vary with time in a linear manner. It is shown that the de- 
parture of the leading edge of the radiation field from a step function 
form is a consequence of diffraction and loss in the finitely conducting 
ground. 


* Manuscript received by the PGAP, February 17, 1956; revised 
manuscript received, October 26, 1956. 
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INTRODUCTION 


T IS THE PURPOSE of the present paper to calcu- 
if late the transient response of the ground wave due - 

to a suddenly applied current at a source located over 
a spherical earth. With no subsequent loss of generality, 
the source can be considered to be an infinitesimal cur- 
rent element whose moment varies with time as a ramp 
function. Furthermore, the height of the vertical current 
source, whether it be a lightning channel or a man-made 
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transient electromagnetic source, can be synthesized by 
a line distribution of dipole sources. Since the observer, 
in most practical instances, is at a distance large com- 
pared to the height of the column, the radiated field is 
proportional to the instantaneous value of the surge 
current integrated over the length of the column. 

The extension of the responses calculated for the ramp 
source to the case of more complicated source functions 
such as lightning strokes can be calculated by super- 
position. As mentioned in the preceding paragraph, the 
analysis will treat only the omnipresent ground-wave 
pulse. The transient signal propagated via the iono- 
sphere, albeit is the most important component at large 
distances, is not considered here. 


THE STEADY-STATE RESPONSE 


It is desirable to first illustrate the radiation from a 
harmonic source at very low radio frequencies. With 
respect to a spherical coordinate system (7, 0, ¢) the 
earth’s surface is defined by r=a and the source is a 
radially oriented dipole of moment P(w) with a time 
factor exp (iwt) located at 9=0 and r=a. The earth is 
assumed to be homogeneous with conductivity o and 
dielectric constant e. It is convenient to express the field 
components in terms of a Hertz vector which has only 
a radial component U as follows: 


a or? a iG 
diera2 r 
Fg (<u) 
r orod \a 
in eng ek ih (1) 
— ed 4 —_— 
i . aoe 


where 6=27/free space wavelength and € is the di- 
electric constant of free space. Van der Pol and Brem- 
mer? have shown that to a good approximation, U, for 
ra, can be expressed in the following form 


U = 2U.W (2) 
where U> is the Hertz potential of the primary field and 
W is the “Watson-type residue series” given by 

Ee enineX 
W = (2m) 1/2X e-i(a/4+9/2) yy eters (3) 
where X = (6a)!/30 and 
a B17/B° 
~ (Ba) [(6.2/6%) — 1)" 


1R. B. Morrison, “The variation with distance in the range 
0-100 km of atmospheric waveforms,” Phil. Mag., vol. 44, pp. 980- 
986; June, 1953. : ; 

2B. Van der Pol and H. Bremmer, “The diffraction of electro- 
magnetic waves from an electrical point source round a finitely con- 
ducting sphere, with applications to radio telegraphy and the theory 
of the rainbow,” Phil. Mag., vol. 24, pp. 141-175, July, 1937; pp. 
825-864, November, 1937; vol. 25, pp. 817-834, June 1937; vol. 27, 
pp. 261-275, March, 1939. 
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and 4=42 X10-7. The coefficients 7, are related to the 
roots of an equation involving Hankel functions of order 
one third. They have been treated numerically in an 
extensive manner by Bremmer.’ Since 6 is not a large 
number, the primary Hertz potential is given by 


P e—iBD 


ae (4) 
Ariwes D 


where D=aé is the distance from the source to the 
observer. Remembering now that £, is a solution of the 
wave equation, it can be expressed as 


e--teiQ] ow 


r sind 


so that the vertical field Z on the surface of the earth is 


Se) (6) 


subject.again to 6 being small. The fact is now utilized 
that W is a slowly varying function of D compared to 
exp (—7i8D)/D, so that 


1pP 4 1 
yi ree E = ——— | 

2xD BD BD? 
For a highly conducting earth, W is near unity at short 
distances and the terms in powers of 1/D, 1/D?, and 
1/D* can be identified as the radiation, induction, and 
static fields, respectively. It is to be noted that W isa 
correction factor which accounts for finite ground con- 
ductivity and diffraction. It modifies not only the radia- 
tion field, but also the induction and static fields at very 
low radio frequencies. 

Using data computed from (7), the amplitude and 
phase of the field Z are plotted as a function of frequency 
in kilocycles per second in Figs. 1 and 2. The conduc- 
tivity of the path is that of sea water ( = 4 mhos/meter) 
and the effective earth’s radius is taken to be four thirds 
that of the actual earth radius, to account for normal 
atmospheric refraction. The amplitude of the fields is 
normalized so that the radiation field at short distances 
is the same. In other words, the quantity that is plotted 
is 2tDE/iuwP. The curves show clearly the effect of the 
induction and static fields at lower frequencies and 
shorter distances. At greater ranges, the relative ampli- 
tude drops appreciably as a result of diffraction by the 
earth’s curvature and the phase lag increases consider- 


eB, (7) 


3H. Bremmer, “Terrestrial Radio Waves,” Elsevier Publishing 
Co., pp. 105-106; 1949. 

4]. R. Wait and H. H. Howe, “Amplitude and Phase Curves for 
Ground-Wave Propagation in the Band 200 cps to 500 Kc,” Natl. 
Bur. of Standards Circular 574, Washington, D. C.; May 21, 1956. 
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Fig. 1—The amplitude of the ground wave as a function of 
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Fig. 2—The phase of the ground wave as a function of frequency. 


ably. If the earth were flat and perfectly conducting, so 
that W=1, the relative amplitude curves would ap- 
proach unity and the phase would approach zero at the 
higher frequencies, 


THE TRANSIENT RESPONSE 


When the moment of the source varies as a function of 
time, the instantaneous response can be expressed, at 
least formally, by a Fourier integral, if the appropriate 
frequency spectrum for the source is employed. For 
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example, if the product of the instantaneous height of 
the column times the instantaneous current is p(¢), then 
the frequency spectrum P(iw) of the source is given by 


Pin) = [pleted (8) 


0 


it being assumed that p(t)=0 for <0. The field re- 
sponse e(t) is then expressed by 


1 +io 7 iy ; 
di) = 5 ee E(ies)e*'d( ies) (9) 
where 
- 
Tet) lg 4 te 
2nD oD ' (iD)? 


and c= (e)—/2=3 X108 meters per sec. The expression 
W(iw) is identical to the W defined in (3) where it is 
emphasized here that it is a function of iw. The principle 
task is to evaluate the contour integral for e(t). The 
difficulty of this task stems from the complexity of the 
function W(iw) which is in the form of a residues series 
whose individual terms themselves are rather compli- 
cated functions of zw. In a previous paper® a method for 
evaluating the contour by operational methods was 
described which was suitable for shorter ranges where 
the earth could be regarded as nearly flat. The formula- 
tion developed there regarded the diffraction by the 
earth’s curvature as a correction to the flat-earth 
Sommerfeld formula, and consequently, W(iw) could be 
expressed in a more convenient form. 

At large distances, where diffraction plays a major 
role in the analytic behavior of W(iw), it seems neces- 
sary to evaluate the contour integral for e(¢) by numeri- 
cal means. A method for doing this has been described 
by the author,® which is particularly suitable if the real 
part of E(iw), when plotted on a log w base, can be ap- 
proximated by a series of straight line segments. To 
apply this method, the function E(iw) is written 


M 1 
E(ia) = > An —— + Rw) + iX(@) 


m=1 (iw) = 


(11) 


where &M is a finite integer and where R(w) and X(w) are 
real and have no poles on the real axis of w. The first 
term contains the singularities of E(iw) on the imaginary 
axis of tw (2.e., on the real axis of w) that can be intro- 
duced by the frequency spectrum P(iw) of the source. 
The coefficients A,, are the residues of these poles. Sub- 
ject to the restriction that there are no poles of E(iw) in 
the right-hand part of the complex plane of iw, it follows 
that 


° J. R. Wait, “Transient fields of a vertical dipole over a homo- 
geneous (gently) curved ground,” Can. Jour. Phys., vol. 34, pp. 27- 
35; January, 1956. 

6 J. R. Wait, “An approximate method of obtaining the transient 
response from the frequency response,” Can. Jour. Tech., vol. 31, 
pp. 127-131; March, 1953. 
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where 
, @ ve 
e(t) = — i} R(w) cos wid. 
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As shown in the earlier paper,’ e’(f) can be approxi- 
mated by 


2 N 
e'(t) ~ — DS Ba[Si(wnat) — Si(wat) Ji? (13) 
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where 


eae f eS 
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is the sine integral and B, is the slope of the segment of 
the R(w) curve between the points w,_; and w, and is 


defined by 
a R(@n) — R(wr_1) 


14 
log (wn/@n—1) tf) 


In applying the above formula it is convenient to take 
@)=0 and w, the angular frequency where R(w) departs 
from its static value so as to make B,~0. On the other 
hand, wy is taken at a sufficiently large value so as to 
have R(wy)-0. The number of intervals is chosen so as 
to obtain the degree of accuracy required. In the present 
instance, JV is taken to be 10. 

To illustrate the fundamental character of the tran- 
sient response, it is considered desirable to choose a 
source that would give rise to a radiation field on a flat 
perfectly conducting earth in the form of step function. 
This is effected by taking the moment p(#) to be a 
linearly rising current or a “ramp” function at ¢=0. 


That is 


p(t) — Po(t/to) fort >0 (15) 


==) tors =< 0: 


The frequency spectrum of the source is then given by 


P(iw) = Po(1/iw)?/to. (16) 
Using the steady frequency response data* of W(iw) 
calculated for frequencies from 200 cps to 500 kc, the 
approximate method described above is used to calcu- 
late the transient response e(#) to the ramp function 
source. For purpose of presentation, the response is 
written 


Popo 
i) = — A(t’ 
i) 2rDi ©) 


(17) 


where t’=—D/c. On a flat perfectly conducting earth, 
W(iw) =1, and consequently 
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That is, for short times (t/&D/c) the response is charac- 
terized by a step function. For longer times, the terms 
proportional to ¢’ and (¢’),? which are the induction and 
static fields, respectively, become significant. The re- 
sponse for ground conductivities of 4, 10-2, and 10-3 
mhos per meter are shown in Figs. 3, 4, and 5, respec- 
tively. The departure of the leading edge of the wave- 
form from a step function is a consequence of diffraction 
and the losses due to the finite ground conductivity. It 
can be seen that for large distances of the order of 1000 
miles, the radiation component of the ground-wave 
pulse takes the order of 100 microsecs to build up to 
one half its final value. Since the induction and static 
fields are present, the response function A(é’) actually 
rises above unity although this is not noticeable at larger 
distances for times less than one millisecond. The rise 
time at short distances in the range from 50 to 100 miles 
is very rapid. At these distances the effect of ground 
conductivity is more significant and diffraction plays a 
negligible role. In fact, at a distance of 75 miles the re- 
sponse for an all sea water path (¢~4) has risen to about 
0.95 at 1.0 microsecond. For the all land paths with con- 
ductivities of 10-? and 10-* the corresponding responses 
at 1.0 microsecond are 0.8 and 0.2 respectively. 

The response function A(t’) for the conductivities of 
10-? and 10-* and for distances of 37 Mi and 75 Mi 
which are included in Figs. 4 and 5 can be compared 
directly with those calculated in an earlier paper® by an 
altogether different method. In the previous work, the 
integrations could be carried out by analytic means for 
these short distances. The agreement between the two 
sets of results is good, 


CONCLUSION 


While the transient responses were calculated ex- 
plicitly for a ramp function source, it is not difficult to 
extend the results to more complicated sources. The 
response 4@(t) to a source function p(t) can be written in 
terms of e(t) the response to the ramp function fot/to as 
follows 


ai) = & fienet — t)dr (18) 


where 
1 +i00 
a if (ius)? P(ieo)e'd (i). 
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In general, the convolution integral representation for 
a(¢) must be evaluated by numerical means. A particular 
type of source that is convenient to handle by this 
method is the “tent” pulse. It may be represented in the 
following way: 
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Fig. 3—The transient response of the vertical electric field at various 


distances over sea water (c =4) from a source whose moment varies 
as a “ramp” function. 
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where u(t)=1 for ¢>0, =0 for ¢<0. In this case the 
moment of the source rises in a linear manner from 0 at 
t=0 to fo at ft, and then decays in a linear manner to 0 
at t=t.. The response @(¢) is then given by 


be 


te — ty 


a) = {el — ett — 4) 


by 
+ e(t — te) ; (20) 


le — ty 


where it is understood that e(t;) =0 for ¢;<0. 

It is believed that the results in this paper shed some 
light on the transient behavior of the ground wave. It is 
hoped that the present analysis can be extended to in- 
clude the reflections from the ionosphere which would 
of course mask the ground-wave pulse at longer times.7 


7J. R. Wait and C. Froese, “Reflection of a transient electro-: 


magnetic wave at a conducting surface,” J. Geophys. Res., vol. 60, 
pp. 97-103; January, 1955. 
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Fig. 4—The transient response of the vertical electric field at various 
distances over land (« =10~) from a source whose moment varies 
as a “ramp” function. 
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Fig. 5—The transient response of the vertical electric field at various 
distances over land (¢ =10~’) from a source whose moment varies 
as a “ramp” function. 
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An Experimental Investigation of the Diffraction of 
Electromagnetic Waves by a Dominating Ridge* 
J. H. CRYSDALE}, J. W. B. DAY}, W. S. COOK{, M. E. PSUTKAf, anp P. E. ROBILLARD+ 


Summary—The diffraction of electromagnetic waves by a domi- 
nating tree-covered ridge north of Ottawa has been the subject of 
an experimental investigation. The measurements were made at 
173, 493, and 1785 mc using both horizontally and vertically polarized 
radiation. The measurement period extended from late in the winter 
until early in the summer of 1955. The bulk of the experimental 

_ data consists of the results of height-gain and azimuth-gain measure- 
ments. Estimates of diffraction loss have been deduced from these 
results for the component of the radiation which was not reflected by 
the terrain between the ridge and the terminals of the experimental 
circuits. These deduced experimental diffraction losses exceed the 
predictions of Fresnel-Kirchhoff knife-edge diffraction theory, the 
discrepancies increasing with increasing frequency. There was no 
evidence of either a pronounced polarization dependence or a pro- 
nounced seasonal dependence. 


INTRODUCTION 
ap HE EXPANDING use of nonoptical communica- 


tion circuits at vhf and uhf has revived interest 
in the diffraction of electromagnetic waves by 
terrain irregularities such as hills, ridges, and moun- 
tains. In general, the geometrical and electrical charac- 
teristics of obstructions of this type preclude rigorous 
analyses and predictions of diffraction loss are often 
based on the Fresnel-Kirchhoff diffraction theory.! While 
this theory has the advantage of relative simplicity, it 
does not take into account such significant parameters 
as the radii of curvature at diffracting edges, the compo- 
sition of the obstruction, and the polarization of the ra- 
diation. Nevertheless, it has been concluded for several 
of the circuits previously investigated that the diffrac- 
tion losses deduced from experimental results are in 
good agreement with the values predicted by the Fres- 
nel-Kirchhoff theory. 
The purpose of this paper is to describe and present 
the results of an experimental investigation of the dif- 


* Manuscript received by the PGAP, August 31, 1956; revised 
manuscript received, December 19, 1956. This paper is based on 
RPL project rep. 34-0-3 (dated December 30, 1955) and represents 
work carried out under project PCC No. D48-38-03-21. Some of the 
results were presented in papers delivered at the 1955 Fall URSI 
meeting, Gainesville, Fla., December 15-17, 1955, and at the Cana- 
dian IRE Cony., Toronto, Ontario, October 1-3, 1956. 

+ Radio Physics Lab., Defence Res. Board, Ottawa, Canada. 

t Federal Electric Corp., Streator, Ill., Formerly at the Radio 
Physics Lab., Defence Res. Board, Ottawa, Canada. ’ 

1 A fairly complete bibliography of the previous theoretical and 
experimental work is contained in the following recent papers. 

R. S. Kirby, H. T. Dougherty, and P. L. McQuate, “Obstacle 
gain measurements over Pike’s Peak at 60 to 1046 mc,” Proc. IRE, 
vol. 43, pp. 1467-1472; October, 1955. (It should be pointed out that 
sections of references 1, 7, and 8 of this paper were criticized by J. H. 
Crysdale, Proc. IRE, vol. 43, p. 627; May 1955. See also Proc. IRE, 
vol. 43, p. 874; July, 1955.) 

K. A. Norton, P. L. Rice, and L. E. Vogler, “The use of angular 
distance in estimating transmission loss and fading range for propa- 
gation through a turbulent atmosphere over irregular terrain,” PRoc. 
IRE, vol. 43, pp. 1488-1526; October, 1955. 

Some later work is presented in the paper by G. W. Swenson, Jr., 
“VHF diffraction by mountains of the Alaska range,” Proc. IRE, 
vol. 44, pp. 1049-1050; August, 1956. 
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fraction of electromagnetic waves by a dominating ridge 
north of Ottawa. This ridge constitutes a portion of the 
southern boundary of the Precambrian Shield and is 
heavily wooded, the trees being primarily of the 
deciduous type. 


EXPERIMENTAL MrtTHop 


The bulk of the experimental data consists of the re- 
sults of height-gain and azimuth-gain measurements on 
two obstructed circuits traversing the ridge and on the 
line-of-sight circuits between a site on the top of the 
ridge and the terminals of the obstructed circuits. These 
measurements were made at three widely separated fre- 
quencies in the vhf and uhf bands (173, 493, and 1785 
mc), using both horizontally and vertically polarized 
radiation. In order to investigate the possibility of sea- 
sonal effects, the measurement period extended from late 
in the winter until early in the summer of 1955. At the 
beginning of the period, the terrain was covered with 
snow. At the end of the period, the foliage was well 
advanced toward maturity. 

The locations of the experimental sites and the 
salient topographical features are indicated in Fig. 1. 
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Fig. 1—Locations of the experimental sites. 


One site was located to the south of the ridge at the 
Radio Physics Laboratory (RPL). Two sites were 
located to the north of the ridge, one near the village of 
Farm Point and the other near the village of Farrellton. 
The site on the top of the ridge was visible from the 
other three sites. Due to the restrictions imposed by the 
nature of the terrain and the availability of roads, it was 
not possible to find suitable sites which would lie 
exactly in a common vertical plane. 

The vertical profiles of the obstructed circuits are 
presented in Figs. 2 and 3.? The terrain to the south of 


2 Figs. 2 and 3 are based on the 31 G/5 West Half and 31 G/12 
West Half 1:50,000 topographical sheets published by the Army 
Survey Establishment, Dept. of Natl. Defence, Ottawa, Canada. 
These profiles have been drawn for an effective earth radius equal 
to 4/3 true earth radius in order to compensate for tropospheric 
refraction. 
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Fig. 2—Vertical profile of the RPL-Farm Point circuit. 
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Fig. 3—Vertical profile of the RPL-Farrellton circuit. 


the ridge is mainly gently rolling farm land. In contrast, 
the terrain to the north is more rugged and relatively 
heavily wooded. 

Antenna height and azimuth angle were continuously 
variable at the sites to the south and north of the ridge 
due to the use of path testing towers.? The RPL tower 
was located on the roof of the laboratory, where its base 
was approximately 110 feet above the level of the 
Ottawa River. At Farm Point and Farrellton, the bases 
of the towers were located at ground level. Height and 
azimuth angle information was transmitted auto- 
matically and recorded simultaneously with the re- 
ceived signal levels. The antennas at the Ridge Site were 
mounted in fixed positions on a 50-foot scaffold tower. 

Directional antennas were used at each of the experi- 
mental frequencies. Yagis, with gains of 10 db relative 
to an isotropic radiator, were used exclusively at 173 
mc. At 493 mc, RCA Victor directional antennas (15 db) 
and a corner reflector (10 db) were employed. The 
antennas for the measurements at 1785 mc consisted of 
dipole-disk fed paraboloidal reflectors (24 db) and 
pyramidal horns (15 and 21 db). Except where stated 
otherwise, each antenna was mounted for horizontal 
beaming of the maximum in the main lobe of the free 
space antenna pattern. The vswr’s at the terminals of 
the antennas were maintained at less than 1.1. Solid 
dielectric cables (RG-8 and RG-17) were used as trans- 
mission lines at the two lower frequencies. At 1785 mc, 
where the attenuation of solid dielectric cables was ex- 
cessive, 2 inch diameter Heliax cable was employed. 

The transmitters for the measurements on the ob- 
structed circuits and on the optical RPL-Ridge circuit 
were located at RPL. These transmitters and the re- 
ceivers at the other sites were modified components of 
commercial communications systems. Signal generators 


3 The towers used were the type described in a paper by R. D. 
Campbell, “Path testing for microwave radio routes,” Elec. Eng., vol. 
72, pp. 571-577; July, 1953. 

4 Manufactured by the Andrew Corp., Chicago, IIl. 
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at the Ridge Site served as transmitters on the Ridge- 
Farm Point and Ridge-Farrellton optical circuits. All 
transmitters were operated unmodulated cw in order to 
facilitate power measurement and calibration. Power 
levels were monitored continuously at the transmitter 
output terminals. In addition, frequent measurements 
of absolute power were made at the output terminals of 
the transmission lines to the antennas. The receiving 
systems were calibrated at the input terminals to their 
transmission lines. All power meters and bridges were 
calibrated against one bridge which served as a 
standard. 


EXPERIMENTAL RESULTS 


Experimental results were obtained on the ob- 
structed RPL-Farm Point circuit and its associated 
optical circuits from Jate in March until May 24, and 
on the obstructed RPL-Farrellton circuit and its asso- 
ciated optical circuits from May 25 until the end of 
June. 

A number of experimental curves which illustrate 
some of the features of the results will now be presented 
and discussed. In general, the curves shown have been 
plotted in terms of a parameter A which represents the 
difference in decibels between experimental transmis- 
sion loss and predicted free space transmission loss, 
assuming the experimental antenna gains for the calcu- 
lation of the latter quantity. The sign of A is defined to 
be such that positive values correspond to signal levels 
in excess of the predicted free space levels. 

In the results which follow, each height-gain curve 
represents the average of at least one run made by 
raising the antenna and one made by lowering the an- 
tenna. Similarly, each azimuth-gain curve is the average 
of at least one clockwise and one anticlockwise measure- 
ment. 


Height-Gain Measurements 


The results of some of the height-gain measurements 
at 493 mc are summarized in Figs. 4 and 5. The curves 
in Fig. 4 represent the results of measurements in which 
antenna height was varied at RPL while the signal level 
was recorded simultaneously at the Ridge and Farm 
Point Sites. Each of these pairs of curves have been 
superimposed for best fit in the region of minimum 
transmission loss. “Antenna height” refers to height 
above the base of the appropriate path testing tower. 
Some experimental curves obtained by varying antenna 
height at Farm Point are presented in Fig. 5. The meas- 
urements of Figs. 4 and 5(a) were preceded by the ad- 
justment of RPL and Farm Point antenna azimuth 
angles for maximum received signal with the RPL 
antenna at the 59-foot level and the Farm Point an- 
tenna at the 112-foot level. For the measurements in 
Fig. 5(b), the Farm Point antenna was reoriented in 
azimuth at the 112-foot level for maximum received 
signal on the Ridge-Farm Point circuit. 
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Fig. 4—Results of some height-gain measurements at 493 mc. In 
these measurements, antenna height was varied at RPL and the 
received signal recorded simultaneously at the Ridge Site and 
Farm Point. The Farm Point antenna was at the 112-foot level. 


It will be noted that each pair of height-gain curves in 
Fig. 4 is characterized by strong similarities but with 
some differences in detail. In particular, the maxima and 
minima were less pronounced on the obstructed circuit 
than on the optical circuit. Among the possible explana- 
tions for this are the following: 

1) It is probable that the greater effective obstruction 
heights associated with the radiation reflected by 
the terrain between RPL and the ridge would re- 
sult in the diffraction losses for this radiation being 
greater than for the unreflected radiation. 

2) Additional regions of stationary phase might be 
expected to result from the nonuniform nature of 
the ridge. Such regions could make substantial con- 
tributions to the total signal level, particularly 
with the RPL antenna at heights for maximum 
destructive interference between the components 
of the direct diffracted signal. 
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Fig. 5—Results of some 493 mc height-gain measurements in which 
antenna height was varied at Farm Point. (a) RPL transmitting. 
(b) Ridge Site transmitting. 


In addition, some of the differences between the curves 
in each pair may be due to differences in the effective 
reflection coefficients resulting from the fact that the 
sites were not located in a common vertical plane. It is 
evident that the effective reflection coefficient of the 
terrain between RPL and the ridge for horizontal 
polarization was greater than for vertical polarization. 

The height-gain curves in Fig. 5(a) were plotted from 
measurements obtained with the RPL antenna at two 
heights which produced maximum signal at Farm 
Point and at one height which gave minimum signal. 
Except for the result obtained using horizontal polariza- 
tion with the RPL antenna at the height for minimum 
signal, the average levels are characterized by a similar 
approximately linear increase with Farm Point antenna 
height. In contrast, average signal levels for the Ridge- 
Farm Point circuit [Fig. 5(b) ] actually decreased slight- 
ly with Farm Point antenna height. From these results 
it has been concluded that the variation of average lev- 
els on the obstructed circuit with Farm Point antenna 
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height was probably due to the ridge itself and not to the 
terrain between the ridge and Farm Point. 

The striking differences in Fig. 5(a) between the shape 
of the curve for horizontal polarization with the RPL 
antenna at the height for minimum signal and the shapes 
of the other curves is probably a consequence of the 
polarization difference of the effective reflection co- 
efficients of the terrain between RPL and the ridge. For 
vertical polarization, the magnitude of this coefficient 
evidently was sufficiently small that the direct diffracted 
signal predominated even when the RPL antenna was 
at a height associated with maximum destructive inter- 
ference between the components of the direct diffracted 
signal. For horizontal polarization, however, it appears 
that the magnitude of this coefficient was sufficiently 
large that other mechanisms, such as propagation via 
secondary regions of stationary phase, made a significant 
contribution to the total signal level when the RPL 
antenna was at this height. 

All of the curves in Fig. 5 are characterized by essen- 
tially similar fine structure. Antenna tilt measurements 
and an analysis utilizing the height separations between 
adjacent maxima and minima suggest that this structure 


was due to interference between the direct radiation and 


radiation reflected by the brow of the hill on which the 
Farm Point Site was located. 

From the results presented in Figs. 4 and 5, it has 
been concluded that the terrain between RPL and the 
ridge made a much greater contribution to the total 
transmission loss for the obstructed circuit than did the 
terrain between the ridge and Farm Point. The corre- 
sponding experimental results indicate that this conclu- 
sion also applies at the other two frequencies. There 
were greater effects due to the terrain to the north of 
the ridge in the results of the measurements utilizing the 
Farrellton Site. 

The height-gain results in Fig. 5 suggest a strong de- 
pendence of the diffraction loss on Farm Point antenna 
height. Such a dependence was evident also in the cor- 
responding results at 173 and 1785 mc. On the other 
hand, it was concluded from the height-gain measure- 
ments on the RPL-Farrellton circuit that the diffraction 
losses for this circuit were not characterized by a marked 
dependence on Farrellton antenna height. Evidence of 
a strong diffraction loss dependence on RPL antenna 
height is not present in any of the height-gain results 
for either of the obstructed circuits with antenna height 
varying at RPL. (Fig. 4, for example.) 


Azgimuth-Gain Measurements 


The azimuth-gain measurements included measure- 
ments in which 1) the azimuth angle of the transmitting 
antenna at RPL was varied and the signal level re- 
corded simultaneously at the Ridge Site and at the site 
in operation to the north of the ridge, and 2) the an- 
tenna at the site in operation to the north of the ridge 
was rotated with RPL and the Ridge Site alternately 
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Fig. 6—Results of some azimuth-gain measurements at 493 mc. For 
each polarization, azimuth angle is measured from the direction 
for maximum signal on the Ridge-Farm Point circuit. RPL an- 
tenna height= 59 feet; Farm Point antenna height~ 112 feet. At 
Farm aes true azimuth of RPL relative to the Ridge Site 
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transmitting. Some results of measurements of the latter 
type are presented in Fig. 6. For these measurements, 
the RPL antenna was at a height which gave maximum 
signal at Farm Point (see Fig. 4). Prior to the measure- 
ments, the antennas at RPL and the Ridge Site were 
oriented for maximum received signal at Farm Point. 
It will be noted that the shapes of the individual curves 
in each pair are quite similar. In addition, the displace- 
ments of the maxima are approximately equal to the 
azimuth of RPL relative to the Ridge Site. Exact equal- 
ity would be expected for a uniform horizontal knife- 
edge if the reflecting surfaces between the terminals and 
the obstruction were uniform and horizontal. 

Some azimuth-gain measurements were made also 
with the RPL antenna at heights which yielded minima 
at the site in operation to the north of the ridge. 
Illustrative results are contained in Fig. 7. Their most 
significant feature is the shift in the azimuth angle for 
maximum signal on horizontal polarization with the 
RPL antenna at the height for minimum signal. This 
shift tends to substantiate the explanation previously 
advanced for the anomalous shape of the corresponding 
height-gain curve in Fig. 5(a). 

The results presented in Figs. 6 and 7 suggest that, on 
493 mc, the direct diffracted signal predominated at 
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Fig. 7—Some results of 493-mc azimuth-gain measurements with the 
RPL antenna at a height for maximum signal at Farm Point 
(59 feet) and at a height for minimum signal (= 72 feet). Farm 
Point antenna height~ 112 feet. 


Farm Point when the RPL antenna was at a height for 
minimum transmission loss. This was not necessarily 
the case, however, when the RPL antenna was at heights 
which resulted in destructive interference between the 
components of the direct diffracted signal. Similar con- 
clusions for the RPL-Farrellton circuit at this frequency 
and for both circuits at 1785 mc were deduced also from 
azimuth-gain measurements. Unfortunately, the situa- 
tion was not so simple at 173 mc. At this frequency, the 
contributions of lateral multipath propagation were 
greater, particularly on vertical polarization. Fig. 8 
illustrates the results obtained. For these measure- 
ments, antenna heights at RPL and Farrellton were 
those which resulted in minimum transmission loss. In 
the results for horizontal polarization, there are only 
minor differences in the shapes of the curves, and the 
displacement of the maxima is approximately equal to 
the true azimuth of RPL relative to the Ridge Site. 
With vertical polarization, there is evidence of consider- 
able distortion in the result for the obstructed circuit. 
This distortion appears to be a flattening relative to the 
free space pattern. Similar 173-mc azimuth-gain meas- 
urements on the RPL-Farm Point circuit at equivalent 
antenna heights also indicated only slight distortion 
with horizontal polarization but considerable distortion 
with vertical polarization. 


Temporal Effects 


Most of the height-gain measurements were repeated 
on different days in order to investigate their reproduci- 
bility and the possibility of seasonal effects. Fig. 9,on page 
208, is a representative example of the diagrams which 
were obtained from the superposition of these sets of ex- 
perimental curves. April 4 was typical of winter condi- 
tions; the Ottawa River was frozen over and the terrain 
was covered with more than a foot of snow. During the 
next few days the snow melted and the ice broke up in 
the river, Shortly afterwards, the buds began to open on 
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Fig. 8—Some results of 173-mc azimuth-gain measurements on the 
RPL-Farrellton circuit. At Farrellton, true azimuth of RPL rela- 
tive to Ridge Site~ —3.2°. 


the deciduous trees and by the end of the first week in 
May the leaves were well advanced toward maturity. 
Undoubtedly some of the differences between the vari- 
ous experimental curves are due to experimental error 
and to atmospheric conditions. However, it seems rea- 
sonable to conclude that there is no evidence of a pro- 
nounced seasonal effect. 

A few short-term fading recordings were made also. 
Some results are given in Fig. 10. For these measure- 
ments, the antenna at RPL was placed alternately at 
a height which produced maximum signal at Farm Point 
and at a height which yielded minimum signal. When 
the RPL antenna was at the former height, fixed at- 
tenuation was inserted between the receiver and its 
transmission line. There was no observable fading in the 
received signal with the RPL antenna at this height. 
The fading with the RPL antenna at the height for 
minimum signal may have been a consequence of the 
antiphased relationship between the direct and ground 
reflected radiation. With this relationship, a small 
change in path phase difference could result in a large 
relative increase in signal level. It is possible also that 
tropospheric scatter may have made a significant con- 
tribution to the total signal level with the RPL antenna 


at this height. 
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Fig. 9—Consolidated plot of 493-mc height-gain results obtained on 
the RPL-Farm Point circuit with antenna height varying at 
RPL and Farm Point antenna height~112 feet. Measurements 
preceded by optimization of azimuth angles with RPL antenna 
height 59 feet and Farm Point antenna height~ 112 feet. 


DIFFRACTION Loss 


Let Zp denote the diffraction loss of the component 
of the radiation which is not reflected by the terrain 
between the terminals of a singly obstructed circuit and 
the obstruction. In this section, estimates of Lp for the 
present circuits are deduced from the experimental re- 
sults and then compared with the values predicted by 
knife-edge theory. 

The condition which must be satisfied in order that 
Lp be essentially independent of antenna characteristics 
is that these characteristics be sufficiently uniform in 
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antenna alternately at a height for maximum signal at Farm 


Point (59 feet) and at a height for minimum signal (~72 feet) ; 
Farm Point antenna height~ 112 feet. 


the directions corresponding to the significant portions 
of the associated regions of stationary phase on the 
obstruction. In general, this condition would be most 
readily satisfied by obstructions which are characterized 
by only one such region. Examples are knife-edges and 
paraboloidal cylinders. 

The widths of the theoretical first Fresnel zones on 
the ridge and the angles which they subtend at the 
terminals have been calculated for the present circuits 
on the assumption that the ridge can be regarded as a 
horizontal knife-edge inclined at an angle of 60° to each 
of the two vertical planes which contained the termi- 
nals. The results of these calculations and the minimum 
antenna beamwidth at each frequency are contained in 
Table I, opposite. It is evident that the angles sub- 
tended by the theoretical first Fresnel zones are con- 
siderably less than the antenna beamwidths. Conse- 
quently, it is concluded that the uniformity condition 
would be satisfied with respect to intensity by the knife- 
edge model if the antennas were beamed toward the cen- 
ters of these zones. It would appear that the uniformity 
condition with respect to polarization was automatically 
satisfied by the essentially linear nature of the antennas. 
The phase characteristics of the antennas were not in- 
vestigated. Nevertheless, it will be assumed that the free 
space equiphase surfaces were spherical over the signifi- 
cant range of directions. 

The rotation measurements have revealed the exist- 
ence of regions of stationary phase which are incon- 
sistent with the knife-edge model. However, except for 
173 mc on vertical polarization, these additional regions 
appear to have made significant contributions to the 
total signal level only when the antenna heights were 
such that there was considerable destructive interference 
between the components of the direct diffracted signal. 
Consequently, although it is possible that the single 
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TABLE I 
THEORETICAL First FRESNEL ZONES ON THE RIDGE AND ANTENNA BEAMWIDTHS 


(iii LiL AA 


Length of First 


Angles Subtended at Terminals by 


Fresnel Zone First Fresnel Zone Minimum 
Ant 
Frequency RPL-Farm Point RPL-Farrellton Heameath 
RPL-_ RPL- (At 3-db 
Farm Point Farrellton At At At At Points) 
RPL Farm Point RPL Farrellton 

173 mc 790 feet 1040 feet 0.64° M25 0.84° 0.60° e 

493 mc 470 feet 615 feet 0.38° 0.75° 0.50° 0.36° ie 

1785 mc 250 feet 325 feet 0.20° 0.39° 0.26° 0.19° 9° 
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regions of stationary phase postulated by the knife-edge 
model were actually composed of a number of regions 
not resolved by the antennas, the rotation measure- 
ments appear to justify reduction of the data in terms of 
the knife-edge model and comparison of the deduced 
experimental diffraction losses with the theoretical 
losses predicted by knife-edge theory, except at 173 mc 
on vertical polarization. 

One of the major problems in the analysis of the ex- 
perimental results was the elimination of the effects of 
the terrain between the terminals and the obstruction. 
In the extension of knife-edge theory which has been 


proposed for the inclusion of these effects,! additional. 


“rays” are postulated, each “ray” being characterized 
by an appropriate diffraction loss and phase shift. How- 
ever, if the differences between the diffraction losses and 
phase shifts associated with the various “rays” are suf- 
ficiently small, the “multi-ray” expression reduces ap- 
proximately to a “one-ray” expression in which the 
effects of the intervening terrain are simply intensifica- 
tion or weakening of the “ray” which would exist in the 
absence of this terrain. If this simplification is valid and 
if the antennas are beamed toward the center of the 
region of stationary phase on the obstruction without 
cross-polarization loss, then 


Ax —Ip+G;+G, (1) 


where A and Lp are as previously defined, and G; and 
G, represent the effective increases in the directivities of 
the antennas toward this region of stationary phase due 
to the terrain between the terminals and the obstruction. 
Eq. (1) is least accurate at antenna heights correspond- 
ing to maximum destructive interference between the 
direct and ground reflected radiation. At.these heights, 
A is most sensitive to the differences between the dif- 
fraction losses and phase shifts associated with the 
various “rays.” 

In the analysis of the experimental data, (1) was as- 
sumed to describe the effects of the intervening terrain 
under the required experimental conditions in the vicin- 
ity of antenna heights corresponding to minimum 
transmission loss. It is unlikely that a less approximate 
extension of knife-edge theory would be applicable at 
other heights because of the increased importance of the 


secondary regions of stationary phase on the ridge. For 
the measurements which were used in the analysis, the 
antennas actually were beamed horizontally instead of 
toward the top of the ridge. However, this probably did 
not cause significant error as the vertical angles involved 
were considerably less than the antenna beamwidths. 
(See Figs. 2 and 3, also Table I; for each circuit, the sum 
of the vertical angles was approximately equal to the 
angle of diffraction.) 

The parameters G; and G, in (1) were identified with 
the appropriate A curves for the optical circuits. Conse- 
quently, some error resulted from the lack of exact co- 
incidence of the reflection areas associated with the 
various obstructed and optical circuits. However, it is 
thought that the ensuing errors in the deduced values of 
Lp are of the order of 1 db at the most. Additional error 
is a likely consequence of the absence of a path-testing 
tower at the Ridge Site. Such a tower would have per- 
mitted the evaluation and elimination of the effects of 
the terrain in the vicinity of this site. Fortunately, the 
shapes of several of the height-gain curves obtained on 
the optical circuits are such that estimates can be de- 
duced of the values of A that would have resulted in the 
absence of reflection by the intervening terrain. (See, for 
example, the appropriate curves in Figs. 4 and 5.) In 
each case, the maximum difference of the deduced value 
from the free space value of A ~0 db is of the order of 2 
db or less. Undoubtedly, some of these differences are 
due partially to experimental error. 

Knife-edge theory predicts that the assumption of 
identical diffraction losses for the direct and reflected 
radiation would yield values of Zp that are too large. 
The error would increase with increased magnitudes of 
the effective reflection coefficients. In principle, it 
should be possible to estimate this error by comparison 
of the depths of the minima in the height-gain curves for 
the optical and obstructed circuits. However, as has 
already been pointed out, lateral multipath was an im- 
portant contributor to the transmission loss for the 
present circuits in the vicinity of the minima. In addi- 
tion, these regions were probably extremely sensitive to 
the differences in the effective reflection coefficients. 
Nevertheless, the height-gain curves have been exam- 
ined on the assumption that the filling of the minima 
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was due entirely to differences in diffraction loss. For 
those results which were susceptible to analysis (Fig. 4, 
for example), it has been concluded that the error in Lp 
would be of the order of 2 db. While this may be a 
reasonable estimate of the error, no attempt has been 
made to correct the values deduced on the basis of (1) 
because of the uncertainty of the relative contributions 
to the minima. 

In the application of (1), the displacement between 
the scales on the left- and right-hand sides of the results 
typified by Fig. 4 was taken to be equal to Lp minus the 
terrain gain G, or G,, as appropriate, for the terrain to 
the north of the ridge. The magnitude of the latter fac- 
tor, which in no case was greater than 1 db, was deduced 
from the height-gain results in which antenna height 
was varied at Farm Point or Farrellton, as appropriate. 
Additional corrections, none of which exceeded 1 db, 
were made to compensate for the differences in the 
height-gain curves obtained on different days. The de- 
duced values of Lp and the values predicted by knife- 
edge theory are contained in Table II. The tabulated ex- 
perimental values of Lp for the RPL-Farm Point circuit 
are respectively the minimum and maximum values for 
the range of antenna heights investigated at the latter 
site. (See discussion in connection with Fig. 5 and the 
comments with respect to the dependence of Lp on an- 
tenna height at RPL and Farrellton.) 


TABLE II 
EXPERIMENTAL AND THEORETICAL DIFFRACTION LOSSES 


Lp(db) 
Fre- 2 
apenry, AY Experimental 
(me) Horizontal | Vertical ae 
Polar- Polar- retica 
ization ization 
173 RPL-Farm Point 27-31 31-36)|* 21 
RPL-Farrellton 25 30 
493 RPL-Farm Point 39-45 41-47 25 
RPL-Farrellton 41 40 
1785 RPL-Farm Point 49->57 48-54 31 
RPL-Farrellton 47 50 


* The significance of these values is conjectural. See discussion in the text. 


The theoretical values of Lp were determined from 
Bullington’s nomogram? with the effective obstruction 
height being taken as 600 feet for the RPL-Farm Point 
circuit and 750 feet for the RPL-Farrellton circuit. (See 
Figs. 2 and 3.) The increases in the theoretical losses 
which result from the assumption of true earth radius 
instead of 4/3 earth radius are approximately 0.2 db for 
the RPL-Farm Point circuit and 0.5 db for the RPL- 


5K. Bullington, “Radio propagation at frequencies above 30 
megacycles,” Proc. IRE, vol. 35, pp. 1122-1136, October, 1947. 
(Fig. 8.) It can be shown that this nomogram may be applied directly 
even though the equivalent knife-edge is not normal to the straight 
line joining the terminals. 
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Farrellton circuit.® It will be noted that at each fre- 
quency knife-edge theory predicts identical diffraction 
losses on both circuits, to the nearest decibel. This is due 
to the fact that the greater distance from the ridge to 
the Farrellton site approximately compensated for the 
larger effective obstruction height. With regard to the 
variation of Lp with antenna height, the derivatives of 
the theoretical knife-edge losses are —0.5 db/100 feet 
at RPL and —1.0 db/100 feet at Farm Point for the 
RPL-Farm Point circuit, and —0.7 db/100 feet at RPL 
and —0.5 db/100 feet at Farrellton for the RPL- 
Farrellton circuit.® 

Consideration was also given to the possibility of 
applying Rice’s theory of diffraction by a parabolic 
cylinder’ to the experimental circuits. In addition to the 
difficulties associated with the boundary conditions, it 
was concluded that the effective radius of curvature at 
the top of the ridge was far greater than can be tolerated 
by the approximate Rice theory. [The restriction 
rh—>>1 associated with Rice’s eq. (2.10) is not satisfied. ] 

Comparison of the deduced experimental diffraction 
losses suggests that there is no evidence of a pronounced 
polarization dependence at either 493 or 1785 me. It is 
unlikely that comparison is valid at 173 mc because of 
the substantial lateral multipath experienced at this 
frequency with vertical polarization. 

From a comparison of the experimental and theoreti- 
cal diffraction losses, it is seen that the predictions of 
knife-edge theory were optimistic. The discrepancies in- 
crease with increasing frequency, ranging from a mini- 
mum of only 4 db at 173 mc to a minimum of 16 db at 
1785 mc. An interesting feature of the results is that, 
in general, at each frequency the differences between the 
experimental values of Zp associated with each circuit 
are relatively small; knife-edge theory predicts that — 
there would be no difference. 


CONCLUSION 


The measurements reported show that simple knife- 
edge theory can result in rather erroneous predictions of 
diffraction losses associated with obstructed circuits. 

Since simple knife-edge theory is of scalar nature, it 
predicts identical diffraction losses for both horizontal 
and vertical polarization. In accordance with this pre- 
diction, the experimental diffraction losses are not char- 
acterized by a marked polarization dependence. 
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6 The asymptotic formula for the Fresnel integral was used in 
these calculations since the theoretical diffraction losses are sufficient- 
ly large that they are given to a good approximation by the use of 
this formula. The asymptotic formula predicts that the derivatives 
of Lp with respect to obstruction and antenna heights are independ- 
ent of frequency. 

7S. O. Rice, “Diffraction of plane radio waves by a parabolic 
cylinder,” Bell Sys. Tech. J., vol. 33, pp. 417-504; March, 1954, 
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A Helical Line Scanner for Beam Steering a Linear Array* 
LOUIS STARK} 


Summary—A new type of antenna beam scanner which uses the 
phasing principle for beam steering a linear array is described. The 
scanner consists of a bank of helical line trombone phase shifters 
which provide a variable delay circuit for each element of the linear 
array. The circuits are ganged mechanically by connecting the slides 


_ of the trombone phase shifters to a lever arm drive. Short helix sec- 


tions concentric with the main lines are used to couple the slide of 


each trombone phase shifter to the main helical lines. 


Measurements of phase characteristics, insertion loss, and vswr 


of the trombone phase shifters have been carried out. Isolation be- 
tween phase shifting circuits has also been measured. Radiation 
studies were conducted using a ten-element scanner in conjunction 
‘with a vertical slot array antenna having ten bays. Patterns showing 
beam scanning out to 52° from broadside are included. 


PHASING PRINCIPLE FOR BEAM STEERING 
Alas DEVICE to be described makes use of the 


phasing principle in antenna theory in order to 

vary the angle of beam formation from a fixed 
linear array. The method is to excite the elements of the 
array with voltages which are phase retarded uniformly 
_ in going from one element to the next. In this way con- 
tributions from the individual elements reinforce on a 
wave front which makes some angle 6 with the array. 
The direction of beam formation is normal to the line of 
reinforcement. The phasing principle is illustrated in 
Fig. 1. It is seen from the figure that the beam angle 6 


BROADSIDE 


DIRECTION OF 
BEAM FORMATION 


WAVEFRONT 


LINEAR ARRAY 


Fig. 1—Phasing principle for beam scanning. 


measured from broadside is given by 


sin @ = D/I 


where 
J=element spacing in inches 
D=incremental delay between elements, in inches of 


free space delay. 


* Manuscript received by the PGAP, August 27, 1956. The re- 
search reported here was supported jointly by the U. S. Army, Navy, 
and Air Force under contract with Mass. Inst. Tech. and performed 
under subcontract with M.I.T. : : 

+ Hughes Aircraft Co., Culver City, Calif. 


PRINCIPLE OF OPERATION OF SCANNER 


The helical line scanner is a ganged arrangement of 
phase shifters providing channels with delays D, 2D, 
3D,+--+, with the number of channels equal to the 
number of radiators. The incremental delay is variable 
By mechanical means. The scanner is placed between a 
suitable power distribution network and the radiating 
elements. 

Each of the phase shifting channels is a helical line 
phase shifter.1 A trombone configuration is employed, 
and the slide is electromagnetically coupled to the heli- 
cal lines by means of coupled helixes as shown schemat- 
ically in Fig. 2. By design, each of the short coupled 


CONCENTRIC COUPLED 
HELIX SECTIONS 


Fig. 2—Schematic drawing of helical line trombone phase 
shifter. 


helixes behaves.as a directional coupler which is unity 
coupled to the main line.1,? Thus, a signal incident from 
terminal 1 is completely transferred to terminal 3, and 
terminals 2 and 4 are not excited. The signal crosses the 
bridge to terminal 7 and is completely transferred to the 
output terminal, no. 5. Similarly, terminals 6 and 8 are 
not excited. This ideal performance is realized very well 
in practice. Ordinarily, only a few helix turns are re- 
quired for the coupled helixes, and performance is good 
over a broad frequency band. 

Displacement of the coupled helix and bridge assem- 
bly varies the phase delay between the input and output 
terminals. If the wavelength on the main helixes is Az 
and if the free space wavelength is A, a displacement of 
A inches changes the electrical length of the channel by 
d=2(No/Aw)A inches of free space delay. The reduced 
wavelength along the main helixes results in a relatively 
rapid delay change with mechanical displacement. The 
helix wavelength is reduced from the free space wave- 
length approximately by the helix wind-up factor. In 
practical designs this factor may range from 10 to 20. 


1L, Stark, “A Helical Line Phase Shifter for Ultra-High Fre- 
quencies,” M.I.T. Lincoln Lab., Tech. Rep. no. 59, February 4, 1954. 

2 R. Kompfner, Bell Telephone Labs., Memo. for File MM-53-150- 
8; March, 1953. 
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Since terminals 4, 2, 6, and 8 are not excited in the 
arrangement of Fig. 2, a complementary channel can be 
provided by bridging across terminals 2 and 6 and using 
terminals 4 and 8 for the second input and output. This 
channel is isolated from the first, and the incremental 
phase delay of one channel increases as it decreases in 
the other. ; 


DESCRIPTION OF AN EXPERIMENTAL TEN- 
ELEMENT SCANNER 


A bank of five complementary phase shifting pairs 
was constructed to form a ten-element scanner. The two 
members of a pair were used to feed radiating elements 
on opposite sides of the center of a ten-element antenna 
array. Each pair of phase shifters was identical with the 
next. A mechanical drive which positioned the slides 
correctly with respect to one another was provided by 
a lever-arm drive pivoted at one end of the bank. The 
scanner is shown in Fig. 3. This design is for the 400—500- 


Fig. 3—Helical line phasing unit assembly. 


mc range. Fig. 4 shows schematically the phase delays 
through the various channels when the incremental de- 
lay is D. Some of the details of construction are shown in 
Figs. 5 and 6. Fig. 5 shows the electrical abstraction of 
a phase-shifting pair. The outer conductors of the coax 


LINEAR ARRAY 


POWER DISTRIBUTION 
NETWORK 


Fig. 4—Delays through phase shifters when incremental 
delay is D. 


Fig. 6—Slide assembly. 
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terminal lines and the shields around the short helixes 
have been removed in this figure. Coupled helixes were 
employed at the input and output terminals to increase 
the breakdown power of the coax-helix transitions.’ The 
coupled helix design allows pressurization of the transi- 
tion region and larger physical spacing between critical 
parts than if a direct connection were made between the 
coax line and inner helix. Fig. 6 shows the construction 
of one of the slides. 


SCANNER TESTS 


Certain tests were performed on the scanner before 
radiation testing was carried out. These tests were con- 
cerned with the following parts of the design and per- 
formance. 


1) Phase delay vs slide displacement characteristic of 
the variable delay channels. 

2) VSWR looking into a channel with matched load 
on output. 

3) Cross coupling between channels. 

4) Channel insertion loss. 


The results of some of these tests are summarized below. 


Phase Characteristic 


A typical measured phase characteristic is plotted in 
Fig. 7. This plot shows the additional path length intro- 
duced into a delay channel when the slide is displaced 
a given amount. The additional path length is given in 
equivalent inches of free space delay. The departure 
from a linear characteristic corresponds to a +5° phase 
error at the measuring frequency, 460 mc. The slope of 
the delay characteristic is 25.5 inches of free space delay 
per inch of displacement. 


Input VSWR of a Delay Channel 


The mismatch looking into a channel was measured 
when this channel and the complementary channel were 
terminated in a matched load. The results are shown in 
Table I. The reflections arise from three sources: the in- 
put terminal, the output terminal, and the bridge junc- 
tion in the slide. The resultant mismatch thus depends 
on both the component reflections and the electrical 
length between them: The component mismatch at one 
of the terminals is shown in Fig. 8. This relatively low 
vswr is the result of some impedance compensation in 
the coax terminal lines. No compensation was used in 
the bridge junctions. 

It is believed that most of the nonlinearity in the 
phase delay characteristic is due to the effect of the 
above reflections. 


3 In Stark, Joc. cit., the breakdown power for a similar type of de- 
sign ae; reported to be 350-kw peak pulse with 6-microsecond pulses. 
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Fig. 7—Phase characteristic. 


TABLE I 
VSWR Looxinc INTO A PHASE-SHIFTING CHANNEL 


VSWR at Bridge Position 
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Frequency, 

mc 0.0 0.25 0.50 0.75 1.0 
inch inch inch inch 

400 1.89 235 1.90 1.45 

420 1.43 ilpges| 1.66 2.10 

440 1.29 1.47 1.63 1.56 

460 1.34 1235 1.42 1532 

480 1.38 1.39 1.20 15 

500 1.64 1.02 1.09 1.67 


VSWR 


FREQUENCY, MCPS 


Fig. 8—Compensated vswr of coax-helix transition at phase- 


shifter terminal. 
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Cross Coupling 


Cross coupling between complementary channels was 
measured over the band 400-500 mc. Matched bolome- 
ters were placed on the output terminal and on the two 
cross-coupled terminals. The relative power levels 
measured are shown in Table II. The numbering of the 

TABLE II 
Cross CoupLinc BETWEEN PUSH AND PULL CHANNELS 
Signal Level Below Input 
Frequency, 
mc At Terminal 1, At Terminal 2, 
db 
400 2552 18.2 
410 30.7 26.2 
420 37.1 30.1 
430 30.0 See) 
440 DS 18.3 
450 30.5 20.1 
460 S56) 23.0 
470 THA 5S) DN 
480 26.7 36.2 
490 28.5 34.0 
500 oes 34.4 


terminals is shown in Fig. 9. Cross coupling is down 
more than 22 db at terminal 1 and more than 18 db at 
terminal 2. Cross coupling between complementary 
channels is attributed primarily to the performance of 
- the coupled helix sections. 

It arises because the energy is not completely trans- 
ferred between the concentric helix sections and because 
the directivity is not infinite. 

Coupling between channels which are not comple- 
mentary but which are physically side by side was also 
measured. The coupling was too small to be measured 
with the apparatus and was less than —50 db. 


TERMINAL NO 2 


TERMINAL WO, | 


Fig. 9—Numbering of terminals for data of Table II. 


Insertion Loss 


The insertion loss through a channel was measured 
with the scanner arm in the broadside position. In this 
position the bridge was 13 inches from both the input 
and the output ends of the main helixes. The loss was 
measured to be 0.5 db. The rate at which the loss 
changed with displacement of the slide was also meas- 
ured, and this result is plotted in Fig. 10. From the slope 
of this curve it can be concluded that substantially all 
the insertion loss is due to attenuation along the helix 
and that radiation loss is small. 

Measurements of attenuation loss along helixes by 
Peter’ show the optimum wire diameter-to-pitch ratio 


4R, W. Peter, J. A. Ruetz, and A. B. Olson, “Attenuation of wire 
helices in dielectric supports,” RCA Review, vol. 13, pp. 558-572; 
December, 1952. 


SLOPE «0.035 db/ INCH 


RELATIVE TRANSMISSION, db 


5.0 7.0 9.0 


BRIDGE DISPLACEMENT, INCHES 


Fig. 10—Rate of loss change with slide displacement. 


to be 0.33. In the design described here this ratio was 
0.28. The conductor loss should be reduced by approxi- 
mately 15 per cent if the optimum ratio were employed. 


RADIATION TESTS 


For antenna radiation pattern studies the scanner was 
connected to an array two elements vertical by ten ele- 
ments horizontal. The elements of the array were 
vertical full-wave slots in a conducting sheet spaced a 
quarter wavelength in front of a reflecting plane. Parti- 
tions were placed between the vertical pairs of slots, 
forming ten cavities. A harness feed connected the 
members of a vertical pair electrically. The horizontal 
spacing of the slots was a half wavelength at 437 mc, and 
the array measured 6 feet by 11.2 feet. 

A branching coaxial line network was used to distrib- 
ute the power from the central input terminal to the ten 
delay channels. A schematic diagram showing the con- 
nections is shown in Fig. 11, opposite. Hybrid-T junc- 
tions were used to obtain the power splits, and by this 
method the power distribution was rendered insensitive 
to the impedance mismatch of the radiating elements or 
of the scanner. The branching network of hybrids pro- 
vided a 3-db drop in excitation at the ends of the array 
for slightly reduced sidelobe levels. 

Some typical radiation patterns are shown in Fig. 12(a) 
and12(b), p.216. Fig. 12(a) shows a runof patterns taken 
at 439 mc with the beam at broadside and at angles of 
13, 27, 43, and 52° off broadside. The displacement of 
the scanner arm corresponding to these angles was 0.60 
inch, 1.20, 1.80, and 2.10 inches measured at the slide 
farthest from the pivot point. The amount of slide 
travel available in this scanner was approximately 
twice as much as was required for these studies. The 
broadening of the beam which occurred as the beam was 
scanned from broadside was in accordance with the re- 
duced projected aperture. 

The above patterns are representative of the results 
obtained over the frequency range 410-480 mc. The 
maximum sidelobe level was less than 14 db below the 
main beam as the beam was scanned from broadside. 
The sidelobe level on broadside was —18.5 db. It was 
observed that the use of hybrid-T junctions in the power 
distribution network had a definite effect in reducing the 
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T- COWNECTOR 


2-BY-10 
FULL- WAVE 
SLOT ARRAY 


HELICAL-LINE 
PHASING UNIT 
(SCANNER) 


1SO- FEED 
CORPORATE STRUCTURE 


HYBRID JUNCTION LEGEND 
SERIES ARM LOAD 


IDE ARMS 


SHUNT ARM 


Fig. 11—Schematic diagram of beam-steerable antenna system. 


sidelobe level when the beam was off broadside. A typi- 
cal comparison of patterns taken with and without the 
effect of the hybrids is shown in Fig. 12(b). For this 
comparison the effect of the hybrids was removed by 
removing the matched loads from the series arms of the 
hybrid 7”s, 

Very little beam steering with changes in frequency 
was noticed. This fact was attributed to a property of 
the helical line delay circuits, which is that the equiva- 
lent length of free space delay in each circuit is inde- 
pendent of frequency. Each phase-shifting circuit be- 
haved like a long-fixed length of coax line. 


CONCLUSION 


A new type of scanner for beam steering a linear array 
has been constructed and tested. Radiation tests demon- 


strating beam scanning out of 52° from broadside were 
carried out. With a feed system giving a —18.5-db side- 
lobe level on broadside, the sidelobe level remained less 
than —15 db in the range to 45° and less than —14 db 
out to 52°. The beam broadening which occurred with 
scan was in accordance with the reduced projected aper- 
ture of the antenna. 

The use of hybrid-T junctions in the feed system was 
found to be an essential factor in maintaining reduced 
sidelobe levels. 
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A Simple Solution to the Problem of the Cylindrical Antenna* 
JESSE G. CHANEY 


Summary—It is recognized that the variational solution of the 
cylindrical antenna does not minimize | Z,, — Z)| , as implied by Storer; 
but yields a minimax of the impedance function. Furthermore, the 
variational solution may be interpreted as simply placing two currents 
in parallel upon the antenna, following which, the solution for the 
magnitudes of the currents and for the driving point impedance may 
be found by elementary circuit analysis. Consequently, a first-order 
solution by the generalized circuit is found by postulating a linearly- 
_ attenuated traveling wave to be superimposed upon the sinusoidal 
standing wave. The results show excellent correlation with other 
theoretical and measured results for both the driving point impedance 
and the current distribution along the antenna. 


INTRODUCTION 
B ROADLY SPEAKING, it may be stated that the 


calculus of variations in the theory of functions 
of real variables serves as a tool for discovering 
certain physical laws. The integrand of a definite inte- 
gral is perturbed by adding to it a parameter times a 
function of the independent variable, subject to the con- 
dition that the function vanishes at both limits of the 
integration. The physical laws are obtained by requiring 
_ the first derivative of the perturbed integral with respect 
to the parameter to vanish. Perhaps the most useful 
laws are those for which the integral is made a minimum 
by the value of the parameter thus obtained. 

Recently, the calculus of variations has been applied 
to many physical problems set up in terms of analytic 
functions of a complex variable. In particular, Storer! 
and Tai? have used it in obtaining a first-order solution 
to the symmetrically driven straight cylindrical an- 
tenna. 

In agreement with Storer,? if 


W = f(z) = U(x, y) + jV(4, 9) 


‘is analytic, and if Wy is a constant, using the Cauchy- 
Riemann equations, it is not difficult to show that 
setting 


fe) 0 
—|W-W,| =0, —|w-W,| =0, 
Ox dy 


is equivalent to setting 


* Manuscript received by the PGAP, April 26, 1956. This paper is 
a revision of Tech. Rep. No. 15, U. S. Naval Postgraduate School, 
Monterey, Calif., January, 1956. 
h Dept. of Ehren U. S. Naval Postgraduate School, Mon- 
terey, Calif. é ’ 

1]. E. Storer, “Variational Solution to the Problem of the Sym- 
metrical Antenna,” Cruft Lab., Harvard Univ., Cambridge, Mass., 
h. Rept. No. 101; 1950. ‘aie 
ar C. T Tai, “A Variational Solution to the Problem of Cylindrical 
Antennas,” Stanford Res. Inst., Stanford, Calif., Tech. Rep. No. 
12; 1950. : ; 

3 Storer, op. cit., pp. 7-8. 


or that the first term of a Taylor’s expansion of the 
function vanishes. 

However, by taking the second derivatives and using 
the Cauchy-Riemann equations in connection with the 
standard tests for maxima and minima of functions of 
two variables,‘ it can be shown that | W—W| is neither 
a maximum nor a minimum when the first derivative 
with respect to z vanishes. That is, if Zo is the true input 
impedance of an antenna, and if Z(e,+je;)=Z(e) is an 
impedance obtained by using I(e, x) =I(«) +en(x) as an 
approximation to the true current distribution function 
I(x), the requirement 


dZ 
@_, 
de 


does not minimize |Z—Z,|, as implied by Storer,5 but 
merely yields a minimax® for the value «=e, found by 
solving the above equation. 


PHYSICAL SIGNIFICANCE OF VARIATIONAL METHOD 


The variational method postulates the generalized 
Kirchoff’s law,” 


qoOkee 
a I(x')G(x’, x)dx’ = ZoI(0)6(x) (1) 
=e 
in which 
1 = the half length of the antenna 
a = the radius of the antenna 
21 
K=— 
nN 
a? etkr(a'—2) 
G(x’, oF | = #)— 
(2, 2) (= 5 r(x’ — x) 


and 6(x) is the Dirac delta impulse function. 
The driving point impedance is then chosen in the 
form 


730 
kI2(0) 


Zo 


[ ference, x)da'dx. (2) 
eit 


It has been shown® that this form, using the approxi- 
mate Green’s function with r= [(«—«’)?+a?]"/? yields 
the same result as that obtained by using the exact 
Green’s function. 


4E. B. Wilson, “Advanced Calculus,” Ginn and Co., New York, 
N. Y., pp. 114-115; 1912. 

5 Storer, op. cit., p. 8. 

6 Wilson, op. cit., p. 115. 

7 Tai, op. cit., p. 4. ; E 

8 C. T. Tai, “A new interpretation of the integral equation formu- 
lation of cylindrical antennas,” IRE TRANs., vol. AP-3; pp. 125-127, 
July, 1955. 
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The current function is perturbed by choosing 
I(«, A) = I(x) + en(x) = Iolf(x) + An(x)]. (3) 


If Ay is a root of dZ)(A)/dA =0, the impedance is 
found to be® 


VirV22 Lie Vie? 


Zo Ao = 610) 4 
me £2(0)V22 — 2f(0)n(0) Vie + 0°(0) Via 4) 
with 
Vu = Lf f poner x)dx'dx 
Vu = £ “Han a)Gla, x)dx’dx 
os Bi ; ! ! ! 
Veg = 5 ff renenee x)dx'dx. 5) 


By simple algebraic manipulation, (4) may be con- 
verted into 
ZZ — 212” 


AON he saregen Pat ea (6) 
Zi + Za. — 2212 


ae 
cae m0 ils f sonore ive 
Zig e ay ip i ONE CLOR DIrEn 
ne aa fe f rerncna, de Coe ekD 


in which it is assumed that neither J,(0) nor J2(0) is zero. 
Eq. (6) is recognized as the solution of the mesh 
equations, 


Z1111(0) + Z1222(0) = Vo 


(8) 
Zo1 (0) + Zo2I2(0) = 
The current distribution becomes 
Vo Z22 —Zi2 
T(z) = I,(x) 
I,(0) ZL 22 at 242" 
Vo Zu-Z 
ee tae: Die Ny oh (9) 
I,(0) ZL ri; 212" 
Suppose that 7.(0) 40 with 7,(0) =0. Let 
TP(a)Zu1 . 7i(x)Zi2 
FRU etn UNL 8 0g 12 = lim ————— - 
z0 = I-91” a0 = Ig4 
Then (6) and (9) become, respectively 
g 
Zin = Zo2 E> Sh (10) 


11 


9 Tai, Rep. No. 12, op. cit., p. 11. 
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and 


I(x) = Vo (= (en £12 


CiZoa — tr? \IR(0) T2(x) - = 1); (11) 


in which Jp;~! is a normalizing factor. 

From (6), (8), and (9), it becomes evident that the 
vanishing of the first term in the Taylor’s expansion of 
the perturbed current simply yields the physical infor- 
mation that two currents have been postulated to exist 
in parallel along the antenna. From (10) and (11), it 
becomes further evident that for singular cases one cur- 
rent becomes the feed current with the other parasitically 
excited. 

This recognition greatly simplifies the variational 
theory of cylindrical antennas. In fact, the statement 
might be ventured that the variational method is not a 
distinct method per se. Any scheme for finding the self 
and mutual impedances of the two postulated current 
distributions may be used. Since | Zin—Zo| is not mini- 
mized, there is no a priori reason for assuming that the 
driving point impedance computed by a method with a 
vanishing first term in its Taylor’s expansion is any 
more accurate than that computed from the parallel 
circuit with the self and mutual impedances being com- 
puted by any recognized method. 


First-ORDER SOLUTION BY THE GENERALIZED CIRCUIT 


For two straight currents, the mutual impedance by 
generalized circuit!’ scheme is given by 


mall io deacon. 
‘tr 


in which the operator Re takes the real part of the prod- 
uct of one current by the complex conjugate of the other 
current. The current functions are normalized at the 
driving point. In case the current functions are different 
and at least one of them is not real, the reciprocity 
theorem requires the use of only the real part of the 
current product. 

It is possible to show that the results obtained with 
the approximate Kernel are the same as those obtained 
when the exact Kernel is used in (12). 

For a relatively low-loss transmission line, Tai! ap- 
proximated the current distribution function by 


I(x) = Io[sin Bl — x) — ja(l — x) cos B(L — x) | 


(12) 


(13) 


from which he chose his pair of currents. The current 
on a lossy line also may be written 


I(x) = Ip[ee* sin B(l — x) — jee sinh a[(l — x) e782) | 
~ Iy[sin B(l — x) — jaev(I — x) e~ #82 ] (14) 
10 J. G. Chaney, “A critical study of the circuit concept,” J. A ppl. 


Phys., vol. 22, pp. 1429-1436; December, 1951. 
uM Tai, Rep. No. 12, op. cit., p. 12. 
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To the writer, this arrangement intuitively seems 
preferable to that of (13). Hence, the current pair is 
chosen as 

I(x) = In sin R(1 — | «| ) 
Io(x) = Ipea(l — | «| )ed#lel, (15) 


The self impedance Z,; is then the well-known zero- 
order solution by the classical induced emf method.” 
For the others, 


30 rt pt y 
a F26 P-KEDG-ED 
| aN Ge a l 
moaeteh (2: oe 
-e/ a|—|2" ene k2 d ad , 
ee “ ) Tie + g°) 
Beef 22a 
212 cant Ne, |) si (i — | x’ | ) 
3? etkrig 
‘cos kx (= + i) dxdx’. (17) 
Ox? 19 


It is not necessary to take the real part in (16) be- 
cause it is a self impedance and the two current func- 
tions are identical. The resulting formulas are listed in 
the Appendix for ka<1, 2a<l. 

The current distribution along the antenna then be- 
comes 


I(x) = 


a Zo —Z 
| sin a0 - cal) 


Z1:1Z02 — Z,,? sin kl 
x 
+ (10 — Z)(1 = =| Jem]. (18) 


ASYMPTOTIC VALUES 


It is of interest to examine the asymptotic values of 
the current and of the impedance. 


21 
For Q = 2in—— « 
a 


60 
Zu =Z2— —j60Qcot kl, Z2—7>—-J 63 2, Zin > Za, 


Sand 
lise Vo sin k(t — | «|) 
ape eS k 1 = x SE pe Res ea Se ee ee 
pe Ce ore 
For k/<1, 


60 
Zu = Ap = Zo — 20(Rl)? rc —2— 2 In 2) 


with Zin =Zi1: 


12S, A. Schelkunoff, “Electromagnetic Waves,” D. Van Nostrand 
Co., New York, N. Y., p. 373; 1947. 


and 


x 


l 


re) 3 *(1 — 
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For kl=nn, n= 17273) 


2 
Zin = Zo. — ga 
Su 
and 
Vo x 
compat 
f1Z22 — £1? uu 1 “ 


— fysin kd — | 2| )|. (19) 


From (19), it is seen that the attenuated traveling 
wave becomes the feed current, and that the standing 
wave becomes parasitically excited. 

Thus asymptotically, this method is equivalent to the 
induced emf method with the removal of the singulari- 
ties. 


RESULTS 


Some values of the driving point impedance computed 
by this method (Tables I and II) were plotted on curves 


TABLE I 

2=10 a=15 
kl 

(ie Lei 
x/2 83.0+741.8 77.5 +942 .3 
Wich 493 +304 360 +7654 
2.6 883 +794 1282 -+71101 
2.9 633 —j7492 2513 —7349 
7 372 —7493 1345 —71447 
3/2 78.3+715.0 92 +32 
5.1 178 + 7104 288 +7287 
526 586 +763 871 +7683 
6.0 340 +7262 1699 — 7130 
Qn 332 +7263 1094 —7904 

TABLE MI 

t kl=7/2 kl=r 

Q 

Zi Zin 
10 83 .0+ 741.8 372 —j493 
12 79 .8+ 742.0 683 —837 
15 77.5 4742.3 1345 —1447 
22 75 .6+742 .3 3797 —73463 


by Tai! (Figs. 1 and 2) for comparison with values com- 
puted by the King-Middleton, the variational, and the 
Schelkunoff methods. They show excellent correlation 
with the other methods, being very close to the results 
of Tai with perhaps about as many of the excursions 
being on the Schelkunoff side of the curves, as on the 
King-Middleton side of the curves. 

For Q=10 and Q=15, the current distribution was 
computed for comparison with the sinusoidal current 
and with the current computed by King and Harrison” 


13 Tai, Rep. No. 12, op. cit., Figs. 4-5. ee 

14 R, King and C. W. Harrison, Jr., “The distribution of current 
along a symmetrical center-driven antenna,” Proc. IRE, vol. 31, 
pp. 548-656; October, 1943. 
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from a first-order solution of the Hallén integral equa- 


tion (Figs. 3 and 4). In view of the known discrepancies 
existing in the King-Harrison method, good correlation 


is indicated. 
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These computed values were also compared with the 
measured values of Barzilai® and of Morita.!® Excellent 
correlation exists between the computed and measured 
results. The curves by Morita showed a comparison of 
the measured values with values computed by the King- 
Middleton first-order solution. The values by the gen- 
eralized circuit first order show better correlation to the 
measured values than do those by the King-Middleton 
first-order method. 


CONCLUSION 


The comparative simplicity of the method of obtain- 
ing the first order generalized circuit solution, together 
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with the excellent correlation with other solutions and 
with published measurements, should justify its intro- 
duction into the already crowded group of antenna 
theories. 

It is believed that the simplicity of the method will 
greatly facilitate the teaching of the theory of the cylin- 
drical antenna. 

The method is being extended to the case of the re- 
ceiving antenna. 

It should be pointed out that this is not a variational 
method, but that the variational method as previously 
presented becomes a special case of this method be- 
cause the current functions were real. 


APPENDIX 


Zu sin? kl = 60cin2kl 30 sin 2k1(2si2kl — si4kl) + 30 cos 2kU(2cin2kl — cin4kl) 
+ 7{60si2k1 — 30 sin 2h1(Q — In 4 + cin4kl — 2cin2kl) + 30 cos 2k1(2si2kl — si4kl)} 


Zi sin kl = 60 cos ki(si4kl — 2si2kl) + 60 sin ki(cin4kl — cin2kl) + j {60 cos ki(in4 — Q9+ 1 — cin4kl + 2cin2kl) 


15 
+ 60 sin Ri(si4ki — si2kl) + a (3 sin kl + 2 sin 2k/ — sin snp 


30 
Zo. = 60cin2kI — 30 + a (sin 2k] — 2si2kl) 


15 
(I)? 


+ 
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15 G. Barzilai, “Experimental determination of the distribution of 
2 en and charge sree cylindrical antennas,” Proc. IRE, vol. 37, 
pp. 825-829; July, 1949. oda 
. 16 T, Morita, “The Measurement of Current and Charge Distribu- 
tions on Cylindrical Antennas,” Cruft. Lab., Harvard Univ., Cam- 
bridge, Mass., Tech, Rep. No. 66, Figs. V-4, V-6; 1949. 
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Investigations with a Model Surface Wave Transmission Line* 
G. GOUBAU} anp C. E. SHARP} 


Summary—This paper summarizes the results of experimental 
studies conducted with a scale model of a surface wave transmission 
line (swt line) for long-distance transmission. The purpose of the 
experiments was to determine the effect of supports, bends, and 
sags on the transmission loss of such a line. Also, the coupling be- 
tween parallel identical lines has been measured and compared 
with the theory. 

Formulas and other information are given which may be applied 
to the design of a swt line for long-distance transmission. 


INTRODUCTION 
Mees FEASIBILITY of employing surface wave 


transmission lines (swt lines) for the transmission 

of wide-band signals over long distances has been 
under study at the Signal Corps Engineering Labora- 
tories for several years. The results indicate that band- 
widths of 100 mc or more are readily obtained with a 
relatively simple and economical installation. The most 
favorable operating frequency range is between 100 
and 300 mc. The practical lower-frequency limit is de- 
termined by the extension of the field, and the upper 
limit by the size of the conductor necessary to obtain 

_ low attenuation. 

The first experiments with a swt line for long-distance 
transmission were made with a line having a segmented 
conductor of the type used for high-tension lines (HH- 
cable manufactured by General Cable Co.) of 1.42-cm 
diameter with a polyethylene coating 0.34-cm thick 
(Fig. 1). This line was installed along a road and sup- 
ported from crossarms of telephone poles by slings of 
nylon cord. The results of these experiments, which have 


~~, 


= POLYETHYLENE 


WAVE CONDUCTOR 


120' ii say ia 


LINE INSTALLATION 


Fig. 1—Experimental surface wave transmission line for 
long-distance transmission. 


* Manuscript received by the PGAP, February 15, 1956; revised 
manuscript received, January 12, 1957. 
{ Signal Corps. Eng. Labs., Fort Monmouth, N. J. 


been published,! demonstrated the feasibility of long- 
distance transmission by swt lines; however, there was 
considerable discrepancy between the loss calculated 
under idealized assumptions, 7.e., without bends, sup-— 
ports, etc. and the actual measured loss. The calculated 
transmission loss at 200 mc, for instance, was 3 db/mile, 
while the measured loss was about 6 db/mile. 

There are several reasons for this’ discrepancy be- 
tween the calculated and the measured loss. First, the 
conductor used had increased resistivity because of its 
segmental structure; the interlocking joints between the 
segments form spirals, and since the current associated 
with the surface wave is axially directed, it is affected 
by the contact resistance of these joints. Second, there 
was an additional loss due to the fact that the line was 
not ideally straight. There were sags between the sus- 
pension points, and bends at the points of support. 
Besides, there were bends in the horizontal plane since 
the line followed the curves of a road. Third, the field 
of the surface wave was intercepted by the line support- 
ing nylon strings, the telephone poles, and the cross- 
arms. The currents induced in such obstacles also result 
in an increase in transmission loss. 

These various loss-increasing factors have been in- 
vestigated individually, not just to explain the dis- 
crepancy between the measured and calculated losses on 
the 2 mile line, but mainly to obtain information leading 
to a better and more efficient construction. Unfortu- 
nately the amount of loss caused by the interlocking 
joints of the line conductor could not be determined — 
since it was found that the loss increased with time. A 
temporary increase in the tension applied to the cable 
would bring the loss back to the lower value. However, 
after a period of several weeks the loss would again in- 
crease appreciably. Apparently the movements of the 
line by the wind affected the joints between the seg- 
ments. Thus, it was concluded that a solid conductor 
should be used in future line constructions. 

Since experiments with the 2-mile line were incon- 
venient for physical reasons, it was decided to continue 
the investigations of the loss factors with a one-tenth 
scale model. The procedures used for the investigations 
with the scale model and the conclusions derived from 
these measurements are given as follows. 


SCALING CONSIDERATIONS 


If, in any electromagnetic transmission system, the 
dielectric constant and permeability of all materials were 


1G. Goubau, C. E. Sharp, and S. W. Attwood, “Investigation of 
a surface wave line for long distance transmission,” IRE TRANS., 
vol. AP-3, pp. 263-267; August, 1952. 
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independent of frequency, and if there were no dielec- 
tric, magnetic, or conductivity losses, then scaling of all 
physical dimensions and the wavelength by the same 
factor would leave the field structure unchanged. The 
loss of such an idealized system would be only radiation 
loss, which is independent of the scaling factor. 

The dielectric constant and permeability of the mate- 
rials used for the construction of a swt line are prac- 
tically independent of the frequency and therefore need 
not be considered in the scaling procedure. Dielectric 
loss factor and conductivity loss follow different scaling 
laws. Fortunately their effect on the field distribution 
around the line is negligibly small. They only cause a 
gradual decrease of the field intensity along the line and 

behave like a uniformly-distributed attenuation which 
is equivalent to a concentrated attenuation inserted at 
one end of a lossless line. Since the line is supported at 
intervals of many wavelengths and these intervals vary 
phasewise at random, the losses caused by each bend, 
sag or support can be considered individually. The total 
_loss in db is therefore the sum of the losses of the in- 
_ dividual discontinuities. The loss at each discontinuity 
is determined only by the cross-sectional field distribu- 
tion and not by the attenuation of the wave along the 
line. Therefore, scaling of dielectric loss factor and con- 
ductivity of the line material is not required. The di- 
- electric losses in the line supporting structures need not 
be considered too, because the polarization currents in 
_ these structures cause primary radiation losses. Thus, a 
geometrically scaled model will suffice to obtain infor- 
mation about the transmission losses which are effected 
_ by bends, sags, and the line supporting structures. 


DESCRIPTION OF THE MODEL LINE 


The model line was constructed for operation in a fre- 
quency band centered at 2000 mc. Assuming the oper- 
ating frequency of the full-size line to be 200 mc, the 
scaling factor with respect to the frequency was there- 
fore 10. A length of about 528 feet simulated a 1-mile 
- section of a full-scale line. The supporting poles were 12 
feet apart corresponding to a scaling by a factor ten of 
the average span of 120 feet between telephone poles. 
- The average height of the model line poles above ground 
was about 5 feet. The height of the line above the ground 
is of no consequence provided it is in the order of wave- 
lengths. The model line poles were 1{ inches square. 
Compared to actual size telephone poles, the thickness 
of the model poles corresponded to a scaling factor of 
about 5 rather than 10. The oversize diameter was cho- 
sen to emphasize the effect of the poles for measuring 
purposes. The crossarms used on the model line poles 
corresponded to the scaling factor 10. The nylon cord 
used for the slings which supported the model line at 
each cross-arm had a diameter of 0.03 inch; with the 
scaling factor 10 applied, this permits a 0.3-inch cord 
for a full-size line having a breaking strength of more 
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than 1000 pounds. Most of the measurements were con- 
ducted with a wire having a conductor diameter of 0.051 
inch (No. 16 B&S) with a coating of brown poly- 
ethylene of 0.034 inch in thickness; the ratio between 
the diameter (do) of the coated conductor and the 
diameter d; of the uncoated conductor was therefore 
2.35. A ratio dy/d; of slightly above two is a good com- 
promise between the copper and dielectric losses of the 
line, and the auxiliary losses due to bends, supports, etc. 


EFFECT OF SUPPORTS 


The simplest and most satisfactory means for sup- 
porting long swt lines from telephone poles are slings of 
nylon cord as shown in Fig. 1. Although this method 
may seem unconventional for supporting transmission 
lines it has considerable merit. The cord is relatively in- 
expensive compared to other kinds of supports con- 
sidered, and in the form of a sling it provides a simple 
way to obtain good alignment of the swt line. The align- 
ment in the horizontal plane is effected by gravity, and 
in the vertical plane by adjusting the length of the nylon 
slings. This type of support also tends to equalize the 
line tension between adjacent spans and therefore the 
line can be kept under a higher average tension, thus re- 
ducing the amount of the undesired sag between sup- 
ports. An additional advantage is that very little ma- 
terial is introduced into the field surrounding the line. 

In order to determine, with the model line, the 
amount of loss due to the nylon cord slings, the support- 
ing cord cross section was varied by placing two or more 
slings in parallel. When the cross-sectional area of a 
support is doubled, then the total polarization current 
is also doubled. The result is that the amplitude of the 
radiated field is increased by a factor two and the 
radiated power by a factor four. Tripling the supporting 
cross section, by placing three slings in parallel, causes 
nine times the radiation loss of a single sling. When two 
slings were added at each point of support the measured 
increase in line loss was approximately 0.2 db. There- 
fore the single sling supports caused § of 0.2 db or about 
0.025-db loss. Since the radiation loss is independent of 
the scaling factor, this same amount of loss is produced 
by the supporting nylon sling of a 1-mile section of a 
correspondingly scaled full-sized line, and of course, this 
loss is negligible. 


EFFECT OF POLES AND CROSSARMS 


Since the poles have a larger cross section than the 
cross arms, and in practice the distance from the latter 
can be made sufficiently large, only the effect of the 
poles was considered. For an approximation, one can 
assume that the loss is proportional to the square of the 
maximum radial electric field strength at the axis of the 
pole. With this simplifying assumption, and the further 
assumption that the diameter of the pole is small com- 
pared to the wavelength, the following formula has been 


224 


derived for the loss L, which is ca sed by a single pole: 


Ly = <5 11,°>(9D) | ab (1) 


ZL. R 


Z =impedance of the line, defined as the quotient 


of the propagated power divided by the square : 


of the current in the conductor. 
H,© = Hankel-function of first order and first kind. 


y= VR —# 
h=2/d,; \, = wavelength of the surface wave. 
k=27/\; \=free space wavelength. 
D =distance of the pole axis from the center of the 
wire. 
F=cross-sectional area of the poles. 
A=empirical constant. The value of A is de- 
pendent on the dielectric constant of the pole. 
The constant A for the model poles (redwood) was 
found to be 1.0105. Since the diameter of the model 
poles was not small compared to the wavelength, A for 
full-scale poles may be somewhat different. 
Fig. 2 shows the measured and calculated loss per 


LOSS PER POLE 
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Fig. 2—Loss due to a pole as a function of the distance D of the line 
from the pole. The cross-sectional area of the model poles is 3 
square inches, the frequency is 2000 mc. Dimensions of the lines: 
Line No. I, Polyethylene insulation, d)=0. 120 inch, d; =0.051 inch 
and Line No. II Polyethylene insulation do=0.130 inch, d;=0.102 
inch. (d)=outer diameter of insulated wire, d;=diameter of wire 
without insulation.) 


pole as a function of the distance (D) from the axis of 
the pole to the center of the conductor (measured in 
wavelengths). The cross-sectional dimensions of the lines 
and the poles used for the model line are given with the 
figure. The curves were calculated with (1) for 
A =1.0X105. The points shown along the curves are the 
measured points. The poles have little effect on the loss 
of line No. I if their distance from the line exceeds 3 X. 
If the data for Line No. I are applied to a full-size line 
with a pole spacing of 120 feet, the loss due to the poles 
will be found to be less than 0.5 db/mile if the poles are 
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at least 2 from the line. Actually the loss would = less 
than supelted from Fig. 2 because the poles of the model 
line were oversize. 


EFrFrEct OF BENDS 


A bend in a swt line radiates like an infinitely long 
V antenna.? The radiation loss from a bend increases 
with the square of the deflection angle, provided this 
angle is small (less than 10°). For larger angles the rela- 
tion between deflection angle and radiation loss is more 
complex. If the bend is not too great, that is if less than 
10 per cent of the energy is radiated, the radiated power 
N, can then be determined with the following formula: 


u2+sin?@ u-+sin@ 
N= 30rl” — n————_ — 2] Watt 
u sin 0 u — sind 
where 
| uw? = h?/k? — cos? 6; 6 = a/2. (2) 


a is the deflection angle at the bend, J is the current in 

the line. The symbols # and & are explained under (1). 

The formula can be greatly simplified if 0=a/2 is very 

small, then 

a2 

1071* 
bv/¢ 


N,= (3) 


Watt (a in radians). 


duv/c is the reduction in phase velocity of the surface 
wave with respect to free space velocity c. 

It is interesting to note that the radiation loss of a 
bend is little dependent on the frequency and mainly 
determined by the ratio between outer diameter do of 
the dielectric coated conductor and the diameter d; of 
the plain conductor. This is illustrated by the curves in 
Fig. 3, opposite, which show the calculated radiation 
loss of bends with small deflection angles as a function of 
do/d:, for \/d;=50 and \/d;=500. The radiation loss 
within the frequency range of 1 to 10 changes by less 
than 10 per cent. 

Fig. 4 shows the radiation loss of bends for two dif- 
ferent lines. The curves were calculated; the points in- 
dicate the measured values. The dimensions of the lines 
are shown in the figure. The agreement between the 
calculated and measured loss is very good for deflection 
angles below 30° that is, as long as the radiation loss is 
considerably less than 1 db. In the case of a line with a 
thinner dielectric layer, the same loss would occur at a 
smaller angle, and the theoretical curve would deviate 
earlier from the measured points. 

The fact that the radiation loss of a bend increases 
approximately with the square of the angle of deflection, 
suggests a method for reducing this loss. If a single bend 
is replaced by two bends of half the deflection angle, 
each of these smaller bends would have one fourth the 
loss of a single bend, and together they would have one 
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Fig. 3— Radiation loss of bends with small deflection angles (a =in 
degrees) for \/d;=50 and \/d;=500 as a function of do/d;. 
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Fig. 4—Measured loss of bends with various deflection angles 
compared with calculated curves. 


half the loss of the single bend. In general, replacing a 


_ single bend by 1 smaller bends each having 1/nth of the 


deflection angle of the original single bend should reduce 


the loss to 1/mth. This reasoning, however, is not con- 


clusive because it implies the assumption that the “sub- 
bends” into which the large bend is divided are radiating 
independently. Actually their radiation fields form a 
resulting radiation field and the corresponding loss may 
differ considerably from that-of independently radiating 
bends. It is obvious that as the distance between the 
sub-bends is continuously decreased, the radiation ap- 
proaches that of a single bend having a deflection angle 
equal to the sum of the deflection angles of the sub- 
bends. 

The experimental study of bends has demonstrated 
that if a single bend is split into two sub-bends, there is 
an optimum spacing between the sub-bends such that 
the loss is a minimum. Fig. 5 is an example of the varia- 
tion in loss when three 30° bends in the line at three con- 
secutive suspension points are divided into pairs of 
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Fig. 5—(a) A 90° turn consisting of three 30° bends each of which 
is split in two sub-bends. The spacing between the 30° bends is 
approximately 24 wavelengths. (b) Loss measured for the 90° 
turn shown in (a) as a function in spacing of the sub-bends. 


equal sub-bends and the spacing between the sub-bends 
is varied. The optimum spacing equals a few wave- 
lengths. The minimum loss is somewhat higher than 
predicted for independently radiating sub-bends. It is 
noted that the optimum spacing is not critical. Spacing 
the sub-bends by two wavelengths is sufficient to ob- 
tain a considerable improvement. A further small im- 
provement is obtained by introducing an additional 
small bend between each pair of sub-bends. The loss of 
the 90° turn is then reduced to about 1 db. This loss is 
actually somewhat smaller than the loss which occurs 
when all the sub-bends are equally distributed around 
the 90° turn. 


EFFECT OF SAGS 


The sags between supports in a swt line also produce 
some radiation loss, mainly due to the bends which are 
introduced at each point of support. Theoretical con- 
siderations show that the loss increases with the square 
of the sag. Assuming that the radiation is due only to 
the bends, the loss can be calculated with (3). But since 
the bends in the line at the supports are not sharp, and 
the line on either side of the bend is curved in the oppo- 
site direction, some compensating effect takes place. 
The measured loss therefore is less than that calculated 
for the bends. Fig. 6 shows an example of the measured 
loss per span as a function of the relative sag (sag s 
divided by the length / of the span). The dimensions of 
the line are given in the figure. The scale on the abscissa 
is proportional to (s/l)?, thus permitting a linear repre- 
sentation of the loss. The dotted curve is calculated for 


4 
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Fig. 6—Loss due to sag as a function of the square of the ratios 
between sag and span. dy/d;=2.35; \/d;=120. 


the bends at the points of support. The measured loss 
is only about 65 to 70 per cent of the calculated loss. 

The measured curve in Fig. 6 can be used for deter- 
mining approximately the loss due to sag of any line. If, 
for instance, the sag of a full-size line is 1 foot at the 
center of a 120-foot span, then (s/l)?=0.7X10-4. As- 
sume d)/d; and \/d; are the same as in Fig. 6; the loss 
due to sag per span is then 1.010-* db. If do/d; is 

_ different, say 1.5, then Fig. 3 may be used to correct 
for the different ratio do/d;. According to Fig. 3, the 
loss in a bend is about twice as large for do/d;=1.5 as it 
is for dp/d;=2.35. The loss per span will thus be about 
2.0X10-? db. The Joss due to sag for 44 spans per mile 
will be about 0.9 db. Since this loss is rather high; the 
ratio dy/d; should be increased to at least 2.0. 

There is an optimum thickness of the dielectric layer 
with respect to the sag. Increased dielectric thickness 
causes increased sag due to increased weight, assuming 
the tension on the line is kept constant. Although the 
thicker dielectric layer tends to reduce the line loss, this 
reduction can be overcompensated by the effect of the 
increased sag. This is illustrated in Fig. 7. In this figure 
a copperweld wire of fixed diameter is considered; the 
dielectric layer is varied while the tension on the line is 
kept constant (40 to 50 per cent of the breaking 
strength). The loss due to the sag is plotted in an arbi- 
trary scale vs the ratio d)/d;. The loss is shown to have 
a broad minimum at do/d;=2. 


CouPpLING BETWEEN Two SWT LINES 


The coupling effects between two swt lines have been 
investigated theoretically by Meyerhoff.* Since this 
theory contains certain simplifying assumptions, it was 
desirable to verify experimentally some of the predic- 
tions of this theory. 

The coupling phenomenon between two swt lines may 
be explained as follows: A system of two parallel lines, 


3A. A. Meyerhoff, “Interaction between surface wave transmis- 
sion lines,” Proc. IRE, vol. 40, pp. 1061-1065; September, 1952. 
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when considered as a unit, is a system which propagates 
two wave modes. If the propagation constants of the 
lines differ and the lines are not very close together the 
two modes are almost identical with the surface waves 
on the individual lines. However, each wave induces a 
small current in the adjacent line, and the reaction 
causes a slight modification of the phase velocity. In the 
case where the propagation constants of the two lines 
are identical, the two modes differ considerably from the 
individual surface wave, independent of the spacing of 
the lines. Each mode has equal currents in both lines. 
The currents are in phase for one mode, and 180° out of 
phase for the other mode. The phase velocities of the 
two modes are different. The difference increases with 
decreasing spacing. If energy is fed to one line only, then 
both modes are excited simultaneously in such a phase 
relation that at the excitation point, the currents add in 
the first line, and cancel in the second line. As the two 
waves propagate along the lines, the phase relation 
varies in accordance with the difference in phase veloci- 
ties. When the mutual phase difference has changed by 
mw, the currents cancel in the first line and add in the 
second line. The total energy is then concentrated 
around the other line. After a further change of the 
mutual phase difference by 7, the energy has returned 
to the first line, and so on. 

In order to obtain complete power transfer the lines 
must have a minimum length which is dependent on the 
spacing of the lines. Assume one line is very long and the 
other is short; either one, say the long line, is excited at 
A (see right-hand side of Fig. 8). If the lines are very far 
apart, the current in the short line is zero. Now, when 
the lines are brought together, a current is induced in 
the second line. For a certain spacing this current has a 
maximum at the termination C. At the same time the 
current at the termination B of the long line is zero. If 
the spacing is further decreased, the current in the short 
line decreases and reaches zero when the coupling is such 
that the energy is transferred twice, first from the long 
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line to the short line, and then back to the long line. 
Further decrease of the spacing brings the energy back 
to the short line, and so on. 

Fig. 8 also demonstrates the energy transfer. The 
curve shows the measured attenuation between the lines 
as a function of their spacing. Zero attenuation means 
the entire energy is transmitted to the short line. The 
measured attenuation between A and C is actually 
somewhat below zero. There is somewhat more energy 
received at the short line than there was received at the 

end of the long line, before the short line was coupled. 
The explanation is that the system of two parallel lines 
has somewhat lower loss than one line alone. 
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Fig. 8—Coupling between two surface wave transmission lines 
with identical wires. 


The calculated spacings for complete power transfer 
are in perfect agreement with the measurements. The 
line sections used in above measurements were carefully 
selected because variations of the dielectric layer of the 
wire by a few per cent will upset the conditions for 
complete energy transfer. 

From the foregoing it is obvious that two identical 
_ long-distance swt lines intended to operate independ- 
ently, should not be strung parallel. However, if the 
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dielectric coats differ by a few per cent the energy ex- 
change becomes very small. The crosstalk between such 
lines still may be considerable if they are only a wave- 
length apart. This crosstalk will be large, however, only 
if the lines transmit signals in the same direction, be- 
cause the coupling is directional. 


CONCLUSION 


SWT lines for long-distance transmission can be sup- 
ported in such a manner that the increase in transmis- 
sion loss effected by poles, crossarms, supporting nylon 
ropes, and the sags between spans is less than a db per 
mile. In order to secure this small loss, the ratio dy/d; 
between outer (conductor plus dielectric) and inner 
diameter (conductor alone) of the line should be at 
least 2. Since the loss of a swt line increases little with 
increasing thickness of the dielectric layer, such a ratio 
is permissible for a low-loss line. If a line for long- 
distance transmission contains curves, the transmission 
loss is further increased by the bends necessary in the 
curves. The amount of additional loss will depend on the 
layout of the line. If do/d; is greater than 2 and radiation 
reducing support methods previously described are used, 
the loss in curves can be kept within reasonable limits, 
say about 1 db for a 90° turn. 

As the result of the information obtained from the 
scale model tests and theoretical studies, it appears en- 
tirely practical to construct a swt line for low-loss long- 
distance transmission with 100 mc or more bandwidth. 
For example, if the conductor diameter is 0.4 inch and 
the outer diameter of the line is 1 inch, and if the di- 
electric layer consists of brown polyethylene, the calcu- 
lated loss is 2.5 db/mile at 100 mc and 4.7 db/mile at 
200 mc. The actual loss will be about 10 per cent higher 
due to imperfections of the conductor surface which, 
are always present. If the loss due to supports is 0.5 
db/mile, then the loss of a straight line section will be 
about 3.3 db/mile at 100 mc and 5.7 db/mile at 200 mc. 
At these frequencies the efficiency of the line is not 
affected by rain. Also dry snow and thin layers of ice 
have little effect on the transmission loss. 
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Antenna-to-Medium Coupling Loss* 
HAROLD STARAS} 


Summary—This paper presents an improved estimate of one of 
the more important systems parameters on a scatter circuit, namely, 
the antenna-to-medium coupling loss. The work presented in this 
paper differs from the earlier work of Booker and de Bettencourt’ in 
that it permits evaluating the coupling loss even when nonconical, 
nonidentical antennas are used at either end of the scatter circuit, and 
it also permits taking the anisotropy of atmospheric turbulence into 
account. An important conclusion is reached that the coupling loss is 
not shared equally by the two antennas (i.e., transmitting and re- 
ceiving) even if they are identical. This is a phenomenon which 
should be fully appreciated by engineers designing a scatter circuit 
or when measuring the coupling loss experimentally. 


INTRODUCTION 


N AWN earlier paper, Booker and de Bettencourt! 
| presented an excellent discussion of the effect of 

very narrow beams on the systems parameters on a 
scatter circuit. In particular, they discussed the phe- 
nomenon of antenna-to-medium coupling loss and 
presented estimates of the loss. However, they re- 
stricted their consideration to identical conical antenna 
beams, a restriction which is removed in this paper. 
Furthermore, in the important region where the antenna 
beamwidths are approximately equal to the medium 
beamwidths, their estimate is obtained by interpolation 
from the case of very broad antenna beams and very 
narrow antenna beams while the estimates presented in 
this paper are evaluated quantitatively in this region. 
In addition, these estimates permit taking into account 
the statistical anisotropy of atmospheric turbulence 
which recent measurements by RCA and other labora- 
tories indicates does exist. Finally, the analysis de- 
veloped in this paper points out clearly that the coupling 
loss on a scatter circuit is not shared equally between the 
transmitting and receiving antennas. An understanding 
of this fact is essential for properly designing a scatter 
circuit and even more so, when conducting experi- 
mental tests to measure the coupling loss. 


Basic FORMULATION 


We begin by referring to the development presented 
by the author? in an earlier paper. This starting point 
can be written for our purposes as follows. 


Poe ha (h — Ho) (fo-f)? 
of ate fe (0? + 9?h?)*(h — Ho/2)? . 


* Manuscript received by the PGAP, May 31, 1956; revised 
manuscript received, January 12, 1957. This work was supported by 
the Signal Corp. under Contract DA-36-039-sc-64555 on a Study and 
Investigation of Tropospheric Scattering. 

+ RCA Labs., Princeton, N. J. 

1H. G. Booker and J. T. de Bettencourt, “Theory of radio trans- 
mission by tropospheric scattering using very narrow beams,” PRoc. 
IRE, vol. 43, pp. 281-290; March, 1955. 

2H. Staras, “Forward scattering of radio waves by anistropic 
turbulence,” Proc. IRE, vol. 43, pp. 1364-1380; October, 1955. 
See, in particular, eqs. (8) and (9) of Sec. V, and the accompanying 
discussion. 


(1) 


where ~ means “proportional to,” (P) is the average 
received scattered power measured relative to free 
space, (fo-f)? represents the antenna power patterns of 
the transmitting and receiving antennas normalized to 
unity along the beam axes, and r is the anistropy 
parameter defined as follows: ‘ 


scale of turbulence in vertical dimension 


scale of turbulence in horizontal dimension 


The other variables x, h, and Hy refer to the coordinate 
system being used and are shown in Fig. 1, which by 


omic BES d/2 d/2 hee 
X 


Fig. 1—Geometry for scatter propagation. 


now has become the symbol of the scatter propagation- 
ist. “a” is the radius of a 4/3 earth. 

Before proceeding to evaluate (1), it may be of in- 
terest to point out that the factor (h—H)/2)- in (1) 
represents the assumption that the intensity of turbu- 
lence varies inversely as the square of the height above 
the earth. This particular assumption leads to a dis- 
tance dependence of d~ for (P). Some investigators be- 
lieve that the distance dependence is more closely ap- 
proximated by d~*. This would imply that the negative 
exponent in the term (h—H)/2) lies between 2 and 3. 
A more detailed discussion of this point is presented by 
Gordon.’ On the other hand, the frequency dependence 
of the average scatter loss, (P), is not affected by the 
term just discussed. A change in the frequency de- 
pendence would imply a change of the scattering cross 
section. The only factor of the scattering cross section 
which appears in (1) is (x?-+r2h2)-2. In order to accom- 
modate a 1/f frequency dependence, for example, this 
factor would be changed to (x*?+r2h?)-/2, This change 
would also increase the distance dependence even with- 
out changing the exponential in the term (h—H)/2). 
Thus, very minor modifications of the integrand in (1) 


3'W. E. Gordon, “Radio scattering in the troposphere,” Proc. 
IRE, vol. 43, pp. 23-28; January, 1955. fit 
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are necessary in order to accommodate different esti- 
mates of the frequency and distance dependence of the 
scatter loss. We will, therefore, attempt to evaluate (1) 
as it stands because the minor modifications necessary 
to explain the distance or frequency dependence will 
hardly change the coupling loss estimate. 

To make progress in evaluating (1), an explicit ana- 
lytical approximation must be made for the antenna 
beam pattern functions. Making a usual assumption 
that the beam shape is approximately Gaussian, then 
each antenna pattern may be considered to have the 
following form in each dimension: 


exp { —0.7(tan B/tan Bo)?} 


where @ is the angle measured from the antenna beam 
axis and fy is the semi-beamwidth, i.¢., the angle between 
the antenna axis and the 3-db down point in the par- 
ticular direction in question. Each antenna may have 
two different Bo, one for each dimension. It was not 
possible to evaluate (1) as it stands even with the 
Gaussian assumption for the antenna patterns; how- 
ever, to a good approximation we can evaluate it as 
follows: First consider antennas which have broad 
beams in the vertical dimension but which may have 
narrow beams in the horizontal dimension and thus 
suffer a coupling loss due to their horizontal beamwidth. 
Then reverse the situation and consider antennas which 
have broad beams in the horizontal dimension but which 
may have narrow beams in the vertical dimension and 
thus suffer a coupling loss due to their vertical beam- 
width. Finally, for antennas with narrow beams in both 
dimensions, the total coupling loss would be approxi- 
mately the db sum of the coupling losses in each di- 
mension. 


EVALUATING THE COUPLING Loss 


With the above approximation, the first integral we 
evaluate (for measuring coupling loss in the horizontal 
_dimension) becomes 


os fe RE 


2 
(x? + H?)? (2) 


—oo 


where H is the parameter rh evaluated approximately 
in the center of the scatter volume. From Fig. 3 and the 
accompanying discussion there; we can estimate that 
H&1.2 rHp. In (2) b: and b; are related to the antenna 
- beamwidths by 6; = (d/2+/0.7) tan Bi, and by = (d/2 0.7) 
tan B:. In the present context (; is the semi-beamwidth 
in the horizontal dimension for one antenna while 2 
would be the semibeamwidth in the horizontal dimen- 
sion for the second antenna. If we define 6-?= b1* +2, 
(2) becomes 
e7 (#/b)2 


0) Deine ie Pt 
P ie Gn) 3 


It is to be noted that in (3), it was explicitly assumed 


(3) 
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that the antenna axes lie in the great circle plane be- 
tween transmitter and receiver. The integral is evalu- 
ated in the Appendix, yielding 


(P) = e@/"T1 — 2(H/0)?][1 — &(/2H/0) | 


2 


where (x) is the probability integral defined by 


(4) 


(0) = = [etna 
0) a eres Gg t. 
ViJd_s 
Eq. (4) has been normalized so that (P)=1 at b= ~, 
1.e., for broad-beam antennas. A graphical presentation 
of (4) is shown on Fig. 2 where the coupling loss is 


14 
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THE ABSGISSA IS @o/. FOR VERTICAL DIMENSION AND 70/1. FOR 
HORIZONTAL DIMENSION 


Fig. 2—Estimate of antenna-to-medium coupling loss. 


plotted as a function of 709/Qo, r being the anistropy 
parameter, 6) being the angular distance between trans- 
mitter and receiver and () being the effective beamwidth 
in the horizontal dimension of the transmitting and re- 
ceiving antennas. The effective beamwidth {2 is defined 
in terms of the individual beamwidths by 


Qo? = Qr? + Qe 


where Q; and Q) are the beamwidths between 3-db points 
of antennas 1 and 2. If one antenna is broad and the 
other narrow in the horizontal dimension, then the 
effective beamwidth is that of the narrow one. If both 
antennas have identical beamwidths, then the effective 
beamwidths is 0.707 that of either antenna. We next 
evaluate the coupling loss due to mismatching antenna 
beamwidths with medium beamwidth in the vertical 
dimension. In this case, we assume that the antennas 
have broadbeams in the horizontal dimension. Inte- 
grating over x in (1) and introducing the Gaussian 
beamwidth assumption with the antenna axes along the 
horizon planes, we obtain 
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(h— Ht) exp 1—-( b jni-a\ ; yt 
(P)~ : : dh. (Sa) 
Hy hh — H/2)? 
Introducing p* 7262 + 0.13694 
ay, Taig Mt + 0.2720? + 0.13694 

Ly FF 

= b? = by? + bo and z = H/2 ’ 


(5a) can be rewritten as 


Pe Vea (Cnc 
Pym mA 


It is to be noted that since these integrals are all normal- 
ized to unity at b= ©, multiplicative factors are always 
neglected. Returning to (5b), it was not possible to 
evaluate it in terms of tabulated functions. However, 
since we are dealing with real positive definite quanti- 
ties, a satisfactory approximation to the denominator of 
the integrand is obtained by the following reasoning. A 
function of the form (z+const)-"Scoe—** for n large. 
This suggests that the denominator of the integrand in 
(5b) can be approximated by an exponential. To fit the 
first two terms of a McLaurin series, the approximation 
is 1/8e-*5*, This approximation was improved by a cor- 

rection factor of the form (1-+1.33z?). A comparison of 
- the given integrand (without the antenna pattern fac- 
tor) and the approximate integrand used is shown in 
Fig. 3.4 It would seem that a quite satisfactory approxi- 
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Fig. 3—A comparison of an approximate integrand 
with the true integrand. 


mation has been made. The integral we then evaluate 
becomes 


(P) ~f ge Horb/2b)°(1 4. 13322) e-8-5edz, (Sc) 
0 


This integral can be evaluated by standard techniques 
yielding 


‘ Both integrands are normalized to have a maximum value of 
unity. 


— Vr/2p(1 + 0.408p? + 0.136p%)er*?[1 — B(p)]} (6) 


where pP60)/0o. - 

This function is also plotted in Fig. 2. As a check, 
(5b) was evaluated for several values of H)/26 by nu- 
merical integration using Simpson’s rule with 18 points. 
The results were within 1 db of those shown in Fig. 2, 
thus confirming the validity of our approximation. 

An important conclusion to be reached from this 
analysis is that the coupling loss is not shared equally 
by the antennas at either end of the scatter circuit. For 
example, utilizing the curves presented in Fig. 2, it can 
be seen that introducing a second narrow-beam antenna 
into a system that has already suffered a coupling loss 
of 6 db, say, in each dimension will not double the 
coupling loss (in db) but will increase it only by about 
1.5 db to 2.0 db in each dimension. It might be of some 
interest to note that when using very high gain pa- 
raboloid antennas, the introduction of a second very high 
gain dish will introduce a coupling loss of about 4.0 db. 
It is to be noted that Booker and de Bettencourt guessed 
that their coupling loss estimate would be too large be- 
cause of their assumptions of isotropic turbulence, a 
scattering parameter constant with height and very 
simplified approximations to the antenna beam shapes. 
A comparison of the Booker-de Bettencourt estimate 
with ours (see Fig. 4) indicates that theirs zs too large. 
For our estimate, it was assumed that the anistropy 
parameter, r, was } and that identical parabolic anten- 
nas were being used at each end of the scatter circuit, Q 
being the beamwidth of each antenna. 

To utilize Fig. 2 most effectively, a satisfactory esti- 
mate for the anisotropy parameter, 7, should be used. In 
an earlier paper already referred to,? it was found that 
r=$ fitted NBS data quite satisfactorily. Using the 
same graphical procedure, a comparison was made be- 
tween Kurihara’s data with theory. It was found that 
7? for the Buffalo-Ithaca path and r&3 for the Weth- 
ersfield-Ithaca path. RCA’s own measurements on a 
153-mile path indicated that an average value for r is #. 
All these estimates are in the same direction, namely, 
r<i (r=1 would represent isotropic turbulence) even 
though the quantitative estimates are fairly variable. At 
the present state of our knowledge, it would seem that 
an estimate of r= should be quite satisfactory. 

It is of interest to compare our theory with some data 


PY. Xurihara, “Trans-horizon microwave propagation over hilly 
terrain,” Proc. IRE, vol. 43, pp. 1362-1368; October, 1955. 


EQUATIONS (4) AND (6) 
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Fig. 4—A comparison of two estimates for coupling loss assum- 
ing identical parabolic antennas at each end. 


that has come to our attention. Recent Lincoln Labora- 
tory data (unpublished) have indicated that on the av- 
erage, no noticeable coupling loss has been observed on a 
618-mile path utilizing a 60-feet dish (freq. =400 mc) at 
the transmitting end while comparing a 60-feet dish 
with a 28-feet dish on the receiving end. For this par- 
ticular situation, our estimates (Fig. 2) would suggest 
that a 2.5-db coupling loss would occur. It is not possible 
to utilize the Booker-de Bettencourt! estimates properly 
for this situation since their estimates are for identical 
antennas at each end of the scatter circuit. However, 
their estimate for the coupling loss between two 60-feet 
antennas on this path is about 13 db. It would appear 
that our estimates are more consistent with Lincoln Labo- 
ratory data and it would furthermore appear that it is 
quite essential to be able to estimate coupling loss for 
dissimilar antennas at either end of the scatter circuit. 
Another source of data of interest in this connection is 
the Bell Laboratory data® obtained on a Labrador path. 
In this experiment, a coupling loss of 2.5 db was meas- 


6K. Bullington, W. J. Inkster, and A. L. Durkee, “Results of 
propagation test at 505 mc and 4,090 mc on beyond-horizon paths,” 
Proc. IRE, vol. 43, pp. 1306-1316; October, 1955. 
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ured (650 hours of data) at 4000 mc when a 28-feet 
transmitting antenna was used and a comparison be- 
tween the received signal on a 5 feet and 28 feet an- 
tenna was made at the receiving end. For this particular 
situation our estimate as obtained from Fig. 2 is 2 db. 
Again it appears that our theory is in substantial agree- 
ment with experiment. 


APPENDIX 


In evaluating an integral of the form 


00 enr% 
Ota peat 
OE rane A: (7) 
the following technique was used. It is to be noted that 
F(q) = : {he tis d 8 
Ti a Tsar eyay i clbganiMgyee 8) 
and that 
1 1 = 
PE a pean ae e-lkIH-ikogp. 
ene LAD TE i een (9) 


Substituting (8) and (9) into (7) yields 


1 is) 
\— fi if einude-sededt 
Whig Be 


Integrating over x by completing the square in the ex- 
ponent of the integrand leaves 


pies i d 
qQ= dp? 


F@ =—- fhe f eBid lb 
d(H?) 2H) _.,, 
ala’ REY Gi gh dee al 
ey a — | oP i4-tdk, (10) 
| d(H?) qHJ, 


Eq. (10) is easily evaluated by completing the square, 
evaluating the integral and then performing the indi- 
cated differentiation recognizing that 


d LT od 
d(H?) 2H dH 
The result for F(q)/F(0) is that given in (10). 
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Precipitation Particle Impact Noise in Aircraft Antennas* 
ROBERT L. TANNER} 


Summary—lIt has been found that precipitation static noise is 
produced in antennas under plastic surfaces, even though these sur- 
faces are covered by conductive films. In this paper, the noise is 
shown to be due to a noise-producing mechanism associated with 
the acquisition of charge by individual precipitation particles upon 
impact in the antenna field region. The existence of such a mecha- 
nism is demonstrated theoretically and substantiated by experiment. 
Production of such noise in antennas is analogous to production of 
shot noise in vacuum tubes. Calculated results indicate that in prac- 
tical antennas this process is capable of generating noise voltages suf- 
ficient to disable lf and mf equipment such as the adf radio compass, 


INTRODUCTION 


HE TYPE of radio interference known as precipi- 
Bi liescec static has long been familiar to aircraft 
radio engineers. Precipitation static is peculiar to 
aircraft in flight, and derives its name from the fact that 
it was first associated with flight through snow storms. 

It is now known that several distinct basic mecha- 
nisms contribute to precipitation static interference. The 
best known of these is associated with what has been 
called “autogenous charging.” Particles of snow, or the 
ice crystals found in high clouds, strike the aircraft, 
_ charging it to a high voltage. The charging continues 
until corona discharges occur at the wing tips and other 
extremities, producing noise which couples into the radio 
receivers. The production of noise by corona discharges 
and the mechanism by which the noise couples into the 
receivers have received considerable attention, and are 
quite well understood.!? Progress has been made in 
alleviating noise from this source, although no com- 
pletely satisfactory method is yet available. The familiar 
wick discharger, incidentally, is directed at overcoming 
this type of noise. 

Another type of precipitation static noise is produced 
in antennas mounted under plastic surfaces such as 
canopies and radomes. The plastic surface is an insulator 
and ice particles striking it charge it to a high voltage 
relative to the surrounding metal portions of the air- 
craft. The charge rapidly builds up until a sparking or 
“streamering” discharge occurs between the plastic sur- 
face and the surrounding metallic structure, with conse- 
quent coupling of noise into the antenna. The solution 
to this type of interference seemed obvious: to coat the 
dielectric surfaces with a film having a conductivity suf- 


* Manuscript received by the PGAP, April 23, 1956; revised 
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ficiently high to prevent the accumulation of static 
charge, but not so high as to act as an rf shield. Con- 
ductive coating materials have been developed, and this 
solution has been tested. Although the application of 
conductive material eliminated the discharges, it failed 
to eliminate the noise. This result made it evident that 
some mechanism besides sparking—one capable of pro- 
ducing noise in spite of the conductive coating—was 
contributing significantly to the observed noise. An 
explanation for the residual noise is proposed in the 
present paper. 


A CouPLING THEOREM 


In discussing the phenomena of precipitation static 
and the manner in which electrical disturbances couple 
into nearby antennas it is convenient to employ a 
coupling theorem developed in the investigation of inter- 
ference from corona discharges.2 The theorem is a 
reciprocity relationship, and can be derived in a manner 
similar to that in which the familiar Lorentz reciprocity 
theorem is derived. The derivation will not be given in 
the present paper. To state the theorem it is necessary 
first to define the two situations illustrated in Fig. 1. 


Region T2 in which charge 
motion, characterized by 


current density ju occurs 


Ue 


Antenna terminals 
region T, 


SITUATION I: Voltage vi is applied to 
terminals T| producing field E! at all 
points of space and in particular in 
the retion T, 


SITUATION II: Motion of charge occurs 
in the region T,- Current density Je 
1s therefore finite in T,- In response 

to the discharge a current i flows in 

the short circuited antenna terminals T, 


Fig. 1—Coupling of electrical effects due to charge motion 
in antenna field. 


For the conditions outlined in these situations the 
coupling theorem states that 


Iw) = 


Be a 
Vio) Jide. (1) 
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Eq. (1) is derived using the Fourier-transformed form 
of Maxwell’s equations, so that, in general, all quantities 


1 ro 
are functions of frequency. The field quantities E! and 


J™ are also, as indicated, functions of the spatial co- 
_ ordinates. To obtain the response as a function of time 
it is necessary to perform the inverse Fourier transform. 
A special condition exists, however, when the electro- 
dynamic situation can be adequately described by the 
quasistatic approximation. Such a condition exists 
when the dimensions of the region in which it is neces- 


sary to define the field E’ are small compared to the 
wavelength of the highest frequencies of interest. 

If the quasistatic approximation is valid, the ratio of 
field to applied voltage is a function of spatial coordi- 
nates only and not of frequency. For such a situation 
(1) can be shown to reduce to 


El 
T() = if Sz (a) Fay thd (2) 


which expresses the terminal current as a function of 
' time directly in terms of the current describing the 
motion of charge, also given as a function of time. The 
presence of conducting material in the field, such as 
conducting films surrounding the antenna, tends to 
render the quasistatic assumption invalid. It can be 
shown in the case of such coatings, however, that if the 
conductivity is not so high as to reduce the sensitivity 
of the antenna no appreciable error results. 


PROCESSES CAPABLE OF PRODUCING NOISE 
IN THE ANTENNA 


It is evident from (1) and (2) that any process involv- 
ing the motion of charge in the field of the antenna will 
produce a short-circuit current at the antenna terminals. 
- Whether such currents have significant components at 
_ frequencies sufficiently high to cause radio interference 
depends upon the time-space distribution of the charges 
_and their motions. In the case of discharge phenomena 
such as corona and over-surface streamering the motion 
_ of the charge is extremely rapid, resulting in short dura- 
_ tion pulses of terminal current having high frequency 
spectral components. Discharge phenomena are ex- 
cluded from consideration in the present investigation, 
however, and we are obliged to search for other 
mechanisms. 

Although previous investigations had established 
quite conclusively that impinging precipitation in gen- 
eral carries negligible charge at the time of impingement, 
the possibility that the motion of charged precipitation 
particles is responsible again suggests itself. Perhaps 
particles charged by impact on the forward portions of 
the aircraft could cause noise upon passing through the 
fields of antennas located further to the rear. An addi- 
tional objection to this possibility, however, is the rela- 
tively low velocity of such particles, which raises some 
doubt as to whether the moving precipitation particles 
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could cause output current pulses sufficiently short to 
contain appreciable spectral components at radio fre- 
quencies. This question can be investigated with the aid 
of (2) which for a single moving point charge can be 
shown to simplify still further to 


ET 
— , 

1) = av (a (3) 
where x’ represents the coordinates of the particle, and 
qg is the particle charge. 

Eq. (3) states that the instantaneous terminal current 
induced by the motion of the particle is proportional to 
the magnitude of the reciprocal field at the point occu- 
pied by the particle. It is evident, therefore, that if the 
field is highly concentrated or has sufficiently rapid 
variation in space, sharp current pulses might be pro- 
duced even without high particle velocities. To obtain a 
quantitative estimate indicating whether conditions 
exist in actual aircraft capable of producing rf noise by 
this mechanism, an investigation was carried out using 
the idealized model illustrated in Fig. 2. The field situa- 


ELECTRIC 
TRAJECTORY FIELD LINES 
OF PARTICLE 
eee ST SS Ae 
= —_— 
\conpucTiNe 
PLATE 


Fig. 2—Field model used in computing terminal current pulses 
due to induction. 


tion illustrated corresponds approximately to that 
which would exist at the metal border surrounding a 
canopy in which an antenna is installed, and has the 
advantage that it can be treated analytically with rela- 
tive ease. 

For the model shown the output current pulse as a 
function of time is given by 


2 1/2 
t+ 4/ * wi (*) (4) 


I(t) = ky ’ 


where k; is a constant which depends upon the velocity 
and charge magnitude of the particle and upon the 
absolute magnitude of the fields at the edge of the plate. 
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The spectrum function of the pulse thus produced can 
also be obtained analytically. It is 


g(w) = <= e~ Wolr) lo], (5) 
je 


Fig. 3 shows the output current pulses produced by 
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Fig, 3—Terminal current pulses produced by induction. 


particles having similar charges but moving in trajec- 
tories spaced different distances from the conducting 
plate. The more peaked pulse corresponds to a spacing, 
Yo, equal to 1 millimeter, while the other corresponds to 
a spacing of 5 millimeters. In both cases a particle 
velocity of 250 meters per second, approximately 500 
miles per hour, is assumed. The corresponding spectral 
functions are shown in Fig. 4. Consideration of the 
spectra of the induction pulses shown, taken in con- 
junction with the fact that in an actual aircraft the field 
areas capable of producing pulses as sharp as those illus- 
trated are extremely small, if any at all exist, indicates 
that this mechanism can probably be neglected as a 
source ‘of noise. It is necessary, therefore, to search 
further for an explanation. 

The situation illustrated by the diagram of Fig. 5, 
together with the result stated in (3), provides the key 
to the production of noise. Eq. (3) states that the in- 
stantaneous value of the terminal current is propor- 
tional to the product of the charge and velocity of the 
particle and of the reciprocal field at the position of the 
particle. In the induction process considered previously 
the field varies with position and the position varies with 
time. The terminal current therefore varies with time 
and thus contains ac components. It is evident that the 
terminal current will also vary with time if the charge 
varies with time, even though the field is constant. This 
is the situation illustrated in Fig. 5. Until the time of 
impact the particle carries no charge and therefore pro- 
duces no reaction in the external circuit. At the time of 
impact, fo, the particle acquires charge instantaneously 
and rebounds with the magnitude of its velocity essen- 
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Fig. 4—Pulse spectra. 
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Fig. 5—Production of terminal current pulses by particle impact 
in antenna field region. 


tially unchanged. This produces a discontinuous or step- 
fronted pulse similar to the one shown in the figure. 

The mechanism of instantaneous charge acquisition 
in the field region appears able to account for the ob- 
served noise. The step-fronted pulses thus produced 
have noise current spectra which at high frequencies 
drop off only inversely with increasing frequency rather 
than exponentially as did the induction pulses. Further- 
more, since the high-frequency components are intro- 
duced by the stepped leading edge, it is obviously not 
necessary to postulate concentrated fields in order to 
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explain the production of high levels of interference at 
relatively high frequencies. 


EXPERIMENTAL CONFIRMATION 


That the step-fronted pulses which have been postu- 
lated do indeed occur and are responsible for the noise 
experienced has been confirmed in an interesting series 
of laboratory experiments. A small draw-down wind 
tunnel, consisting of a 150-gallon vacuum tank fitted 
with an inlet orifice and valve, was placed in a large re- 
frigerated chamber in which ice fogs could be produced 
by introduction of steam followed by seeding with dry 
ice. After air was pumped from the tank, ice fogs were 
seeded in the chamber and the inlet valve opened. A 
stream of ice fog was drawn into the tank at speeds 
approaching Mach 1 and allowed to impinge upon an 
instrumented surface, permitting observation of the 
_ effects by oscillographic and other electronic techniques. 

The crystals occurring in ice fogs are very much 
smaller and more numerous than those found in clouds, 
and consequently have much smaller individual charg- 
ing effects. This fact in addition to the great numerical 
density of the particles—the stream of air flowing into 
the tank is only $ inches in diameter but carries approxi- 
mately 15 million particles per second—creates con- 
siderable difficulty in distinguishing the effects of indi- 
vidual particles. By designing an electrode arrangement 
in which the reciprocal fields are large and highly con- 
centrated in an area of approximately 0.5 square milli- 
meter, but negligibly small over the remainder of the 
surface upon which the airstream impinges, it was 
_ possible to obtain oscillograms of pulses due to indi- 
vidual particles corresponding exactly to those pre- 
dicted. One of these is shown in Fig. 6. The pulse as 
recorded has a finite rate of rise rather than the instan- 
taneous step pictured in Fig. 5. The discrepancy was 
shown, however, to correspond exactly to that which 
would be introduced by unavoidable time constants 
associated with the measuring equipment. The short 
_ duration of the pulses is due to the particles passing out 
of the region of high coupling fields. As explained, the 
region of high fields was made very small in order to 
isolate the effects of single particles. 

The data obtained in the ice-fog tests indicate that 
the charge deposited by the individual particles aver- 
ages approximately 10-" coulombs. It is believed that 
the charging effect per particle of cirrus particles or 
snow will average at least 100 times greater than that 
observed for ice-fog crystals. 

In addition to the oscillographic investigation of in- 
dividual particle effects, measurements were made with 
apparatus devised to permit spectral analysis of the 
noise generated when the stream of ice fog was allowed 
to impinge upon a simulated antenna. The rms noise 
voltage per unit bandwidth was observed to have a dis- 
tribution which fell close to a curve representing a dis- 
tribution proportional to the inverse of the frequency. 
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Fig. 6—Oscillograms of pulses produced by ice-fog particle impact. 


It can be shown that such a spectral distribution would 
be expected if the basic noise-generating mechanism 
were a step-fronted pulse of the type postulated. 


FLIGHt TEST PROGRAM? 


Toconfirm that the mechanism involving instantane- 
ous acquisition of charge in the field region is responsible 
for the noise observed in flight, a series of flight tests was 
carried out at Wright Air Development Center using a 
specially instrumented B-47. The results of the flight 
tests substantiate very satisfactorily the theoretical pre- 
dictions and, in addition, provide much additional infor- 
mation on particle numerical densities and particle 
charge distributions in typical precipitation conditions. 

The flight test data show that in both cirrus clouds 
and in upper levels of cumulo-nimbus clouds particle 
densities ranged between 104 and 105 particles per cubic 
meter. Individual particle charges ranged between 1 and 
40-uy coulomb, the highest values being encountered in 
thunderheads. Maximum particle charging encountered 
in cirrus clouds were under 20-yy coulombs. Average 
charge values were usually between 5 and 10-yy 
coulombs although for one cirrus flight the average 
particle charge was only 3.1-uy coulombs. Distribution 
of particle charges was such that the rms charge values 
exceeded the average values by 15-20 per cent. 

Oscillographic records of the electrical effects pro- 
duced by individual particle impacts reveal information 
concerning the actual course of events during impact. 
In most cases the impulse produced by the particle has 
a leading edge which is not a single large step, as shown 


3 A more detailed account of the results of the flight tests and of 
other aspects of the program is given in the report of Tanner and 
Nanevicz, footnote 5. 
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in Fig. 5; but consists instead of a sequence of small 
steps. This is interpreted as being caused by particle 
shattering. As the particle collides with the surface it 
grinds itself to pieces, and each fragment, as it flies off, 
carries a part of the total charge, thereby producing an 
incremental step. 

Flight test data also revealed some interesting effects 
regarding the dependence of aircraft charging upon 
speed. Earlier investigators found charging rate to be 
approximately proportional to the third power of air- 
craft speed,! while more recent unpublished results ob- 
tained by Wright Air Development Center personnel 
indicate that at higher speeds charging rate increased 
approximately as the sixth power of the speed. It has 
been assumed previously that the particles impinging 
on the aircraft are those in the volume swept by out it, 
so that the number increased in proportion to the speed. 
From this it was inferred that the charge per particle 
must increase as the second to fifth power of the speed. 
The results of the flight tests, while not entirely con- 
clusive, indicate that if particle charge is dependent on 
speed the dependence is certainly much weaker than the 
fifth power, and probably much weaker than the second 
power. Other data were obtained, however, which in- 
dicate that the particle trajectories are strongly de- 
pendent upon aerodynamic factors which vary rapidly 
with speed, so that the number of particles impinging 
can vary as a rather high power of speed. This result is 
corroborated by theoretical work performed at the 
Lewis Laboratory of the NACA by R. J. Brun, R. G. 
Dorsch and their co-workers.’ In connection with the 
problem of aircraft icing, the NACA investigators have 
calculated the trajectories of spherical particles imping- 
ing upon various aerodynamic shapes. While the aero- 
dynamics of ice crystals producing precipitation static 
are much more complicated than the aerodynamics of 
the spherical droplets assumed in the NACA investiga- 
tions, the NACA results apply qualitatively to the 
problem at hand and indicate that the number of parti- 
cles impinging can quite probably increase as rapidly as 
the sixth power of the speed. 


COMPUTATION OF NOISE IN IDEALIZED ANTENNA 


The production of noise in an antenna by the mecha- 
nism just described is similar to the production of shot 
noise in a vacuum tube but considerably more compli- 
cated in detail. Whereas in shot effect the contribution 
of each electron is exactly the same as every other one, 
the contribution to the total “particle impact noise” of 
each precipitation particle depends on its size, its point 
of impact in the antenna field, and its trajectory subse- 
quent to impact. To accurately include all of these 
factors for any practical antenna configuration is beyond 
the capability of available analytical techniques. Valua- 
ble insight and good approximations can be achieved, 


4 NACA, Washington, D. C., Tech. Notes, TN 2476, TN 3147, 
TN 3153, TN 3099, TN 3410, TN 3587. 


however, by calculations using simplified analytical 
models of the antenna together with particle statistics 
determined in the flight tests. Such calculations have 
been made® and give surprisingly good quantitative 
agreement with measured noise values. For a typical 
canopy-mounted adf sense antenna flight test data and 
calculations indicate that particle impact effects pro- 
duce equivalent noise fields in excess of 100 yv/m ina 
2-kc bandwidth at a frequency of 300 kc. 


CONCLUSION 


From the foregoing it is clear that the particle impact 
effect is capable of accounting well for the magnitude 
and spectral characteristics of the noise observed in 
antennas under dielectric surfaces coated with conduc- 
tive films, leaving little doubt this is the responsible 
mechanism. It is important, therefore, to consider the 
implications of this fact for the design of aircraft lf and 
mf antennas. One of the more significant conclusions is 
that there is no simple remedy which will eliminate noise 
in existing antenna designs. So long as appreciable an- 
tenna fields exist in areas exposed to particle impinge- 
ment, noise will be produced. Furthermore, the fact that 
particle impact noise is closely analogous to shot noise— 
of a nonimpulsive nature—precludes the use of special 
tricks of circuitry such as interference blankers which 
have been under development as an antiprécipitation 
static measure. The only measure which appears to hold 
reasonable hope of overcoming particle impact noise is 
to develop antenna designs which do not expose the an- 
tenna field regions to particle impact. 

Another very significant result of the present investi- 
gation is assessment of the relative effectiveness of con- 
ductive coatings in reducing noise. While noise produced | 
by particle impact is generated whether a conductive 
film is present or not, and while this type of noise may be 
high enough to disable radio equipment in conditions of 
heavy precipitation, flight test data as well as subse- 
quent theoretical investigations indicate that impact 
noise is less by about 40 db than noise produced by sur- 
face discharge. Thus, with a properly applied conductive 
coating there should be many circumstances—light pre- 
cipitation or high signal—where equipment should re- 
main operable. Without a conductive coating, the same 
precipitation conditions would almost certainly result in 
complete disablement of the equipment regardless of sig- 
nal level. This suggests that a renewed effort to develop 
satisfactory conductive coatings is warranted. 
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mediate vicinity of a linear antenna of any arbi- 


pas 
The Exact Solution of the Field Intensities from 
a Linear Radiating Source* 
RABINDRA N. GHOSE+ 
HE EVALUATION of field intensities in the im- ——- EH 
| IP(p 552) 


trary length is often of some interest. In this region 
the electric and magnetic field intensities consist of con- 
tributions from the radiation, induction, and electro- 
static fields. An exact solution of these field intensities 


_ which will be particularly suitable for the immediate 


neighborhood of the antenna can be obtained both for 


_ the electric and the magnetic type of linear radiating 


sources with sinusoidal current or voltage distribution 


_ along the axis of the source. 


The generalized solution! of the field intensities for 
any source can be expressed as 


E= — grad V — pA — Curl F 
i= — grad 6 — ef + Curl A (1) 


where A, F, and ® are the magnetic vector potential, the 
electric vector potential, and the magnetic scalar poten- 
tial, respectively. 

For a thin filamentary line source? with sinusoidal 
current distribution along the s or zg direction (Fig. 1), 


k= —(= +6") 4. (2) 


jae \Oz* 


* Received by the PGAP, November 29, 1956. 

+ Ramo-Wooldridge Corp., Los Angeles, Calif. A 

1E. C. Jordan, “Electromagnetic Waves and Radiating Systems, 
Prentice-Hall, Inc., New York, N. Y.; 1950. 5 

2S. A. Schelkunoff, “Advanced Antenna Theory,” John Wiley and 
Sons, Inc., New York, N. Y.; 1952. 
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Fig. 1—Radiation from a linear electric dipole with 
sinusoidal current distribution. 
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Upon expanding sin@(H—s) and sin6(H-+s), one ob- 
tains 


T { H E188 g—1By/ p+ 2-8)” 
panel ae pee ov te, fF 
Oy GW tes ers)" 


H e7iBx/ P+ (2-8)? 
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Making use of the following transformation: 


2—s = psinha, (5) 


the magnetic potential A, can be expressed as 


I () a 
Ae = —— | ~e#arn df e- Bee “day -- eB (H—z) i eBoy 
b b 


8x1 


b b 
= gwen [ eBee tidy + ewan f core da (6) 
a a 


where 
; fed Fs ; g 
@ = sinh7! 5) b = sinh! —, 
p p 
and 
z+H 
@.= sinh! 
p 


Expanding the integrand into binomial series and inter- 
changing the order of integration and summation, one 


finds 
ene FOP) ==) 
p 


k-k! 


A;= | cosh (is —1Bz—k sinh 


471 k=0 
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z 
—2 cos BH cosh (ise k sinh-! =) 
i p 


aH 
Leal (ien-+i8s+ sinh =") (7) 
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Once the magnetic potential is determined, the electric 
and magnetic field components can be determined 
uniquely by the following relations: 


0A, i Vota 
UO es ote 
Op jwe Oz0p 
pe? 
H,=0 E, = —|— + 8) A;,. 
joe \O2? 
He = 0 Es — 0. (8) 


It may be remarked that as no approximation is in- 
volved in determining A,, the field evaluated by (8) will 
be exact. In particular, the numerical computation of 
the field will be much simpler for very small (Gp); that 
is, for a region very close to the antenna, as only a few 
terms of the infinite series need be considered. 

A direct application of this analysis may be helpful to 
determine the region, about a very high power broad- 
cast or television antenna, which may not be safe for 
maintenance personnel because of the harmful effect of 
radiation on human bodies. Another application of the 
analysis will be the determination of self-impedance of 
an antenna or the mutual impedances of parallel an- 
tennas of unequal heights which are often encountered 
in broadcast-antenna phasing problems. 

The problem of determining the field from a slot an- 
tenna with sinusoidal distribution of E along the axis of 
the slot can be solved as a dual problem. If J is replaced 
by £, the electric field at the center of the slot, in the. 
expression of A,, one may obtain the electric vector po- 
tential F; for the slot antenna. The remaining field com- 
ponents can be determined by usual vector operations 
ou Je 
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Experimental Measurement of the Absorption of 
Millimeter Radio Waves Over Extended Ranges* 
C. W. TOLBERT} anp A. W. STRAITON} 


radio waves has been made using path lengths up 
to 61 miles for a wavelength of 8.6 millimeters 
and up to 7.2 miles for a wavelength of 4.3 mm. 
Although some uncertainty exists as to mean meteor- 
ological conditions over the paths, it is possible to draw 
_ some general conclusions as to the water vapor and oxy- 
gen absorption of the atmosphere. 
The approximate values obtained from these meas- 
urements are listed in Table I. 


\ SERIES of measurements of the attenuation of 


TABLE I 
Oxygen Water Vapor 
Wave- Average Tinat ine 
agar Plea Loss in | Breadth | Loss in Breadth 
ee ont get db/km Const. b/ kn Const. 
Cm— (g/m*) Cm— 
4.3 850 0.5 0.09 0.024 0.3 
8.6 850 —- — 0.010 0.28 
8.6 6,500 — — 0.011 0.32 
8.6 14,000 as — 0.0020 0.045 23 
to to WF 
0.0027 0,06 ™ 
INTRODUCTION 


As part of the millimeter propagation research spon- 


sored at The University of Texas jointly by The Office 


of Naval Research and the Signal Corps, measurements 
of the loss of 4.3 and 8.6-mm radio waves have been 
made over a number of paths. These included 3.5, 7.2, 
and 12.1-mile paths in the vicinity of Austin, Texas, a 
50-mile path from Cheyenne Mountain in Colorado to 
‘the plains east of the mountains, and a 61-mile path be- 
‘tween Pikes Peak and Mount Evans in Colorado. 

The measurement technique and the details of each 
set of measurements are described in a series of techni- 
‘cal reports! prepared for the sponsoring agencies. This 
paper is for the purpose of summarizing the results of 
these measurements as they relate to the energy absorp- 
tion of the radio waves by oxygen and water vapor in 
the atmosphere. 


4.3-MILLIMETER MEASUREMENTS 


Paths in the vicinity of Austin, Texas, with lengths of 
3.5 and 7.2 miles were used at two slightly different 
wavelengths of approximately 4.3 mm. The frequencies 
were 69.45 and 70.10 kilomegacycles. Calibration was 


* Manuscript received by the PGAP, February 23, 1956; revised 


manuscript received, May 14, 1956. : 
+ Elec. Eng. Res. Labs., Univ. of Texas, Austin, Texas. 
1 Elec. Eng. Res. Lab., Univ. of Texas, Tech. Reps. Nos. 63, 69, 


70, 73, and 77. 


made by reference to a 1000-foot path and the median 
signal level was determined for a number of 5-minute 
sampling periods. 

Sling psychometer readings were made at each end of 
the path and the moisture content of the air as deter- 
mined from thesesoundings was found to be ingood agree- 
ment with that obtained from the local Weather Bu- 
reau. All of these measurements were made in the ab- 
sence of precipitation. 

The attenuation of these signals below the free space 
values was determined in decibels per mile for each sam- 
ple and the results plotted as a function of the water 
vapor content of the air in Fig. 1. The line shown in this 
figure was plotted to minimize the square of deviation 
of the points from the line. 


rics 
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*— 69 45 Kmc — 7 MILE DATA 
O— 69.45 Kmc — 3 5MILE DATA 
x —70.10 Kmc —3.5 MILE DATA 


ATTENUATION PER MILE IN db RELATIVE TO FREE SPACE 


6 8 10 12 14 16 16 20 
WATER VAPOR CONTENT IN GRAMS PER CU. METER 


Fig. 1—Loss of radio signal per mile vs water-vapor content. 


The accuracy of the data did not justify attempting 
to separate the effects of the two frequencies and all of 
the data were given equal weight. 

The oxygen loss is determined from the y intercept to 
be 0.8 db/mile or 0.5 db/km. This value agrees within 
the limit of accuracy with that predicted by Van Vleck? 
based on laboratory measurements and corresponds to 
the loss that would be predicted for a line-breadth con- 
stant of 0.09 cm7}. 

The water-vapor absorption in grams per cubic meter 
is determined from the slope of the line in Fig. 1 and 
found to be 0.038 db/mile/g/m* or 0.024 db/km/g/m*. 
This value is approximately three times that predicted 
by Van Vleck and corresponds to a line-breadth con- 
Stamu-or 0.0 cm. 


2D. E. Kerr, ed., “Propagation of Short Radio Waves,” M.I.T. 
Radiation Lab. Series, No. 13, sec. 8; 1951. 
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8.6-MM MEASUREMENTS ON 3.5, 7.2, AND 
12.1-MILE PATHS 


Loss measurements were made at 8.6 mm on paths in 
the vicinity of Austin, Texas, with lengths of 3.5, 7.2, 
and 12.1 miles. The attenuation below the free space 
values are plotted in Fig. 2 on a per mile basis as a func- 
tion of the water-vapor content of the atmosphere. The 
losses at 8.6 mm are distinctly reduced over those at 4.3 
mm and the total magnitude of the absorption over 
these paths is small. For this reason, there is a very 
large scatter of the points plotted in Fig. 2. These data 
are inadequate in themselves to separate the oxygen and 
water-vapor effects. The oxygen loss is, however, small 
at 8.6 mm and is known with sufficient accuracy to as- 
sume a zero water-vapor loss value in Fig. 2 of 0.03 db 
as predicted from a line-breadth constant of 0.02 cm™'. 
With this value as the y intercept, a line was so drawn 
that half of the points fell on either side of it. It is felt 
that the accuracy of the data does not justify more 
elaborate methods of drawing this line. 


ATTENUATION IN db/MILE 


Fig. 2—Water-vapor absorption over 3.5, 7.2, and 12.1- 
mile path at Austin, Texas. 


Using the slope of the line in Fig. 2, a water-vapor ab- 
sorption of 0.17 db/mi/g/m’ or 0.10 db/km/g/m is 
found. This value is that which would be predicted from 
line breadth constant of 0.28 cm7}. 


8.6-MM MEASUREMENTS OVER A 50-MILE PATH 


The data obtained at 8.6 mm over a path from Chey- 
enne Mountain to the plains east of Colorado Springs, 
Colo., were treated in the same manner as that dis- 
cussed in the previous section. In this case, however, 
the 8.6-mm signals were compared to a 3.2-cm level with 
a small height-gain correction to this reference level. 

The data for this test are shown in Fig. 3 and the line 
is drawn on the same basis as used in the previous sec- 
tion except that a small correction in the oxygen absorp- 
tion for its change in the density-temperature ratio is 


made in accordance with the approximation proposed 
by Van Vleck.? 
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ATTENUATION IN db 


Fig. 3—Water-vapor absorption over 50-mile path in Colorado. 


The slope of this line is found to be 0.18 db/m/g/m# 
or 0.11 db/km/g/m?. This value is that which would be 
predicted from a line-breadth constant of 0.3 cm for 
water-vapor absorption. 


8.6-MM MEASUREMENTS OVER A 61-MILE 
PEAK-TO-PEAK PATH 


Measurements were made on a path 60.9 miles in 
length between Pikes Peak and Mount Evans at a mean © 
elevation of approximately 14,000 feet. Frequent pre- 
cipitation was encountered on this path, whereas the 
previous measurements were made at time free of rain, 
snow or sleet. The results of the measurements are shown 
in Fig. 4 as a function of the average of the water-vapor 
content measured at the two ends. The loss of 8.6-mm 
radio waves is again compared to a 3.2-cm signal re- 
corded simultaneously. In Fig. 4, there are seven sam- | 
ples of millimeter data with losses of less than 3 db and 
the remaining samples had losses greater than 5 db. The 
lower-loss group contains all of the cases when it was 
reasonably certain that no precipitation was occurring 
on the path and the samples with the higher losses con- 
tain all of the cases where there was evidence of precipi- 
tation. For this reason, the lower-loss points were as- 
sumed to be those associated with no precipitation. 

The oxygen loss was again calculated, using the ap- 
proximate correction for altitude, and a line drawn di- 
viding the lower loss points equally. From the slope of 
this line it is found that water-vapor.absorption is 0.032 
db/mi/g/m* or 0.002 db/km/g/m', corresponding to a 
line-breadth constant of 0.045 cm7}. 

There is some question as to whether the water-vapor 
content as measured on the peaks is representative of 
the whole path. Radiosonde data taken by the Air Force 
at Lowry Field near Denver were analyzed to obtain the 
water-vapor content at 14,000 feet. The data taken at 
0800 were compared to data taken on the peak at 0930 
and they were found to be approximately 25 per cent 
less than humidity measurements taken on the peaks 
on four days when comparison was possible. The mete- 


L 


ATTENUATION oF 0.86 mm RADIO WAVES IN db 
RELATIVE TO 3.2cm SIGNALS 


WATER VAPOR CONTENT IN g/m? 


Fig. 4—8.6-mm attenuation as a function of water-vapor 
content of the atmosphere. 


, orological situation was such that lower water-vapor 
content would be expected over Lowry Field at the ear- 
lier hour, but it is doubtful whether the difference would 
_ be as great as noted. If the water-vapor contents are all 
decreased by 25 per cent of those used in Fig. 4, the 
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slope of the water-vapor line would be increased by 4. 
The resulting absorption by water-vapor would then be 
0.0043 db/mi/g/m* or 0.0027 db/km/g/m*. The asso- 
ciated line breadth constant is 0.06 cm}. It is felt that 
absorption calculated with and without taking into ac- 
count the radiosonde data provides approximate limits 
to the loss values and the average of the extremes is 
taken as the best estimate of the loss values. 

‘There is considerable uncertainty in the results due 
to lack of complete moisture content data along the 
paths and to the interference from other propagation 
factors such as refraction or reflection. Some general 
conclusions can, however, be drawn from the measure- 
ments concerning the absorption of millimeter radio 
waves by oxygen and water-vapor in the atmosphere. 

Only the 4.3-mm measurements provided data on ox- 
ygen absorption. These measurements agreed within the 
limits of accuracy with the losses predicted from labora- 
tory measurements.? 

The water-vapor losses for the measurements at 850- 


foot and 6500-foot elevations for both 4.3 and 8.6 mm 


gave approximately the same line-breadth constant. 
This value was, however, many times larger than recent 
predictions* based on laboratory measurements. 

On the other hand, the measurements made between 
two 14,000-foot peaks gave losses distinctly less than 
those predicted. Becker and Autler® indicated a de- 
pendence on the line-breadth constant on the water- 
vapor content, but of a much smaller value than pre- 
dicted by these measurements. 

It is hoped that additional measurements will resolve 
the discrepancy between the measured and theoretical 
results. 


3G. E. Becker and S, H. Autler, “Water vapor absorption of elec- 
tromagnetic radiation in the centimeter wave length range,” Phys. 
Rev., vol. 70, September, 1946, 
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1924-1925, taught in high schools of Texas, 
1925-1928, and in the Agricultural and Me- 
chanical College of Texas from 1928 until 
April, 1942. He was on duty as a Naval Re- 
serve Officer during World War II and is a 
Commander, U.S.N.R. (Ret.). 

Since October, 1944, he has been profes- 
sor of electronics at U. S. Naval Postgradu- 
ate School, Monterey, Calif., first as an of- 
ficer, and later as a civilian. 
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Professor Chaney is a member of the 
American Physical Society, Pi Mu Epsilon, 
Pi Kappa Delta, and Pi Gamma Mu. 
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William S. Cook was born in Moncton, 
New Brunswick, Canada, on March 16, 1932. 

He attended Mount Allison University 
in Sackville, New 
Brunswick, for four 
years in fulfillment of 
the requirements for 
the B.Sc. degree. 

From 1954 until 
the fall of 1956, he 
was affiliated with 
the Defence Re- 
search’s Radio Phys- 
ics Laboratory in 
Ottawa, Canada. At 
the present time, Mr. 
Cook is affiliated 
with the Federal Electric Corporation, lo- 
cated in Streator, Ill. 
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John H. Crysdale was born on April 12, 
1925, in Toronto, Ontario, Canada. He re- 
ceived the B.A. degree in physics and chemis- 
try from the Univer- 
sity of Toronto in 
1948 and the M.Sc. 


degree in physics 
from McGill Univer- 
sity in 1953. 


With the excep- 
tion of a year spent 
at McGill Univer- 
sity, he has been as- 
sociated with the 
Defence Research 
Board’s Radio Phys- 
ics Laboratory in Ot- 
tawa, Canada, since graduation from the 
University of Toronto. At the present time 
he is a member of the Laboratory’s Upper 
Atmospheric Physics Section. 

Mr. Crysdale is a member of the Cana- 
dian Association of Physicists. 
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John W. B. Day (M’56) was born on 
December 14, 1931, in Vancouver, British 
Columbia. He received the B.A.Sc. degree in 
engineering physics 
from the University 
of British Columbia 
in 1954. 

Since his gradua- 
tion, he has been as- 
sociated with the Ra- 
dio Physics Labora- 
tory of the Defence 
Research Board in 
Ottawa, Canada. At 
the present time, he 
is with the Lower At- 
mosphere Research Group at the Laboratory. 
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Georg Goubau (A’49-SM’56-F’57) was 
born in Munich, Germany, on November 
29, 1906. He received the Dipl. Phys. de- 
gree in 1930, and 
the Dr. Ing. degree 
in 1931, both from 
the Munich Technical 
University. From 
1931.to 1939, he was — 
employed-in research 
and teaching in the 
physics department 
of the same univer- 
sity. During this time 
he was_ principally 
concerned with iono- 
spheric research. He 
established the first German Ionospheric 
Research Station (Herzogstand/Kochel) 
and was in charge of research work there. 

In 1930, Dr. Goubau became professor 
and director of the department of applied 
physics at the Friedrich-Schiller University 
in Jena, Germany. Since 1948, he has beena — 
consultant at the Signal Corps Engineering 
Laboratories, Fort Monmouth, N. J. ; 

Dr. Goubau is senior author of the vol- 
umes on electronics of the “FIAT Review of 
German Science,’”’ organized by the Mili- 
tary Government of Germany. He also is 
the senior author of a book on microwave 
circuits published in 1955 in Germany. In 
1957, Dr. Goubau received the ‘Harry 
Diamond Award” of the IRE. 
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Isadore Katz (SM’56) was born in Phila- 
delphia, Pa., on October 21, 1916. He re- 
ceived the B.S. degree with majors in phys- 
ics and mathematics 
from Temple Uni- | 
versity in 1937. From 
1938 until 1941 he 
worked for the U.S. 
Weather Bureau. In 
1942 he joined the 
staff of M.I.T. Radi- 
ation Laboratory 
where he was engaged 
in research in radio- 
meteorological prob- 
lems. In 1946 he 
joined the Naval Re- 
search Laboratory staff in Washington, 
D.C., where he carried out research on 
microwave propagation and atmospheric 
physics. Since 1952 he has been with the 
Applied Physics Laboratory of the Johns 
Hopkins University where he continues the 
study of various phases of radar propagation 
and atmospheric turbulence. 

He is a member of the American Physical 
Society and a professional member of the 
American Meteorological Society. 


I. Katz 


G 


.> 


* 


Robin I. Primich (S’48-A’50—M’55) was 
born in Johannesburg, South Africa, on June 
17, 1927. He received the B.S. degree in elec- 
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trical engineering from the University of the 
Witwatersrand in 1949, 

From 1950 to 1954 he studied at Impe- 

: rial College, London, 
Eng., on a scholar- 
ship from the Witwa- 
tersrand University 
and with financial as- 
sistance from the De- 
partment of Scientific 
and Industrial Re- 
search, England. 

In 1954 he was 
awarded the diploma 
of Imperial College 
and the Ph.D. degree 
from the University 
of London. Since the beginning of 1955, he 
has been with the Radio Physics Labora- 
tory, Defence Research Board, Ottawa, 


R. I. Primicu 


__ Can., engaged in microwave work. 
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Michael E. Psutka was born on October 


__ 23, 1920, in Kitchener, Ontario, Canada. 


In 1939, he was graduated from Kitch- 
ener-Waterloo Col- 
legiate, Kitchener, 
Ont. He then spent 
four years in the Ca- 
nadian Army, where 
he specialized in ra- 
dar direction finding. 

Resuming his ed- 
ucation, he attended 
the Radio College of 
Canada in Toronto, 
graduating in 1947. 

Since that time, 
Mr. Psutka has been 
a member of the Radio Physics Laboratory 
of the Defence Research Board in Ottawa, 
Canada. 
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Paul E. Robillard was born on September 
26, 1932 in Belle River, Ontario, Canada. 

He was graduated from the Ryerson 
Institute of Tech- 
nology in Toronto, 
Canada, in 1954... 

Upon his gradua- 
tion from the Insti- 
tute of Technology, 
he joined. the Radio 
Physics Laboratory 
of the Defence Re- 
search Board in Ot- 
tawa, Canada. 

At the present 
time, Mr. Robillard 
is associated with the 
Microwave Measurements Group at the Ra- 
dio Physics Laboratory. 
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Chester E. Sharp (A’51) was born on 
June 2, 1906, in Long Branch, N.J. In 1940, 
he joined the Field Radio Section of the 


_ Squier Signal Laboratory at F ort Monmouth, 
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N. J., where he was responsible for the design 
and development of portable frequency- 
modulated field radio equipment. In 1943, 
Mr. Sharp was as- 
signed to the Radio 
Communication Re- 
search Branch, Coles 
Signal Laboratory, 
for investigation work 
on microwave receiv- 
ers and components. 
Later, asa project en- 
gineer with the Radio 
Relay and Microwave 
Section he was respon- 
sible for the develop- 
ment of man-pack 
portable and SHF radio relay equipments. 
He received War Department commenda- 
tion from the Office of the Chief Signal Of- 
ficer in 1946. 

From 1950 to 1955, while with the Radio 
Communication Research Section, he inves- 
tigated improved designs of microwave filters 
and transmission lines. 

Since 1955, Mr. Sharp has been assigned 
to radio wave propagation research and con- 
sultation with the Communication Division 
Research Staff, Signal Corps Engineering 
Laboratories, Fort Monmouth, N. J. 
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Harold Staras (M’50-SM’53) was born 
in New York, N. Y., on December 24, 1922. 
Upon receiving the B.S. degree from City 
College, New York, 
in January, 1944, he 
joined the Langley 
Field Laboratories 
of the National Ad- 
visory Committee for 
Aeronautics as a jun- 
ior physicist where 
he did development 
work on _ electronic 
instruments for use in 
aeronautical research. 

In 1945, he en- 
tered the U. S. Army 
and served as a Counter-Intelligence Corp 
officer. After leaving military service in 
1947, Mr. Staras returned to full time grad- 
uate study, receiving the M.S. degree in 
physics in 1948 from New York University. 
He then joined the Central Radio Propaga- 
tion Laboratory of the National Bureau of 
Standards where he was engaged in the study 
and analysis of tropospheric propagation 
phenomena. During this time, Mr. Staras 
continued this graduate studies in physics at 
the University of Maryland, obtaining the 
Ph.D. degree in 1955. 

In March, 1954, he joined the Advanced 
Development section of RCA, Camden, 
N. J., where he was engaged in the analysis 
of radar performance in presence of target 
scintillation and ground clutter, as well as in 
the performance of communication systems 
ona scatter circuit. 

In 1956, Dr. Staras transferred to the 
Research Laboratories of RCA where, in ad- 
dition to continuing to work on scatter sys- 
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tems, he has become involved in bandwidth 
compression studies of television signals. 

Dr. Staras is a member of the American 
Physical Society, Sigma Xi, and of Com- 
mission II, URSI. 
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Louis Stark (M’54) was born on April 
5, 1926 in Detroit, Mich. He enrolled at 
Massachusetts. Institute of Technology, 
Cambridge, Mass., 
in 1944 and after 
completing one se- 
mester there, he 
served with the U. S. 
Navy for two years. 
He then resumed his 
studies at M.I.T. and 
received the B.S. de- 
gree in 1950 and the 
M.S. degree in 1952, 
both in electrical 
engineering. From 
1950 to 1952, he was 
a research assistant in the M.I.T. Research 
Laboratory of Electronics. 

In 1953, he joined the staff of M.I.T.’s 
Lincoln Laboratory and worked on antennas 
and microwave components for radar sys- 
tems. In June, 1954, he joined the Hughes 
Aircraft Co., Culver City, Calif., where he 
has been doing antenna work. 

Mr. Stark is a member of Sigma Xi and 
Eta Kappa Nu. 
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Robert L. Tanner was born in Idaho 
Falls, Idaho, on December 4, 1921. He at- 
tended Stanford University where he re- 
ceived the B.S. de- 
gree in 1944, the 
M.S. degree in 1947, 
and the Ph.D. degree 
in 1953. 

From 1944 to 
1946, he served in 
the Signal Corps, 
working with pulse 
modulated commu- 
nication equipment 
and with voice cod- 
ing equipment. In 
1947, Dr. Tanner 
joined the electrical engineering faculty of 
the University of Washington and remained 
there until 1950. 

During 1950, he was employed by the 
Boeing Airplane Company where he was en- 
gaged in Antenna research. 

Since January, 1951, Dr. Tanner has 
been with Stanford Research Institute where 
his interests have been largely in the field 
of aircraft antennas and related systems 
problems. 

He presently holds the position of Group 
Head in the Antenna Research Group of the 
Antenna Systems Laboratory. 

Dr. Tanner is a member of Phi Beta 
Kappa, Tau Beta Pi, Sigma Xi, and the 
Scientific Research Society of America. 
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J. R. Wait (SM’56) was born in Ottawa, 
Canada, in January, 1924. He attended Me- 
Gill University for a brief period before en- 
listing in the Canadi- 
an Army in 1942. By 
the end of the war he 
was a foreman in a 
radar workshop at 
Kingston, Ontario. 

He received the 
B.A.Sc. and M.A.Sc. 
degrees in engineer- 
ing physics from the 
University of To- 
ronto in 1948 and 
1949 respectively. At 
this time he was em- 
ployed as a junior research engineer at the 
Hydro Electric Power Commission of On- 
tario where he assisted in the development 
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of an infra-red bolometer. Returning for 


further graduate work to the University of 
Toronto, he obtained the Ph.D. degree in 
electromagnetic theory in 1951. 

From 1949 to 1952, Dr, Wait was associ- 
ated with Newmont Exploration Ltd. of 
Jerome, Ariz. where he conducted theoreti- 
cal and experimental research in electrical 
prospecting. From 1952 to 1955 he was a 
section leader in the Defence Research Tele- 
communications Establishment in Ottawa 
where he was mainly concerned with theo- 
retical problems in radiation. He is now a 
consultant with the Central Radio Propaga- 
tion Laboratory of the National Bureau of 
Standards in Boulder, Colo. 

He is a member of the Research Society 
of America, the Society of Exploration Phys- 


icists, Commissions III and VI of URSI, and . 


the Canadian Association of Physicists. 


William Elwyn Williams was born in~ 


Llangefni, Anglesey, Wales, in July 1931. 


He received B. Se. degree in mathematics — 


from Manchester Uni- 
versity in 1951, the 
M.Sc. degree in 1952 
and Ph.D. degree in 
1954. From 1954-1955 
he was visiting re- 
search associate with 
the Institute of Math- 
ematical Science, 


been concerned with 
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acoustic and electro- 


problems and radio — 


propagation. Since 1955, he has been em- 
ployed by the English Electric Company at 
Luton, Eng. 
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ANDREW CORPORATION, 363 E. 75th St., Chicago 19, Illinois 
Antennas, Antenna Systems, Coaxial Transmission Lines, Design, Development, Production 


ANTLAB, 4950 North High St., Columbus 14, Ohio 
Antenna Pattern Range Measuring Systems 


COLLINS RADIO COMPANY, Cedar Rapids, Iowa 
Antenna Design and Propagation Research Related for Airborne and Ground Communication Systems. 3 


DEVELOPMENTAL ENGINEERING CORP., 1001 Conn. Ave, N.W., Washington, D. C. and Leesburg, Va 
Research, Development, Installation of Antennas and Antenna Equipment for Super Power Stations 


DORNE AND MARGOLIN, INC., 30 Sylvester Street, Westbury, L. I., New York 
Antenna Research and Development—Radiation Pattern Measuring Services. 


THE GABRIEL LABORATORIES, Div. of the Gabriel Co., 135 Crescent Road, Needham Heights 94, Mass. 
Research and Development of Antenna Equipment for Government and Industry. 


HUGHES AIRCRAFT COMPANY, Culver City, California 
Research, Development, Mfr.: Radar, Missiles, Antennas, Radomes, Tubes, Solid State Physics, Computers. 


JANSKY & BAILEY, INC., 1339 Wisconsin Ave. N.W., Washington 7, D.C. 
Radio & Electronic Engineering; Antenna Research & Propagation Measurements; Systems Design & Evaluation 


MARYLAND ELECTRONIC MANUFACTURING CORPORATION, College Park, Md. 
Antenna and System Development and Production for Civil and Military Requirements. 


THE RAMO-WOOLDRIDGE CORPORATION, Los Angeles 45, California 


TRANSCO PRODUCTS, INC., 12210 Nebraska Ave., Los Angeles 25, Calif. 
Antenna Design and Production—Radiation Pattern Measuring Services 


WHEELER LABORATORIES, INC., 122 Cutter Mill Road, Great Neck, New York 
Consulting Services, Research and Development, Microwave Antennas and Waveguide Components. 


The charge for an Institutional Listing is $25.00 per issue or $75.00 for four 
consecutive issues. Application may be made to the Technical Secretary, 
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